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Chapter 1 
Introduction 



1.1 What are algebraic calculations? 

• comparison of numerical to algebraic, or symbolic, computation 



Numerical 


Symbolic 


2/6 -» 0.333333 


2/6 -» 1/3 


2 + 3-5-5 


a; + 2a; — > 3x 


cos(3. 14159) -» -0.999999 


COs(7r) — > — 1 




sin(2a;) — > 2 sin a; cos a; 




4v ^ 2a; 


Jo 1/2 ^r da; ^ - 1438 


f a: , ln|^ 2 -l| 

J z 2 -l ' 2 




a 2 -6 2 ->(a + 6)(a-&) 


numerical 


algebraic 


evaluation 


simplification 



• typical features of algebraic computation 

— computation with arbitrary precision numbers — no rounding 

— computation with symbols and variables (e.g., x,y) 

— computation with functions (e.g., sin, cos, /) 

— manipulation of formulas 

— symbolic operations ( (see chapter2.8) , (see chapter2.10) , (see chapter2.12) , etc.) 

• algebraic computation or symbol manipulation or computer algebra is the field of scientific computation which 
develops, analyzes, implements and uses algebraic algorithms 



1.2 Why do algebraic calculations? 



• in many research fields, one often needs to process very large algebraic expressions (of perhaps hundreds of 
thousands or more terms, e.g., (see chapter5.5.2), computations from the (see chapter5.4.1), (see chapter5.1) 
calculations) or perform long analytical operations 

• compared to a human, a computer does not make errors (assuming the programming is correct! This unfortu- 
nately is not always the case ...) 

• there exist algorithms which cannot easily be performed by a human with pencil and paper (e.g., (see 
chapter2.12) , (see chapter2.13) ) 

• many others: 



— algebraic solutions are usually more compact than a set of numerical solutions; algebraic solution gives 
more direct information about the relationship between the variables than figures 

— algebraic solutions are always exact, numerical solutions will normally be approximations; this can arise 
from rounding and truncation errors, further errors can creep in when the user interpolates data given in 
tabular form 

— computer algebra can save both time and effort in solving a wide range of problems; much larger problems 
can be investigated than by using traditional methods 

— computer algebra reduces the need for tables of functions, series and integrals; symbolic computation 
using computers have highlited many errors in such materials 

— traditional teaching of applied mathematics has to involve much time in teaching techniques of solution; 
computer algebra systems tend to produce solutions quickly and without errors, so they enable more time 
to be devoted to studying the properties of the solution 

— using of computer algebra allows also very effective construction of numerical algorithms and their semi- 
automatic programming by code generation, the effectiveness of the work and reliability of the results 
can be strongly increased 

1.3 History 

• the field originated from the needs of physics researchers 

• the first programs dealing with formulas were written by physicists in order to save them from performing 
long, tedious, error prone calculations by hand 

• 1955 - first programs for formula derivation 

• 1965 - first general purpose computer systems working with algebraic expressions 

• 1975 - new research field with its own conferences, journals, etc. 

• 1990 - general spreading of (see chapter3) into almost all branches of science 

1.3.1 Chronology of computer algebra systems 

• table by Brian Evans jevans@eedsp.gatech.edu,;, 



System name 


Year 


Related systems 


e-mail address, tel. 


ALPAK 


1964 


ALTRAN 


(Bell Labs) 


ALTRAN 


1968 




(Bell Labs) 


FORMULA (Algol) . 








FORMAC 




FORMAC (PL/I) 


(IBM) 


FORMAC (PL/I) 






(IBM) 


MATHLAB (DECUS) 


1968 


MACSYMA 


(DEC) 


CAMAL 








REDUCE 


1968 




reduce-netlibOrand . org 


MACSYMA 


1970 


Symbolics Macsyma, 
VAXIMA, DOE-Macsyma, 
ALJABR, ParaMacs 


(See Below) 


SchoonShip 


1971 




archive.umich.edu (FTP) 


muMath 


1979 


Derive 




VAXIMA 


1980 




(312) 972-7250 


SMP 


1982 


Mathematica 


NOT ON MARKET 


Symbolics MACSYMA 


1983 


macs 


yma-service@symbolics . com 


DOE-Macsyma 


1984 


ALJABR 


gcookSllnl . gov 


Maple 


1985 




wmsiSdaisy . Waterloo . edu 



MathCAD 


1985 (?) Mathcad 


Powermath 


1985 


REDUCE/PC 


1986 


Derive 


1988 


Mathematica 


1988 


Theorist 


1988 


PARI 


1988(?) 


FORM 


1989 


MACSYMA/PC 


1989 


ALJABR 


1991 


Mathcad 


1991 


SymbMath 


1991 


Axiom 


1991 


ParaMacs 


1991 


SIMATH 


1992 



1-800-MATHCAD 
NOT ON MARKET 
reduce-netlibSrand . org 
(Soft Warehouse Inc.) 
inf oOwri . com 
(415) 543-2252 (Prescience Corp) 
ftp to math.ucla.edu 
form@can.nl 
macsyma-serviceSsymbolics . com 
al j abrSf pr . com 
1-800-MATHCAD 
chenOdeakin .oz.au 
(708) 971-2337 
lphOparadigm. com 
marcSmath . uni-sb . de 



Chapter 2 

Algorithms for algebraic computation 



2.1 Algebraic structures 

• basic requirements 

— precise representation of algebraic structures 

— precise arithmetic with algebraic structures 

— other analytical operations with these structures (e.g., (see chapter2.8) , (see chapter2.10) ) 

2.1.1 Number domains 

• {Z,+,—, xjintegers with operations of addition, subtraction, multiplication; see the examples of (see chapter4.6.1) 

• {Z m , + m , — m , x m }integers under modular arithmetic where m is a positive integer 

• {Q, -I-, — , x , /jrational numbers with operations of addition, subtraction, multiplication and division; see the 
examples of (see chapter4.6.1) 

• {{x = a + ib, a,b € Z}, +, — , xjGaussian integers, i.e., complex numbers which have integer real and imaginary 
parts; see the (see chapter4.6.1) examples 

• {<9(a), +, — , x , /},p(a) = Oalgebraic extension field; the algebraic number ais defined by the polynomial pwith 
integer coefficients which can be represented precisely; the algebraic number ais a root of the polynomial p, 
e.g., \/7is represented by the polynomial a 2 — 7and the algebraic number fccould be defined by the polynomial 
36 2 — 56+1; see the examples of (see chapter4.6.1) 

• {R], +, — , x,/}floating point numbers with a precision of ndecimal digits where nis an arbitrary positive 
integer — it can be lOOor 1000; see the examples of (see chapter4.6.1) 

2.1.2 Algebraic expression domains 

• ring of polynomials 0[x\, ...,x n ]in nvariables with operations of addition, subtraction, multiplication and 
exponentiation by a nonnegative integer; see the examples of (see chapter4.6.2) ; polynomial coefficients can 
be numbers from a (see chapter2.1.1) 

• power series, other kinds of series 

• rational functions K (x i, ...,a;„) (extension of polynomials by the operation of division) with operations of 
addition, subtraction, multiplication, division and exponentiation by an integer; see the examples of (see 
chapter4.6.3) 

• extension of rational functions by radicals (rational exponents), with operations of addition, subtraction, 
multiplication, division and exponentiation by a rational number 



• algebraic functions {yi,pi(X,yi, ...,y m ) = 0,i = 1, ...,m}implicitly defined by polynomials with integer coef- 
ficients pjwhich depend on algebraic functions y«and variables from X = {x\,X2, ...,a;„}, e.g., the algebraic 
function ^defined by the polynomial x 2 + y 2 — 1 

• elementary transcendental functions exp,ln, extension of rational functions by elementary transcendental func- 
tions; if we have only one variable a;and an expression contains exp x, In x, we can denote y = exp x, z = In a;and 
work with a rational function of x,y,z — the extension is given only by rules like x = exp In a;, In a; 2 = 2 In a;, 
etc. 

• transcendental functions: e.g., sin, cos, rmerf; extension of rational functions by transcendental functions 

• matrix rings 

• differential fields (K, 1 ) 

• finite groups 

• a user can use algebraic expressions from an arbitrary domain for the most part; the program will decide which 
domain the expression belongs in and use an appropriate algorithm 

2.2 Representation of algebraic structures 

• to work on a computer with algebraic structures, we need to represent them by some sort of data structures 

• representation is very important because often the effectiveness of an algorithm will depend on the represen- 
tation that is used 

2.2.1 Representation of integers 

• integers are standardly represented by one computer word of nbits; in such a representation, the size of an 
integer ais limited by a < 2™ _1 ; for example, 

if n = 16, a < 32768, 

if n = 32, a < 2147483648and 

if n = 64, a < 9223372036854775808 

• however, we also need to represent arbitrarily big integers; e.g., see (see chapter2.5.3) 

• one possibility to to represent a big integer aby an array Awhere all elements of the array are limited by 
A t < 2™- 1 and 



a = J2A i (2^Y 



i=0 



n-\ 



• thus, the big integer is expressed in a number system with base 2 

• of course, memory is finite and so we can store only a number of a limited size, however using this representation 
we can work with quite large integers, e.g., in 1 kB of memory we can store an integer of size 10 2588 

2.2.2 Representation of polynomials 

• representation by a string of characters is not advantageous as it is not a dynamical structure and algorithm 
implementation would be difficult 



Prefix representation 

• polynomials are represented by lists using e.g., the prefix operators PLUS for addition, DIFFERENCE for 
subtraction (or unary MINUS), TIMES for multiplication and EXPT for exponentiation; the first element of 
a list is the prefix operator and the rest of the elements are its arguments 

• the polynomial 

4a; 3 + 2x - 5 

• is represented by the list (PLUS (TIMES 4 (EXPT X 3) ) (TIMES 2 X) (MINUS 5) ) 

• this representation can be used for arbitrary algebraic expressions 

• algorithms using this representation are not particularly fast 

Dense representation 

• a polynomial in one variable x 



^a(f)a 



is represented by the array of all its coefficients a(i),i = 0, ..., n 

• in this representation, one needs for the polynomial 

2* 1000 - 1 

to store lOOlcoefHcients even though only 2nonzero coefficients are really necessary; for polynomials in more 
variables, the situation is even worse 

• the time complexity of the algorithm for adding two polynomials of degree nin this representation is of the order 
0(n); polynomials can have high degree and only a few terms so that most of the operations are unnecessary 
(such as the addition of two zero coefficients) 

• most polynomials which we encounter in real life are sparse, that is, most of their coefficients are zero 

Sparse representation 

• a polynomial in one variable x 



J2 a ^ xi 



i=0 

is represented by pairs of corresponding exponents and coefficients (i, a(i))for each term of the polynomial, 
i.e., all coefficients a(i)'m this representation are nonzero 

• in this representation, the polynomial 

2a ,iooo _ j 

would be represented by the list ( (1000 2) (0 -1) ) 

• the time complexity of the algorithm for adding two polynomials with nterms is 0(n) 

• to increase the effectiveness of algorithms using this sparse representation, a rule is usually applied which 
specifies the order of pairs in the list based on the exponents, e.g., in the example above, the lists are sorted 
in descending order of exponents 



Recursive representation 

• recursive representation is sparse representation with exponent ordering used to represent multivariate poly- 
nomials 

• for recursive representation, an ordering of variables has to be chosen, e.g., alphabetically 

• the variable which is chosen first in the ordering is called the main variable of the polynomial 

• coefficients of the powers of the main variable will be polynomials in the other variables 

• for simplicity, we will consider a polynomial in 2variables x,y, with variable ordering x > yso that the main 
variable is x 

n m n m 

'Y^^2 i a{i,3)x l y :i =~Y^c{i)x l , where c(i) = J^ a(i, j)y J 

• such a polynomial is represented by a list of pairs (i, c(i))where each coefficient c(i)is a polynomial in yand is 
represented by the list of pairs (j,a(i,j)) 

• to distinquish a polynomial in x, yfrom a polynomial in a, c, we replace each pair of (exponent, coefficient) by 
the triple (variable, exponent, coefficient) 



the polynomial 



\y 2 +8y)-x 2 (7y + 3)-5xy + 4 



is then represented by the list ((x 3 ((y 2 1) (y 1 8))) (x 2 ((y 1 7) (y 3))) (x 1 ((y 1 5))) 
(x 4)) 

Recursive representation in Macsyma 

• Macsyma Canonical Rational Expression (CRE): 



x 7 +4x 3 -2x + 11 

(4zy 2 + 5)x 3 - Ix 

Recursive representation in Reduce 

• Reduce standard form of a polynomial 



((x 7 1 3 4 1 -2 11) . 1) 

((x 3 (y 2 (z 1 4) 5) 1 -7) . 1) 



(((* . 7) . 1) ((x . 3) . 4) ((x . 1) . -2) . 11) 
(((* . 3) ((y . 2)((z . 1) . 4)) . 5) ((* . 1) . -7)) 



x 7 + 4x 3 - 2x + 11 - 
(4zy 2 + 5)x 3 -7x — 

2.2.3 Representation of Expressions 

• general prefix representation (external form vs internal form): 

sin(aa; + b) — > (sin (+ (x a x) b)) 

v v ' 

• internal representation in terms of linked lists: 



sin 










X 
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2.3 Arithmetic 

• by arithmetic, we understand this to mean operations of addition, subtraction, multiplication and division 

2.3.1 Numeric vs symbolic arithmetic 

numeric arithmetic symbolic arithmetic 

• normally inexact • always exact 

• usually produces some sort of answer (may • frequently cannot produce an answer (prob- 
suffer loss of accuracy, fail to converge, etc.) lem is intractable, expression swell, etc.) 

• arithmetic relatively cheap • arithmetic often expensive 

• zeroes usually not given special treatment • zeroes can greatly simplify calculations by 
except to avoid computational singularities making some of them unnecessary 

• zero detection: number sufficiently small? • zero detection: expression equivalent to 

zero? 



• 



maximize size of pivots to promote stability 



• minimize complexity of 

pivots to stunt expression growth 



2.3.2 Arithmetic of integers 

• algorithms for addition, subtraction and multiplication of integers are analogous to procedures used by humans 
for hand calculations; the difference is that with hand calculations, we calculate in the decimal number system 
while these algorithms use the number system with base 2 n_1 

• there exists more effective algorithms for multiplication of integers 

• the standard procedure for division of integers requires an initial estimate — there exist algorithms that do this 
estimate with reasonable precision 

2.3.3 Arithmetic of polynomials 

• polynomial in one variable in (see chapter2.2.2) 

• (see chapter2.3.3) 

• algorithm for the addition of two polynomials 

PLUSP0L(a, b) := 

if a = ( ) then return b 

else if b = ( ) then return a 
else 

ea := first first a 

eb := first first b 

ca := second first a 7,a=((ea ca) ...) 

cb := second first b °/b=((eb cb) ...) 

return ( 

if ea > eb then 

cons(first a, pluspol(rest a, b)) 
else if ea < eb then 

cons(first b, pluspol(rest b, a)) 
else if ca + cb = then 
pluspol(rest a, rest b) 
else cons(list (ea, ca + cb) , pluspol(rest a, rest b)) 
fi fi fi) 
fi fi 



• for example, adding the two polynomials 4a; 2 — 3a;and 5x + 7would be performed (after converting them into 
sparse representation) by PLUSP0L[( (2 4) (1 -3)),((1 5) (0 7))] => ((2 4) (1 2) (0 7)) which is the 
sparse representation of the sum 4x 2 + 2a; + 7 

• algorithm for the multiplication of two polynomials a, b; if a = a\ + G^and b = b\ + 62 then ab = a\b\ + a\bi + 
a<ib\ + a 2 &2 

TIMESPOL(a, b) := 

if a=()orb=() then return ( ) 
else 

ea := first first a 

eb := first first b 

ca := second first a '/,a=((ea ca) ...) 

cb := second first b °/b=((eb cb) ...) 

return(cons(list(ea + eb, ca * cb) , 

pluspol(timespol(list first a, rest b) , 
timespol(rest a, b)) 
) 
) 
fi 

• for example, multiplying the two polynomials Ax — 3a 5a; + 2would be performed (after converting them into 
sparse representation) by TIMESP0L[((1 4) (0 -3)),((1 5) (0 2))] => ((2 20) (1 -7) (0 -6)) which 
is the sparse representation of their product 20a; 2 — 7a; — 6 

Lists 

• a list is an ordered set of elements (a 1 b . . . ) 

• an element of a list can be a number, an identifier or another list, e.g., (a (2)) 

• empty list ( ) 

• selection operations 

— first returns the first element of a list: first [(a (2))] => a 

— second returns the second element of a list: second [(a (2))] => (2) 

— rest returns everything but the first element of a list: rest [(a (2))] => ((2)) 

• constructive operations 

— cons adds an element to the beginning of a list: cons[l, (a b)] => (1 a b) 

— list builds a list from its arguments: list [a, b] => (a b) 

2.4 Simplification 

• most operations in computer algebra are some form of simplification 

• here is a question on what representation to use for algebraic expressions: which form is simpler? 

(a + b) 2 <S> a 2 + 2ab + b 2 

• Sis a canonical simplification operator if and only if 

Vt, S(t) ~ t 
Vti,t 2 h ~ t 2 =► S(h) = S(t 2 ) 



• the canonical simplification operator gives us an unique form for equivalent formulas: it defines a unique 
representative for each class of equivalent algebraic expressions 

• most often non-canonical simplification for general algebraic expressions is used in computer algebra programs — 
typically, a combination of canonical simplification and pattern matching 

2.4.1 Canonical simplification on algebraic domains 

• for polynomials, it is necessary to define a term ordering — often the ordering of polynomial variables is sufficient 
and terms are ordered according to the variables present and their degrees 

• for rational functions 

— canonical simplification of numerator and denominator 

— dividing numerator and denominator by their greatest common divisor 

— ensuring that the leading coefficient (coefficient of the first term of a polynomial in a given ordering) 
satisfies some condition (e.g., making sure that the leading coefficient of the denominator of a rational 
function is positive) 

• for rational functions extended by rational exponents (radicals) 

— unnested radicals (a radical cannot be inside another radical): there exist theories for simplifying these 
forms 

— nested radicals: the theories are very complicated — it is possible to make use of the same methods as for 
algebraic functions 

• algebraic functions: simplifications are performed "modulo" the Groebner basis of the system of polynomials 
defining the algebraic functions 

• elementary transcendental functions and transcendental functions: there exist structural theorems — these are 
very complicated 



2.4.2 Complexity of expressions 

• which is simpler? 

(x - l) 10 = x w - 10a; 9 + 45a; 8 - 120a; 7 + 210a; 6 - 252a; 5 + 210a; 4 - 120a; 3 + 45a; 2 - lOx + 1 



a; 10 -l 
• how about here? 



(x - l)(x + l)(a; 4 - x 3 + x 2 - x + l)(a; 4 + x 3 + x 2 + x + 1) 



Form 


Expression 


complexity 


factored sum 

expanded 

factored 


(a; + l) lu -(a;-l) lu 

20a; 9 + 240a; 7 + 504a; 5 + 240a; 3 + 20a; 

4a;(a; 4 + 10a; 2 + 5) (5a; 4 + 10a; 2 + 1) 


13 
24 
25 



one classification scheme: 



Symbol 


Classification 


Complexity 


R 

F 
B 

E 


rational numbers (integers are consid- 
ered to have a denominator of 1) 

floating point numbers 

extended precision floating point 
numbers (bigfloats) 

expressions that are not simple 
numbers 


sum of the absolute value of the nu- 
merator plus the denominator 

absolute value of the number 

absolute value of the number 

size of the expression (the number of 
operators and atomic operands) 



2.5 Greatest common divisor 

• abbreviated by GCD 

• very important algorithm — it is used by many other algorithms 

2.5.1 GCD of integers 

• for two integers a,6for which a > b, we define their quotient quot(a,6)and remainder rem(a, 6)as two integers 
for which 

a = quot(a, b)b + rem(a, b) 

where < rem(a, b) < b 

• if we denote the greatest common divisor of a, foas gcd(a, 6)and if rem(a, 6)is nonzero, then 
gcd(a, b) = gcd(6, rem(a, b)) 

• gcd(a, 6)divides both aand b; it also divides quot(a,6)6so it also has to divide rem(a, b) 

• the Euclidean algorithm for calculating the greatest common divisor is based on this identity 

• Euclidean algorithm for calculating the GCD of two integers 



gcd := GCDI(a, b) := 
[a, b are integers 
algorithm used: 

rem(a, b) - remainder after dividing integer a by integer b] 
1. if a < b then 





r := a 






a := b 






b := r 






fi 




2. 


while b != 


do 




r := rem(a, 


b) 




a := b 






b := r 






od 




3. 


return a 





7, ! = means not equal 



2.5.2 GCD of polynomials with rational coefficients 

• here we will consider polynomials from Q[a;]in one variable arwith rational coefficients 

• the degree of polynomial a(a;)will be denoted by deg a(x) 

• for two polynomials a(x),b(x)from Q[a;]for which deg a(x) > deg b(x), we define their quotient quot(a(x), 6(a;))and 
remainder rem(a(a;),6(a;))as polynomials from Q[a;]which satisfy 

a(x) = quot(a(x),b(x))b(x) +rem(a(x),b(x)) 

where < deg(rem(a(a;),6(a;))) < deg(fo(a;)) 

• as was the case for integers, the following identity for the greatest common divisor (GCD) of polynomials holds 
gcd(a(x),b(x)) = gcd(b(x),rem(a(x),b(x))) 

• Euclidean algorithm for calculating GCD of two polynomials with rational coefficients 



gcd := GCDPQ(a(x), b(x)) := 
[a, b are polynomials in Q [x] 
algorithms used: 

deg(a) - degree of polynomial a 

rem(a, b) - remainder after dividing polynomial a by polynomial b] 
1. if deg(a) < deg(b) then 



r := a 




a := b 




b := r 




fi 




while b != 


do 


r := rem(a, 


b) 


a := b 




b := r 




do 




return a 





• this algorithm requires many greatest common divisor calculations of two integers when performing calculations 
with rational numbers and frequently, the integers can become rather large 

• these calculations are time consuming so the algorithm is not so efficient 

• we can avoid calculations with rational numbers by doing the computations in the domain of polynomials with 
integer coefficients 

Example of GCD in Q[x] 

• we want to calculate gcd(a, 6)where 

a = x 8 + x 6 - 3x 4 - 3x 3 + 8x 2 + 2x - 5 

b = 3x 6 + 5x 4 - Ax 2 - 9a; + 21 

• by using the algorithm GCDPQ, we get the following remainders r 

T\ 

r-i 

Ti 543589225 

• so the polynomials aand bare co-prime since their greatest common divisor is 1 

2.5.3 GCD of polynomials with integer coefficients 

• now we consider polynomials from Z[x]m one variable a;with integer coefficients 

• pseudo-division of polynomials a(x),b(x)for which dega(a;) > deg6(a;)is defined by 

lcof(6(a;)) m -™ +1 a(a;) = pquot(a(a;), b(x))b(x) + prem(a(x),6(x)) 

where m = dega(x),n = deg6(a;)and lcof(fo(a;))is the leading coefficient of the polynomial b(x), i.e., the 
coefficient of the n-th power of x, so that the pseudo-quotient pquot(a(a;), 6(a;))and the pseudo-remainder 
prem(a(a;),6(a;))are also from Z[x] 



5 4^ 

~9 X + 


1 2 1 

t:X — - 
9 3 


117 o 

-IT 


441 


233150 
6591 X 


102500 
2197 


1288744821 



• then < degprem(a(a;),&(a;)) < deg&(a;) 

• the polynomial a(x)is primitive if all its coefficients are mutually co-prime 

a(x) = conta(a;) ppa(a;) 

• conta(a;)is the greatest common divisor of all coefficients of the polynomial and ppa(a;)is the primitive part of 
the polynomial a, hence, a(x) = conta(a;)ppa(a;) 

• Euclidean algorithm for polynomial reminder sequence 

gcd := GCDPRS(a(x), b(x)) := 

[suppose that degree of polynomial a is greater or equal to the 
degree of polynomial b, i.e., deg(a) >= deg(b) 
algorithms used: 

prem(a, b) - pseudo-remainder of polynomial a with polynomial b 
pp(a) - primitive part of the polynomial a 
gcdi(j, k) - gcd of two integers j, k] 

1. A := pp(a) 
B := pp(b) 

2. while B != do 

r := prem(A, B) 
A := B 
B := r 
od 

3. return gcdi(cont (a) , cont(b)) pp(a) 

• coefficients and partial results grow very quickly — this algorithm is also not efficient 

• it is possible to calculate the primitive part of the remainder in each step, however, the calculation of the 
primitive part requires the calculation of many greatest common divisors of integer coefficients which can often 
be large 

• modular approach: perform calculations modulo a prime p 

— apply homomorphism to coefficients 

— calculate the greatest common divisor of the transformed polynomials modulo p 

— use the Chinese remainder algorithm for reconstructing the coefficients of the greatest common divisor 
back in the integers 

Example of GCD in Z[x] 

• we want to calculate gcd(a,6), where 

a = x s + x 6 - 3a; 4 - 3a; 3 + 8x 2 + 2x - 5 



b = 3x 6 + 5a; 4 - 4a; 2 - 9x + 21 

using the algorithm GCDPRS, we get following pseudo-remainders r 

pseudo-reminder primitive part 

n =-15x*+3x z -9 -5a; 4 + x 1 - 3 

r 2 = 15795a; 2 + 30375a; - 59535 13a; 2 + 25a; - 49 

r 3 = 1254542875143750a; - 1654608338437500 4663a; - 6150 

r 4 = 12593338795500743100931141992187500 1 



2.6 Resultant 

• let a, bbe polynomials with coefficients from the domain O 

• the domain Ocan be e.g., integers or polynomials in other variables 

n m 



i=0 



i=0 



the Sylvester matrix of polynomials a, bis the matrix 



/ a„ a n -i ■ ■ ■ a\ a 

a„ a n -i ■ ■ ■ a\ a 



\ 





• • • a n a n -\ ■ ■ ■ a\ ao 

• • • a„ a„-i ■ ■ ■ a\ ao 

b m 6 m _i • • • 6i 6 • • • 

b m 6 m _i • • • bi b • • • 




V o 



b m b m -i ■■ ■ foi bo 
b m 6 m _i • • • foi bo J 



• where mrows are constructed from the coefficients of the polynomial aand nrows from the coefficients of the 
polynomial b 

• the resultant of the polynomials aand fcwith respect to the variable x, denoted by Res(a;, a, b), is the determinant 
of the Sylvester matrix 

• using the special structure of the Sylvester matrix, we can calculate the resultant by the following recursive 
algorithm 

• Algorithm 

res := RES(x, a, b) := 
[a, b are polynomials 
algorithms used: 

rem(a, b) - remainder after dividing a by b 
deg(a) - degree of the polynomial a 

lcof(a) - leading coefficient of the polynomial a, i.e., the 
coefficient of x"deg(a)] 

1. n := deg(a) 
m := deg(b) 

2. if n > m then res :=(-l)~(n m) RES(x, b, a) 

else lc := lcof(a) 

if n = then res := lc"m 
else r := rem(b, a) 

if r = then res := 
else p := deg(r) 

res := lc"(m - p) RES(x, a, r) 
fi 
fi 
fi 



• this algorithm gradually replaces rows in the Sylvester matrix derived from the coefficients of the polynomial 
foby rows derived from the coefficients of the remainder r = rem(b, a); these new rows are computed as linear 
combinations of the rows of the original matrix 



• after the above replacement, the matrix has only zeroes below the diagonal in its first m — pcolumns while 
the diagonal has the values lc; the matrix block consisting of the last n — prows and n — pcolumns is the 
Sylvester matrix of the polynomials a,r; thus, in general case of the algorithm, we can use the formula 
Res(x,a,b) := lc m ~ p Res(x,a,r) 

2.7 Solving polynomial equations 

• solving linear systems is easy as only the inversion of a matrix is required 

• system of polynomial equations 

/i(xi,...,x„) = 

Jmt^l ; ■ ■ ■ ; %n) — U 

• Grobner basis according to lexicographic ordering of variables x\ < X2 < ■■■ < a;„(very high time complexity) 

• the system of polynomial equations and its Grobner basis have the same solution 

• example 



r 2 a 2 

/j = xz - xy - 4x - - 

h = y 2 z + 2x + ^ 

f 3 = x 2 z + y 2 + -x 



Grobner basis of these polynomials with variables ordering x < y < z 



64 , 432 „ 168 2 354 8 

91 = Z+ 65 X -~65 X + ~65 X ~-^ X+ l 

2 8 4 54 o 8 2 17 

92 = V ~13 X + 13 X ~13 X + 26 X 

C Z i A O Zi. o 

,3 = x"--x * + 2x 3--x 2 +, + - 

• this is in "triangular" form 

• a;can be determined by solving the 3rd equation; after substituting x = a;,into the 2nd equation, we can 
determine y, etc. 



one approximate solution is (-0.128475,0.321145,-2.356718) 



• solving (even numerically) of a polynomial equation in one variable is much simpler than solving a system of 
polynomial equations 

• solving a polynomial system is transformed into the successive solution of equations in one variable; see the 
example (see chapter4.6.2) 



2.8 Differentiation 

• differentiation algorithm is very simple 

• for differentiation, the following rules are sufficient: 

(a + b)' = a' + b' 
(a - b)' = a' - b' 
(a b)' = a'b + a b' 
(a/b)' = (a'b - a b')/b~2 
f(g(t))' = f (g(t)) g'(t) 

rules for differentiation of functions like In, exp, sin, cos, /(these can be implemented using a table look up) 

• (see chapter4.6.2) 

2.9 Summation 

• given the sequence a\, a-i, ■■■, a n , ..., what is the sum 



S(n) = ^ai 



i=l 



• discrete analogue of integration 

2.9.1 Simple example 

• let 



then we can express the ratio 



• where we have denoted 



™ 1 1 

S(n) = V] ,i.e.,Qj = 

^i(i + 2) i(i + 2) 



(n — l)(n + 1) q(n)p(n) 



a n -\ n(n + 2) p(n — l)r(n) 



p(n) = n + 1 
q(n) = n — 1 
r(n) = n + 2 



the idea is to try to express the sum 5(n)as 



q(n + l) 

S{n) = --— a„/(n) (2.1) 

p(n) 

n X ft \ 



(n + l)n(n + 2)' 
• where /(n)is a polynomial in n; note that 

a n = S(n) - S(n - 1) (2.2) 



• substituting (2.1) into (2.2), we obtain a recurrence relation for /(n): 

n + 1 = nf{n) - (n + 2)/(n - 1) 

• to solve the recurrence relation, we need to know the degree of the polynomial f(n) 

• we can rewrite (2.3) as 



(2.3) 



n + 1 = (n- (n + 2)) 



/(„)+/(„-!) 



+ {n + {n + 2)) 



/(„) - f( n - 1) 



(2.4) 



• introducing 



f(n) = c k n k + c k -m k 1 + . . . + c 



and substituting this formula into (2.4), we obtain 



n + 1 = (k - 2)c k n k + 0(n k - 1 ) 

• from which it follows that k < 2 (if k > 2, then the previous equation implies that c k = 0) 

• therefore, /(n)is a polynomial of at most degree two 

f(n) = C2n 2 + c\n + c 

• the solution of (2.3) is then 



CO = 


= x e R 


C\ = 


5 + Qx 



Cl 



4 
3 + a; 



where ads an arbitrary real parameter 
• the value of ads obtained from the initial condition 5(0) = 0, which gives x = Oand hence, 



/(") 



3n 2 + 5n 



• the final solution is 



S(n) = V 



3n + hn 



■^ i(i + 2) 4(n 2 + 3n + 2) 



2.9.2 Gosper algorithm 

• given the sequence a,, what is the sum 



S ( n ) = ^2 a i 



i=l 

• we know that a, = S(i) — S(i — 1) 

• let us assume that S(n)/S(n — l)is a rational function 

• the quotient of two succesive terms of the sequence can be expressed as 

a n S(n) — S(n — 1) 



a n -\ S(n — 1) — S(n — 2) 

S(n)/S(n - 1) - 1 



1 - S(n - 2)/S(n - 1) 

• thus, a„/a„_iis also a rational function 

• on the basis of the following lemma, the quotient can be always written as 

a n _ p(n)q(n) 
a„_i p(n - l)r(n) 

• where p,q,r&re polynomials for which gcd(q(n),r(n + j)) = lfor all j > 0( gcdis the (see chapter2.5) ) 

• Lemma A rational function a(n)/6(n)can be always written as 

a(n) p(n)q(n) 

b(n) p(n — l)r(n) ' 

where p, q, rare polynomials in nand 

gcd(g(n),r(n + j)) = l,Vj > 0, 

• Theorem Let S(n)/S(n — l)be a rational function as above. Then 

p(n) 



f(n) = S(n) 



q(n + l)a n 



is a polynomial for which 



Gosper algorithm 



p(n) = q(n + l)f(n) — r(n)f(n — 1). 



(S(n),b):=G0SPER( a(n)) 

[assumes that a(n)/a(n-l) is a rational function, 
if S(n)/S(n-1) is also a rational function then b=true 
and S(n) is computed, otherwise b=false 
algorithms used: 

num(a) - numerator of the rational function a 

den(a) - denominator of the ratinal function a 

Res(x, p, q) - resultant of polynomials p, q with respect to the 

variable x 
gcd(p,q) - greatest common divisor of the polynomials p and q 



deg(p(n)) - degree of the polynomial p(n) 

coef(p(n), i) - coefficient of the ith power of n in the polynomial 
p(n)] 

1. b := true 

2. if a(n) = then S(n) := 0; return fi 

3. p(n) := 1 

q(n) := num(a(n) / a(n-l)) 
r(n) := den(a(n) / a(n-l)) 

4. while (Res(n, q(n) , r(n+j)) has nonnegative integer root j = jO) do 

g(n) := gcd(q(n), r(n+j0)) 
q(n) := q(n)/g(n) 
r(n) := r(n)/g(n-j0) 

p(n) := p(n) g(n) g(n-l) ... g(n-j0+l) 
od 

5. lp := deg(q(n+l) + r(n)) 
lm := deg(q(n+l) - r(n)) 

if lp <= lm then k := deg(p(n))-lm 
else 

kO := 2 (-lp coef(q, lp) - coef(q, lp-1) + coef(r, lp-1))/ 
(coef(q, lp) + coef(r, lp)) 

if (kO is integer) then k := max(k0, deg(p(n) )-lp+l) 

else k := deg(p(n) )-lp+l 

fi 
fi 
if k < then b := false; return fi; 

6. solve recurrence relation p(n)=q(n+l)f (n)-r (n)f (n-1) 
with initial condition f (l)=p(l)/q(2) 

for f(n)=ck n~k + ... + cO 

if (solution does not exist) then b := false; return fi; 

7. S(n) := q(n+l) a(n) f(n)/p(n); 
return 

2.9.3 Examples using the Gosper algorithm 



V 



„ i .-, 1 Tim (b^ + Cj + c) 



^n;=i(^+ci+c) c -d 



mx m+2 - (m + l)x m+1 + x 
nx" = 



£< 






(see chapter4.6.5) 



2.10 Integration 

• as in the summation problem, the basic structure of the algorithm is the following: 

— a structural theorem gives the general form of the solution 

— need to determine the degree bounds on the polynomial appearing in the solution 

— finally, determine the coefficients of these polynomials 



2.10.1 Integration of rational functions 

• Theorem (Rothstein) Let A(x),B(x)he polynomials with B(a;)square free (i.e., there are no squared or higher 
degree factors in its factorization) and suppose that the degree of Ais less than the degree of B( Aeg(A) < 
deg(fl)). Then 



(A{x) ^ 



where c\, ...,c„are the distinct roots of the polynomial R(c) = Res(x,A(x) — cB'(x),B(x))and Vi = gcd(A(x) — 
CiB'(x),B(x)) for i = l, ...,n( Res is the (see chapter2.6) ). 

2.10.2 Integration of elementary transcendental functions 

• extend the rational functions by the natural logarithm and exponential function (see (see chapter2.1.2) ) 

• Theorem (Liouville) Let Kbe a differential field and /be from K. Then, an elementary extension of the field 
K, which has the same field of constants as K&nd contains an element gsuch that g' = /, exists if and only if 
there exist constants c\, ..., c„from KsmA functions u,u\, ...,u„from Ksuch that 



f = v ! + Y t cp- 

Z J q, . 



i=\ 



i.e. 



f = u + ^2ci\ni 



i=\ 



• Risch 1968-1969 - first decision procedure for the integration of elementary transcendental and algebraic func- 
tions; the procedure determines if the integral exists within a given class of functions and if so, then calculates 
the value of the integral 



2.10.3 Integration examples 



/: 



7 ^dx=\Mx 1 -2) 



I 



x 3 +2 



arctan 



+ - 



ln(x 2 - </2x + </2 2 

C = 77= 

6^2 2 
(iNwf) j ln(a;+ff2) 



^/2 2 V3 



3S/2 2 



2*e* a ln(*) + ^ + HX) - 2 



+ 



x (ln(a;) 2 + x) 2 
21n(x)/x + l/x + l' 



ln(a;) 2 + x 



d x 



e x In (a;) 



ln(x) 



ln(a;) 2 + x 



+ ln(ln(a;) 2 + x) 



(see chapter4.6.5) 



2.11 Ordinary differential equations 

• Kovacic algorithm for differential equations of the type 

ay" + by 1 + cy = 
where a, b, care polynomials in ^(implemented in Macsyma) 

• Singer algorithm for differential equations of the type 

y (n) + c n - 1 y (n - 1) + ... + ay' + c y = 
where c„_i, ...,coare rational or algebraic functions (not implemented until now due to its massive complexity) 

• algorithms that look for Liouville type solutions (i.e., rational functions extended by algebraic functions, 
exponential functions and integrals) or decide if such a solution exists 

• (see chapter4.6.5) 

2.12 Polynomial factorization 

• consider polynomials with integer coefficients (it is simple to convert a polynomial with rational coefficients 
into a problem in this domain) 

• factorization splits a given polynomial a(a;)mto the product of polynomials cij(x): 

n 

a(x) = Y[aj(x) 

• a polynomial a(x)is "square free", i.e., it does not have a factor which is the second or higher power of a 
polynomial, if and only if gcd(a(x) , a' (x)) = 1, that is, there does not exist a polynomial other than the 
constant lthat divides the polynomial a(a;)and its derivative 

• if the polynomial a(x)is not "square free" then 

a(x) = b(x) 2 c(x) 

a'(x) = 2b(x)b'(x)c(x)+b(x) 2 c'(x) 
and gcd(a(x) , a' (x)) = b(x)d(x) 

• Berlekamp algorithm for factorization of " square free" polynomials in one variable modulo a prime p 

• Berlekamp-Hensel algorithm for factorization of " square free" polynomials in one variable with integer coeffi- 
cients 

• Kronecker algorithm for factorization of polynomial in more variables 

• these algorithms exhibit exponential complexity 

• there exists algorithms with only polynomial complexity 

• uses of factorization: solving of polynomial equations, integration of rational functions, etc. 

• example 



8a; 3 y + 40x 2 y 4 + 2x 2 y 2 z + 4x 2 y 

+ 4x 2 z 2 + 10xy 5 z + 20xy 4 + 20xy 3 z 

+ xyz 3 + 2xz 2 + 5y 4 z 3 + 10y 3 z 2 

= (2xy + z 2 )(4:X + yz + 2)(x + 5y 3 ) 

(see chapter4.6.2) 



2.13 Quantifier elimination 

• quantified formula in prenex form over a real field: 

G = (Qf+iXf+i) . . . (Q r x r )F(xi,. ..,x r ) 

where Qiis a general quantifier (for all) or existential quantifier (there exists) and Fis a logical combination of 
polynomial equations and inequalities in the variables x\, ...,x r 

• /variables are free (not quantified) and r — /variables are quantified 

• there exists a quantifier free formula H'm free variables x\, ...,2;/ which is equivalent to the formula G 

• quantifier elimination is the procedure which transforms the quantified formula Ginto the quantifier free formula 
H 

• around 1930 Tarski proved that quantifier elimination is possible 

• first quantifier elimination method was proposed by Tarski 1951, however its complexity cannot be bound by 
any tower of exponentials 

• cylindrical algebraic decomposition method by Collins 1975 made a breakthrough having "only" double expo- 
nential complexity 

• improvements in the method of partial cylindrical algebraic decomposition: Collins and Hong 1991, programm 
QEPCAD 

• example: the quantified formula 



(VSi) (VS 2 ) (C> 0) A 

[(0 < Si < 1 A < S 2 < 1) => 
-2SiS 2 C 3 + 3SiS 2 C 2 + S1C 2 
-2SiC + S 2 C 2 - 2S 2 C + 1 > 0] 



• is equivalent to the quantifier free formula 

Example with QEPCAD program 

• formula is quantified for all TGX and TGY; 
inputs , outputs 



o<c<l 



Quantifier Elimination 

in 

Elementary Algebra and Geometry 

by 

Partial Cylindrical Algebraic Decomposition 

Version 10 (Interactive) 
June 1992 

by 

Hoon Hong 
(hhongOrisc . uni-linz . ac . at) 
Research Institute for Symbolic Computation 



<<QEPCAD>> Enter an informal description between '[' and ']': 
[Wendr M2, 1. factor;=0] 

<<QEPCAD>> Enter a variable list: 
(A,B,TGX,TGY) 

<<QEPCAD>> Enter the number of free variables: 
2 

<<QEPCAD>> Enter a prenex formula: 

(A TGX) (A TGY) 

( A + B ;= 1 / A + B i= -1 / A - B ;= 1 / A - B i= -1 

/ ( A~2 TGX-2 (TGY~2 + 1) + 2 A B TGX TGY ( - TGX TGY + 1) + A TGX~2 ( 

TGY-2 + 1) + B~2 TGY2 (TGX"2 + 1) - B TGY2 (TGX"2 + 1) ;= 

))• 



<<QEPCAD>> Before Normalizations 
finish 



An equivalent quantifier-free formula: 

( A <= A B >= A B + A - 1 <= A B + A + 1 >= A \\ 
B-A+1>=0/\B-A-1<=0) 



Chapter 3 

Integrated mathematical systems 



3.1 Computer algebra systems 

• specialized systems 

- TRIGMAN 1970, celestial mechanics 

- SCHOONSCHIP 1971, quantum physics 

— CAMAL 1975, celestial mechanics, general theory of relativity 

— SHEEP 1977, general theory of relativity, very fast 

• general purpose systems 

— REDUCE 1968, classical system, delivered with sources 

- MACSYMA 1970, classical system, mainly in USA 

- MAPLE 1985, fast and good system 

— DERIVE 1988, on PC, even on pocket computers 

- MATHEMATICA 1988, nice graphics and user interface (notebook) 

— AXIOM 1991, former Scratchpad II, modern "object oriented" concepts, most complete system 

• (see chapterl.3.1) 

• for other computer algebra systems and packages, including also public systems, consult the overview at CAIN 
Netherlands or the collection of symbolic software at the University of Berkeley 

3.2 Peculiarities of programming in computer algebra systems 

• computer algebra system (CAS) can be used as an algebraic calculator 

• one can program in a CAS using a high level language 

• one can get VERY LARGE algebraic expressions quite easily (see chapter??) 

• it is hard to guess the memory and CPU time requirements for a given calculation 

• often calculations need to be done only once — once a result is known, it does not need to be calculated again 

• efficiency questions (algorithmic) 

• questions of representation which can influence the efficiency as well 

• experience is very important 

29 



3.3 Expression Swell 

EXPRESSION SWELL is a common phenomenon of exact computations in which the size of numbers and 
expressions involved in a calculation grows dramatically as the calculation progresses. 

INTERMEDIATE EXPRESSION SWELL is an important special case of expression swell in which, during 
the middle stages of a calculation, intermediate expressions can expand substantially, but the final results of 
the calculation are comparatively simple. 

3.3.1 An Example of Expression Swell 

• consider the Hankel matrix 



-1 1 1-1-1 1 1 -1 -1 \ 
I -1 1 1-1-1 1-1 1 1 ) 



• it has the following characteristic polynomial: 

c(A) = A 9 + A 8 - 40A 7 - 24A 6 + 240A 5 + 144A 4 
= A 4 (A + 6) (A 4 - 5A 3 - 10A 2 + 36A + 24) 

• the Hankel matrix has four zero eigenvalues, one eigenvalue is -6, and the other four eigenvalues are roots of 
the seemingly simple looking quartic polynomial 

A 4 - 5A 3 - 10A 2 + 36A + 24 

• here is one root of this quartic polynomial 



\ 



(18 ( ^VWni + 25i36\ 2 ' s _ 465 (ie^pm + 25i36y> s + ls56 \ \ s _ 2 J_l V 3 2 JLJ +999 ( ieVWni + 2s 

\ V d 21 I V i *> 1 ) \ I / 16 V807Ti t 2R13R Y' 3 V 3 ' 

, , . 1/3 I SB ( l6vf"' + M) 2/3 +465 ( l6V807Ti + ^Ufl) 1/3 +371 T 

I 16 VWTT i I 25136 \ I/3 / V 3 2 ~ I V 3 2 ~ I 

^ 3 27 > y ( 16V807TJ , 25136 y S 



6\/2 



/ 36 ( !6v|07Ti + 25jl6^/ 3 +465 ^6VWn i + 2S13^y^ + 3712 



I 16 V8071 ; i 25136 *\ 
^ 3 T 27 j 



12 



5 

+ 4 



the other three roots are similar in structure. 



3.3.2 An Example of Intermediate Expression Swell 

• the left-hand side of the (tensor) Bianchi identity for a symmetric connection is 



Kj hk\p + Kj kp\h + Kj ph\k i 

where K is the Riemann curvature tensor 



• expanding in terms of Christoffel symbols of the second kind, one obtains 

_ r I r #19 r #27 r t r #25 r l r #19 r #22 _ r t r #20 

L #19 h L j #27 L p k~ rL jh,#25 L p k " r x #19 #22^ j i 1 ? A; x j #20,h L p * 

, F * r #18 , r * r #10 r #16 _ r * r #10 r #12 _ T l r #11 

"•" A j #18, A; 1 p ft "T" A #10 k L j #16 i p ft i #10#12 i j k L p ft i jA,#ll i p ft 

, r ( r #9 , r * r #1 r #7 r I r #1 r #3 r £ r #2 r i r #18 

"I - ^ j #9,^^ A: p + i #lft i j #7 L k p i #l#3 i j h L k p L jh,#2 L k p + L j P,#18 L k ft 

, r * r #10 r #15 _ r * r #10 r #13 _ r t r #11 _ r #20 r i 

1~ i #10#15 i j p L k ft i #10p i j #13 i i; ft i j#ll,p i A ft A ft k L j p,#20 

, r * r #9 , r l r #1 r #7 _ r * r #1 r #4 r i r #19 r #27 

"•" L #9 A.fc 1 j p" l " i #lft i #7 A 1 ^ p 1 #lt 1 #4 ft 1 ^ p "T" A #19 ft 1 #27 k L j p 

r t r #20 _ r * r #2 r i r #19 _ r * r #26 r #19 

" 1 ~ i #20A: i j p,ft L #2h,k L j p" l " i #19ft i i p,A A #19 #26 i ft A^j p 

r #26 r i r #19 _ r #26 r * r #19 r i r #19 _ r I r #18 

" l " i #19 ft 1 #26 A: 1 j p A #19 A; 1 #26 ft 1 j p "•" L #19 ft.A: 1 j p i #18ft i j p,k 

_ r * r #10 r #16 _ r i r #10 r #15 r * r #10 _ r #15 r i r #10 

1 #10 A; 1 #16 h L j p L #10k L j p,ft "•" L #10 ft 1 #15 A: 1 j p L #10 A; i #15ft i j p 

_p £ r #10, r #9 r ^ _ r ^ r #9 _ r i r #l r #7 

i #10A;,ft i i p^ft p L jk,#9 i #9p i j A, ft L #1 ft 1 #7 p 1 j * 

_ r * r #19 r #27 r i r #25 r * r #19 r #23 _ r * p #20 

L #19 ft 1 #27 p L j k ~ rL #25 h,p L j k " r x #19 p 1 #23 ft^j A; i #20p,ft i j A; 

, r * r #10 r #13 r * r #11 r i r #10 _ r #12 r * r #10 

" ri #10p i #13 ft^j k ~ rL #11 h L j k,p " r x #10 p 1 j A,h i #10 h L #12 p L j * 

r #12 r * r #10 r * r #10 _ r * r #1 r i r #6 r #1 

" l " i #10 p L #12 h L j A "T" i #10 p.ft 1 ^ A; i #lft i j k,p ~*~ L #1 #6 i ft p 1 j A; 

_ r #6 r * r #i , r #6 r ^ r # J _ r * r* 1 o.r <■ r #1 r # 4 

X #l ft i #6p i j k^ L #l p i #6ft i j A; i #lft,p i i t +1 #U L #4 p 1 j ft 

_ r i r #25 _ r ^ r #19 r #23 r i r #2 _ r ^ r #19 

^#25 k L j h,p i #19p i #23 k L j ft " r ^ #2 p L j h,k L #19 p L j h,k 

, r #22 r <? r #19 _ r #22 r <! r #19 _ r I y #19 , r I y #18 
+ i #19 A; i #22p i j ft i #19 p 1 #22 k L j ft i #19p,A; i i ft "+" L #18 p.A; 1 j ft 

, r ^ F #10 r #16 _ r l r #10 r #13 _ r ^ r #11 r ^ r #1 

"I" ^ #10 A; ± #16 p L j ft i #10p i #13 A; 1 ^ ft i #llA;,p i i ft-|- i #lA; i i ft,p 

, r #3 r I r #1 r #3 r <! r #1 , r ^ r #1 r I r #4 r #1 

+ i #l A; i #3p i i ft J^*! p i #3A; i i h + L #1 k,p L j ft J^ #1 k L ft p 1 j #4 

r #25 r i r l r #23 r #19 _ r #2 r i 

"I" 1 ft k L j #25,p "T" A #19 p L h k L j #23 A ft p L j#2,k 

• this sum contains 72 terms, each of which is a product of 2 or 3 Christoffel symbols, for a total of 180 Christoffel 
symbols 

• however, upon simplifying this expression by consistently renaming the dummy indices, the simple result of 
zero is obtained, which verifies the identity 

3.3.3 Expression Size 

• Integer: (number of base {3 digits comprising n) 



• Rational number: 



• General expression: 



H [v, n = 



A/>(^) =M (m)+M (n) 



number of terms 

number of operators and atomic operands 

number of characters 

etc. 



3.3.4 Another Example of Intermediate Expression Swell 

• stages in the expression swell analysis of the computation of the characteristic polynomial of a 5 x 5 general 
matrix containing "prime" 4-digit rational numbers (both the numerator and the denominator consist of 4 
digits) 

1. Initial matrix. 



9533 


9547 


9587 


9613 


9623 


9539 


9551 


9601 


9619 


9629 


9631 


9649 


9677 


9689 


9719 


9643 


9661 


9679 


9697 


9721 


9733 


9743 


9767 


9781 


9791 


9739 


9749 


9769 


9787 


9803 


9811 


9829 


9839 


9857 


9871 


9817 


9833 


9851 


9859 


9883 


9887 


9907 


9929 


9941 


9967 


9901 


9923 


9931 


9949 


9973 



The sizes of the entries in the upper Hessenberg matrix H that is similar to A. For example, the (5,4) 
entry is a rational number with a 309-digit numerator and a 313-digit denominator! 



/ 4 33 105 175 4 \ 

I A 32 104 175 4 \ 



MoW 



A 
A 




Vo 



33 
32 

33 
32 

53 
56 






105 
104 

106 
105 

113 
117 

155 
158 





175 
175 
175 
175 
187 
190 
229 
231 
309 
313 



± 

A 

13 
16 
69 
72 

143 
146 



The sizes of the coefficients of the characteristic polynomial cjj(A) = c^(A). For example, the coefficient 
of A 3 is a rational number with a 98-digit numerator and a 100-digit denominator. The numbers involved 
involved have decreased in size, but are still large! 

,. , ,.„ x5 21 4 98 a3 95 a2 92 A 88 
Mo(c*(A)) = A* + -A< + —A' + -A* + -A + - 



3.3.5 Expression Swell Analysis 

• bounds on expression size provide estimates of 

— memory needed 

— CPU time required 

to perform a given calculation. 

• bounds can be determined 

— empirically (statistical survey) 

— theoretically (worst case and best case analyses) 



Chapter 4 

Basic possibilities of integrated 
mathematical systems 

4.1 Axiom 

• inputs 

outputs 

4.1.1 Number domains 

In Axiom 

For comparison with other CAS choose from: Derive (see chapter 4.2.1) Macsyma (see chapter 4.3.1) Maple (see 
chapter 4.4.1) Mathematica (see chapter 4.5.1) Reduce (see chapter 4.6.1) 

Big integers 

In Axiom 

For comparison with other CAS choose from: Derive (see chapter 4.2.1) Macsyma (see chapter 4.3.1) Maple (see 
chapter 4.4.1) Mathematica (see chapter 4.5.1) Reduce (see chapter 4.6.1) 

• integers of arbitrary size 
23**12 

21914624432020321 

Type: Positivelnteger 

factorial 60 

8320987112741390144276341183223364380754172606361245952449277696409600 
000000000000 

Type: Positivelnteger 

bi:=23**4*37*59*101 

61700183203 

Type: Positivelnteger 



33 



• factorization of integers 
factor bi 

4 
23 37 59 101 



Type: Factored Integer 



bia:=23*ll**6 



40745903 



Type: Positivelnteger 



• integer greatest common divisor 
gcd(bi,bia) 

23 



Type: Positivelnteger 



Rational numbers 

In Axiom 

For comparison with other CAS choose from: Derive (see chapter 4.2.1) Macsyma (see chapter 4.3.1) Maple (see 
chapter 4.4.1) Mathematica (see chapter 4.5.1) Reduce (see chapter 4.6.1) 

• exact calculation with rational numbers 
1234567890/98765432 

617283945 

49382716 

Type: Fraction Integer 



rn: = l/2+2/15-64/47 

1027 

1410 

Type: Fraction Integer 



Complex numbers 

In Axiom 

For comparison with other CAS choose from: Derive (see chapter 4.2.1) Macsyma (see chapter 4.3.1) Maple (see 
chapter 4.4.1) Mathematica (see chapter 4.5.1) Reduce (see chapter 4.6.1) 

• exact calculation with complex numbers 
cn:= (2+3*%i)*(15-6*%i)+2/(2-4*%i) 



63 + 115'/.i 
2 + V.± 



• real and imaginary part 
real(cn) + %i*imag(cn) 

+ + 

167\|- 1 + 241 



Type: Fraction Complex Integer 



Type: Expression Integer 



Algebraic numbers 

In Axiom 

For comparison with other CAS choose from: Derive (see chapter 4.2.1) Macsyma (see chapter 4.3.1) Maple (see 
chapter ??) Mathematica (see chapter ??) Reduce (see chapter 4.6.1) 

sqrt2:= rootOf(sqrt2**2-2) 

sqrt2 

Type: AlgebraicNumber 

l/(sqrt2+l) 

1 



sqrt2 + 1 



Type: AlgebraicNumber 



(x**2+2*sqrt2*x+2)/(x+sqrt2) 
x + sqrt2 



Type: Fraction Polynomial AlgebraicNumber 



on gcd; 

(x**3+(sqrt2-2)*x**2-(2*sqrt2+3)*x-3*sqrt2)/(x**2-2); 

2 
x - 2x - 3 

- sqrt2 - 4 



x + 



2sqrt2 + 1 



Type: Fraction Polynomial AlgebraicNumber 



normalize 



2 2 

(2x - 4x - 6)sqrt2 + x - 2x - 3 

(2x - l)sqrt2 + x - 4 



Type: Expression Integer 



sqrt(x**2-2*sqrt2*x*y+2*y**2) 

+ + 

I 2 2 

\ I - 2x y sqrt2 + 2y + x 



• multiple algebraic extensions 
sqrt5:= rootOf(sqrt5**2-5) 

sqrt5 



cbrt3:= rootOf(cbrt3**3-3) 

cbrt3 



cbrt3**3 

3 



sqrt5**2; 
5 

cbrt3; 

cbrt3 



Type: Expression Integer 



Type: AlgebraicNumber 



Type: AlgebraicNumber 



Type: AlgebraicNumber 



Type: AlgebraicNumber 



Type: AlgebraicNumber 



S qrt(x**2+2*(sqrt5-cbrt3)*x+5-2*sqrt5*cbrt3+cbrt3**2) 

+ + 

I 2 2 

\|(- 2cbrt3 + 2x)sqrt5 + cbrt3 - 2x cbrt3 + x +5 

Type: Expression Integer 



Big floating point numbers 

In Axiom 

For comparison with other CAS choose from: Derive (see chapter 4.2.1) Macsyma (see chapter 4.3.1) Maple (see 
chapter 4.4.1) Mathematica (see chapter 4.5.1) Reduce (see chapter 4.6.1) 

rn:= - 1027/1410 



1027 
1410 

cn:=(167*%i + 241)/5 
241 + 167'/.i 



• computation with floating point numbers 
rn :: Float 

- 0.7283687943 2624113475 



en :: Complex Float 

48.2 + 33.4 '/,i 

%pi :: Float 

3.1415926535 897932385 



Type: Fraction Integer 



Type: Fraction Complex Integer 



cos(%pi :: Float) 
- 1.0 

sin 1.0 

0.8414709848 0789650665 



• computation with an arbitrary number of digits 
digits 50; 



Type: Float 



Type : Complex Float 



Type: Float 



Type : Float 



Type : Float 



Type: Positivelnteger 



%pi :: Float 



3.1415926535 8979323846 2643383279 5028841971 693993751 



cos %pi :: Float 

- 1.0 



should be cos(pi/6) = sqrt(3)/2 
cos(%pi/6 :: Float) 

0.8660254037 8443864676 3723170752 9361834714 0262690519 



Type : Float 



Type : Float 



0.75 



Type : Float 



Type: Float 



digits 20; 



Type: Positivelnteger 



• no underflows appears 
exp(-100000.1**2) 

0.1184406313 2021703038 E -4342953504 



• complex functions 
tan(1.0 + 1.0*%i) 

0.2717525853 1951171653 + 1.0839233273 386945435 '/.i 



log(1.0 + 1.0*%i) 

0.3465735902 7997265471 + 0.7853981633 9744830961 '/.i 



Type : Float 



Type : Complex Float 



Type : Complex Float 



4.1.2 Polynomials 

In Axiom 

For comparison with other CAS choose from: Derive (see chapter 4.2.2) Macsyma (see chapter 4.3.2) Maple (see 
chapter 4.4.2) Mathematica (see chapter 4.5.2) Reduce (see chapter 4.6.2) 



Basic operations 

In Axiom 

For comparison with other CAS choose from: Derive (see chapter 4.2.2) Macsyma (see chapter 4.3.2) Maple (see 
chapter 4.4.2) Mathematica (see chapter 4.5.2) Reduce (see chapter 4.6.2) 

• by default, parentheses are expanded 
pol:= (a+b+c)**4 

4 32 223 223 

c + (4b + 4a) c + (6b + 12a b + 6a )c + (4b + 12a b + 12a b + 4a )c 
+ 

4 3 2 2 3 4 

b +4ab +6ab +4ab + a 

Type : Polynomial Integer 

• differentiation 
dpol:= D(pol, a) 

3 22 23223 

4c + (12b + 12a) c + (12b + 24a b + 12a )c + 4b + 12a b + 12a b + 4a 

Type : Polynomial Integer 



D(D(pol, a), b, 2) 

24c + 24b + 24a 

Type : Polynomial Integer 

• integration 

integrate (dpol, a) 

3 22 223 32234 

4a c + (12a b + 6a )c + (12a b + 12a b+4a)c+4ab +6ab +4ab+a 

Type: Polynomial Fraction Integer 



• polynomial greatest common divisor 
(a**2-b**2)/(a**2-2*a*b+b**2) 

- b - a 

b - a 



Type: Fraction Polynomial Integer 



g := 34*x**19-91*x+70*x**7-25*x**16+20*x**3-86 

19 16 7 3 
34x - 25x + 70x + 20x - 91x - 86 

Type : Polynomial Integer 



fl:= g * (64*x**34-21*x**47-126*x**8-46*x**5-16*x**60-81) 

79 76 67 66 63 61 60 54 

- 544x + 400x - 1120x - 714x + 205x + 1456x + 1376x - 1470x 
+ 

53 50 48 47 41 37 35 
2176x - 2020x + 1911x + 1806x + 4480x + 1280x - 5824x 
+ 

34 27 24 21 19 16 15 

- 5504x - 4284x + 1586x + 1150x - 2754x + 2025x - 8820x 
+ 

12 11 9 8 7 6 5 3 

- 3220x - 2520x + 11466x + 9916x - 5670x + 4186x + 3956x - 1620x 
+ 

7371x + 6966 

Type : Polynomial Integer 



f2:= g * (72*x**60-25*x**25-19*x**23-22*x**39-83*x**52+54*x**10+81) 

79 76 71 68 67 63 61 
2448x - 1800x - 2822x + 2075x + 5040x + 1440x - 6552x 
+ 

60 59 58 55 53 52 46 

- 6192x - 5810x - 748x - lllOx + 7553x + 7138x - 1540x 
+ 

44 42 41 40 39 32 30 

- 850x - 1086x + 625x + 2002x + 2367x - 1750x - 1330x 
+ 

29 28 26 25 24 23 19 17 
1836x - 500x + 545x + 2150x + 1729x + 1634x + 2754x + 3780x 
+ 

16 13 11 10 7 3 

- 2025x + 1080x - 4914x - 4644x + 5670x + 1620x - 7371x - 6966 

Type : Polynomial Integer 



gcd(fl,f2) 



19 16 7 3 
34x - 25x + 70x + 20x - 91x - 86 

Type: Polynomial Integer 



Factorization 

In Axiom 

For comparison with other CAS choose from: Derive (see chapter 4.2.2) Macsyma (see chapter 4.3.2) Maple (see 
chapter 4.4.2) Mathematica (see chapter 4.5.2) Reduce (see chapter 4.6.2) 

• factorization is the transformation of a polynomial into a product of polynomials 
factor(a**2-b**2) 

- (b - a)(b + a) 

Type: Factored Polynomial Integer 



factor(a**2+b**2, [rootOf(i~2 + 1)]) 

1 1 

(b - - a)(b + - a) 
i i 

Type: Factored Polynomial AlgebraicNumber 



fa:= (x**2*z+y**4*z**2+5)* _ 
(x*y**3+z**2)* _ 

(-x**3*y+z**2+3)* _ 
(x**3*y**4+z**2) 

48 27 38 7 35 4 6 
yz +xz +(xy +xy - x y +3y + 5)z 
+ 

5 4 3 3 5 2 5 
(xy +xy -xy+3x)z 
+ 

4 11 69 4 38 7 34 3 3 4 

(x y - x y + (- x + 3x )y + 3x y + 5x y + 5x y - 5x y + 15)z 
+ 

67 85 6 54 333 

(x y - x y + (- x + 3x )y + 3x y )z 
+ 

7 12 4 11 47 65 4 34 32 

(- x y + 3x y + 5x y - 5x y + (- 5x + 15x )y + 15x y )z 
+ 

98 67 78 47 

(-xy +3xy)z-5xy + 15x y 

Type : Polynomial Integer 



factor fa 

23 2 3234422 

(z -xy+3)(z +xy)(z +xy)(yz +xz+5) 

Type : Factored Polynomial Integer 



Grobner bases 

In Axiom 

For comparison with other CAS choose from: Derive (see chapter ??) Macsyma (see chapter 4.3.2) Maple (see 
chapter 4.4.2) Mathematica (see chapter 4.5.2) Reduce (see chapter 4.6.2) 

polys := [45*p + 35*s - 165*b - 36, _ 

35*p + 40*z + 25*t - 27*s, _ 

15*w + 25*p*s + 30*z - 18*t - 165*b**2, _ 

- 9*w + 15*p*t + 20*z*s, _ 

w*p + 2*z*t - ll*b**3, _ 

99*w - ll*s*b + 3*b**2, _ 

b**2 + 33/50*b + 2673/10000] 



[35s + 45p - 165b - 36, 40z + 25t - 27s + 35p, 

2 3 

30z + 15w - 18t + 25p s - 165b , 20s z - 9w + 15p t, 2t z + p w - lib , 
2 2 33 2673 

99w - lib s + 3b , b + -- b + ] 

50 10000 

Type: List Polynomial Fraction Integer 



vars := [w, p, z, t, s, b] 
[w,p,z,t,s,b] 



Type: List Symbol 



groebner (polys) 

49 1143 19 1323 37 27 5 9 

[z + — b + , w + b + , t - — b + , s - - b , 

36 2000 120 20000 15 250 2 200 
31 153 2 33 2673 

p - — b , b + — b + ] 

18 200 50 10000 

Type: List Polynomial Fraction Integer 



• solving a system of polynomial equations by Grobner bases 
solve(polys, vars) 

[ 

- 9500b - 3969 3100b + 1377 - 24500b - 10287 1850b - 81 

[w= , p= , z= , t= , 

60000 1800 18000 750 

500b +9 2 

s= , 10000b + 6600b + 2673= 0] 

200 
] 

Type: List List Equation Fraction Polynomial Integer 



(see chapter4.1.4) 



4.1.3 Rational functions 



In Axiom 

For comparison with other CAS choose from: Derive (see chapter 4.2.3) Macsyma (see chapter 4.3.3) Maple (see 
chapter 4.4.3) Mathematica (see chapter 4.5.3) Reduce (see chapter 4.6.3) 

rf:= (3*a*b**2-5*a**2*b)/(a**4-2) 

2 2 

3a b - 5a b 

4 
a - 2 

Type: Fraction Polynomial Integer 



• integration 

integrate (rf, a) 

+ + 

|2 2 4 

(\|- 12'/.'/.Gl - 8'/.'/.G0 '/.'/.Gl - 12'/.'/.G0 + 9b - 27.7.G1 - 2°///.G0) 
* 

log 

(24a T/.Gl + 24a '/,'/,G0 - 60b - 50a) 
* 

+ + 

|2 2 4 

\|- 12'/.'/.Gl - 8'/.'/.G0 '/.'/.Gl - 12'/.'/.G0 + 9b 
+ 

2 2 

48a '/.'/.Gl + (120b + 100a) '/.'/.Gl + 48a '/.'/.GO + (120b + 100a) '/.'/.GO 
+ 

4 

- 27a b - 250b 
+ 

+ + 

|2 2 4 

(- \|- 12'/.'/.Gl - 8'/.'/.G0 '/.'/.Gl - 12'/.'/.G0 + 9b - 2'/.'/.Gl - 2'/.'/.G0) 
* 

log 

(24a '/.'/.Gl + 24a '/.'/.GO - 60b - 50a) 
* 

+ + 

|2 2 4 

\|- 12'/.'/.Gl - 8'/.'/.G0 '/.'/.Gl - 12'/.'/.G0 + 9b 
+ 

2 2 

- 48a '/.'/.Gl + (- 120b - 100a) '/.'/.Gl - 48a '/.'/.GO 
+ 

4 
(- 120b - 100a) '/.'/.GO + 27a b + 250b 
+ 

2 4 

4'/.'/.Gl log (96a '/.'/.Gl + (240b + 200a) '/.'/.Gl - 27a b + 250b) 
+ 

2 4 

4'/.'/.G0 log (96a '/.'/.GO + (240b + 200a) '/.'/.GO - 27a b + 250b) 
/ 
4 

Type: Union(Expression Integer,...) 



• partial fraction decomposition 

(10*x**2-ll*x-6)/(x**3-x**2-2*x) 

2 

lOx - llx - 6 

3 2 

x - x - 2x 



Type: Fraction Polynomial Integer 

padicFraction( _ 

partialFraction(numerator( 

factor (denominator ( 

Factored UnivariatePolynomial(x, Fraction Integer))) 

3 2 5 

x x - 2 x + 1 

Type: Part ialFract ion UnivariatePolynomial(x, Fraction Integer) 



4.1.4 Solving equations 

In Axiom 

For comparison with other CAS choose from: Derive (see chapter 4.2.4) Macsyma (see chapter 4.3.4) Maple (see 
chapter 4.4.4) Mathematica (see chapter 4.5.4) Reduce (see chapter 4.6.4) 

Linear systems 

In Axiom 

For comparison with other CAS choose from: Derive (see chapter 4.2.4) Macsyma (see chapter 4.3.4) Maple (see 
chapter 4.4.4) Mathematica (see chapter 4.5.4) Reduce (see chapter 4.6.4) 

solve([2*xl+x2+3*x3-9, xl-2*x2+x3+2, 3*xl+2*x2+2*x3-7], [xl, x2, x3]) 

[[xl= - l,x2= 2,x3= 3]] 

Type: List List Equation Fraction Polynomial Integer 



Nonlinear equations 

In Axiom 

For comparison with other CAS choose from: Derive (see chapter 4.2.4) Macsyma (see chapter 4.3.4) Maple (see 
chapter 4.4.4) Mathematica (see chapter 4.5.4) Reduce (see chapter 4.6.4) 

solve(x**8-8*x**7+34*x**6-92*x**5+175*x**4-236*x**3+226*x**2-140*x+46, x) 

8 7 6 5 4 3 2 
[x - 8x + 34x - 92x + 175x - 236x + 226x - 140x + 46= 0] 

Type: List Equation Fraction Polynomial Integer 



solve(log(acos(asin(x**(2/3)-b)-l))+2, x) 

[] 

Type: List Equation Expression Integer 



Nonlinear systems 

In Axiom 

For comparison with other CAS choose from: Derive (see chapter ??) Macsyma (see chapter ??) Maple (see 
chapter 4.4.4) Mathematica (see chapter 4.5.4) Reduce (see chapter 4.6.4) 

solve ( _ 



[ alpha * cl - beta * cl**2 - gamma*cl*c2 + epsilon*c3, _ 
-gamma*cl*c2 + (epsilon+theta)*c3 -eta *c2, _ 
gamma*cl*c2 + eta*c2 - (epsilon+theta) * c3], _ 
[c3, c2, cl]) 

>> Error detected within library code: 

system does not have a finite number of solutions 

You are being returned to the top level of the interpreter. 



• (see chapter4.1.2) 

4.1.5 Analytical operations 

In Axiom 

For comparison with other CAS choose from: Derive (see chapter 4.2.5) Macsyma (see chapter 4.3.5) Maple (see 
chapter 4.4.5) Mathematica (see chapter 4.5.5) Reduce (see chapter 4.6.5) 

Limits 

In Axiom 

For comparison with other CAS choose from: Derive (see chapter 4.2.5) Macsyma (see chapter 4.3.5) Maple (see 
chapter 4.4.5) Mathematica (see chapter 4.5.5) Reduce (see chapter 4.6.5) 

limit(sin(x)/x, x = 0) 

1 

Type: Union (OrderedComplet ion Expression Integer,...) 



limit((3*sin(%pi*x) - sin(3*%pi*x))/x**3, x = 0) 

3 

47„pi 

Type: Union (OrderedComplet ion Expression Integer,...) 

limit((2*x+5)/(3*x-2), x = 

2 

3 

Type: Union (OrderedComplet ion Fraction Polynomial Integer,...) 



Taylor series 

In Axiom 

For comparison with other CAS choose from: Derive (see chapter 4.2.5) Macsyma (see chapter 4.3.5) Maple (see 
chapter 4.4.5) Mathematica (see chapter 4.5.5) Reduce (see chapter 4.6.5) 

)set streams calculate 4 

series(%e**x, x = 0) 

12 13 14 5 

l+x+-x +-x + — x + 0(x ) 
2 6 24 

Type: UnivariatePuiseuxSeries (Expression Integer, x,0) 



)set streams calculate 2 

2 3 

1 + 2x + 2x + 0(x ) 



Type: UnivariatePuiseuxSeries (Expression Integer, x,0) 



Summation and Products 

In Axiom 

For comparison with other CAS choose from: Derive (see chapter 4.2.5) Macsyma (see chapter 4.3.5) Maple (see 
chapter 4.4.5) Mathematica (see chapter 4.5.5) Reduce (see chapter 4.6.5) 

sum(n**2*x**n, n) 

2 3 2 22n-l 

((n - 2n + l)x + (- 2n + 2n + l)x + n x)x 



3 2 

x - 3x + 3x - 1 



sum(cos((2*r-l)*%pi/(2*n+l)), r) 



— + (2°/.A - l)°/,pi 

> cos( ) 

— + 2n + 1 

U 



product(%e**(sin(n*x)), n) 



++-++ sin(7,A x) 

I I °/.e 
I I 

n 



Type: Expression Integer 



Type: Expression Integer 



Type: Expression Integer 



for all n,m such that fixp m let factorial(n+m)=if m l then factorial(n+m-l)*(n+m) else facto- 
rial(n+m+l)/(n+m+l); 

sum(n*2**n/factorial(n+2), n) 



n °/.A 

--+ °/.A 2 



-+ (7.A + 2) ! 



7.A 



Type: Expression Integer 



Integration 

In Axiom 

For comparison with other CAS choose from: Derive (see chapter 4.2.5) Macsyma (see chapter 4.3.5) Maple (see 
chapter 4.4.5) Mathematica (see chapter 4.5.5) Reduce (see chapter 4.6.5) 

integrate(x**2*(a+b*x)**p, x) 

32 3 33 22 22 2 3p log(b x + a) 

((b p + 3b p + 2b )x + (a b p + a b p)x - 2a b p x + 2a )'/,e 

3 3 3 2 3 3 

b p + 6b p + lib p + 6b 

Type: Union(Expression Integer,...) 



integrate (x**2*log(x**2+a**2), x) 

3 2 2 3 x 3 2 

3x log(x + a ) + 6a atan(- -) - 2x + 6a x 

a 

9 

Type: Union(Expression Integer,...) 



integrate(x*d**x*sin x, x) 

3 2 2 

(x log(d) - log(d) + x log(d) + l)sin(x) - x cos(x)log(d) 
+ 

2cos(x)log(d) - x cos(x) 
* 

x log(d) 
7.e 
/ 

4 2 

log(d) + 21og(d) + 1 

Type: Union(Expression Integer,...) 



integrate(x*sqrt(a+b*x)**p, x) 

+ + 

2 2 2 2 p log(\|b x + a ) 

((2b p + 4b )x + 2a b p x - 4a )'/,€ 

2 2 2 2 

b p + 6b p + 8b 

Type: Union(Expression Integer,...) 



integrate(2*x*%e**(x**2)*log(x)+%e**(x**2)/x+(log(x)-2)/(log(x)**2+x)**2+ 

((2/x)*log(x)4-(l/x)4-l)/(log(x)**2+x),x) 

2 2 

2 2 x 3 x 

(log(x) + x)log(log(x) + x) + y.e log(x) + (x '/,e - l)log(x) 



2 
log(x) + X 

Type: Union(Expression Integer,...) 



Ordinary differential equations 

In Axiom 

For comparison with other CAS choose from: Derive (see chapter ??) Macsyma (see chapter 4.3.5) Maple (see 
chapter 4.4.5) Mathematica (see chapter 4.5.5) Reduce (see chapter 4.6.5) 

y:= operator('y); 

Type: BasicOperator 

solve(D(y(x), x) + y(x) * sin x/cos x - 1/cos x = 0, y, x) 

[particular= sin(x) ,basis= [cos(x)]] 
Type: Union(Record(particular : Expression Integer , basis : List Expression Integer),...) 



• Bernoulli equation 

solve(x*(l-x**2)*D(y(x), x) + (2*x**2 -l)*y(x) - x**3*y(x)**3 = 0, y, x) 

5 2 4 2 

- 2x y(x) + 5x - 5x 



2 
5y(x) 

Type: Union(Expression Integer,...) 



solve(D(y(x), x, 2)+4*D(y(x), x)+4*y(x)-x*exp(x) = 0, y, x) 

- 2x x 3 
(3x - 2)*/,e C/.e ) - 2x - 2x 

[particular= ,basis= [°/e ,x °/,e ]] 

27 
Type: Union(Record(particular : Expression Integer , basis : List Expression Integer),...) 



Substitutions - pattern matching 

In Axiom 

For comparison with other CAS choose from: Derive (see chapter ??) Macsyma (see chapter 4.3.5) Maple (see 
chapter ??) Mathematica (see chapter ??) Reduce (see chapter 4.6.5) 

sincosRules:= rule _ 

(cos(x)*cos(y) == (cos(x+y) + cos(x-y))/2; _ 

cos(x)*sin(y) == (sin(x+y) - sin(x-y))/2; _ 

sin(x)*sin(y) == (cos(x-y) - cos(x+y))/2; _ 

cos(x)**2 == (1 + cos(2*x))/2; _ 

sin(x)**2 == (1 - cos(2*x))/2 ) 

V.X cos(y + x) + 7,X cos(- y + x) 

{°/X cos(x)cos(y) == , 

2 



'/.Y sin(y + x) - '/,Y sin(- y + x) 

°/,Y cos(x)sin(y) == , 

2 
- °/,Z cos(y + x) + °/Z cos(- y + x) 

°/,Z sin(x)sin(y) == , 

2 
2 cos(2x) +1 2 - cos(2x) + 1 

cos(x) == , sin(x) == } 

2 2 

Type: Ruleset (Integer , Integer , Expression Integer) 

sincosRules (al*cos(wt) + a3*cos(3*wt) + bl*sin(wt) + b3*sin(3*wt))**3 

3 3 2 2 2 2 

b3 sin(3wt) + (3bl b3 sin(wt) + 3a3 b3 cos(3wt) + 3al b3 cos(wt) )sin(3wt) 
+ 

2 2 
3b 1 b3 sin(wt) + (6a3 bl b3 cos(3wt) + 6al bl b3 cos(wt) )sin(wt) 
+ 

2 2 2 2 

3a3 b3 cos(3wt) + 6al a3 b3 cos(wt)cos(3wt) + 3al b3 cos(wt) 
* 

sin(3wt) 
+ 

3 3 2 2 2 

bl sin(wt) + (3a3 bl cos(3wt) + 3al bl cos(wt) )sin(wt) 
+ 

2 2 2 2 

(3a3 bl cos(3wt) + 6al a3 bl cos(wt)cos(3wt) + 3al bl cos(wt) )sin(wt) 
+ 

3 3 2 2 2 2 

a3 cos(3wt) + 3al a3 cos(wt)cos(3wt) + 3al a3 cos(wt) cos(3wt) 
+ 

3 3 
al cos(wt) 

Type: Expression Integer 



int:= operator('int); 

Type: BasicOperator 

intRules:= rule _ 

(int(x + :y, z) = = int(x, z) + int(y, z); _ 

int(k*x — freeOf?(k, z), z) == k*int(x, z); _ 

int(y — integer? y, z) == y*z; _ 

int(x**(?p — D(p, x) = 0), x) == x**(p+l)/(p+l) ) 

{int(y + x,z) == ' int(y.z) + 'int(x,z), int(k x,z) == k'int(x,z), 

P + 1 
P x 
int(y,z) == y z, int(x ,x) == } 

P + 1 
Type: Ruleset (Integer , Integer , Expression Integer) 



intRules int(a**2*b+a**b+3*a-5, a) 



b + 1 2 2 2 2 

a + (a b + a )int(b,a) + (3a - 5a)b + 3a - 5a 



b + 1 



Type: Expression Integer 



intRules int(a**(a+l), a) 

a + 2 



a + 2 



Type: Expression Integer 



4.1.6 Matrices 

In Axiom 

For comparison with other CAS choose from: Derive (see chapter 4.2.6) Macsyma (see chapter 4.3.6) Maple (see 
chapter 4.4.6) Mathematica (see chapter 4.5.6) Reduce (see chapter 4.6.6) 

xx:= matrix([[all, al2], [a21, a22]]) 

+all a!2+ 



+a21 a22+ 



yy:= matrix([[yl], [y2]]) 

+yl+ 

I I 
+y2+ 



determinant xx 

all a22 - al2 a21 



Type: Matrix Polynomial Integer 



Type: Matrix Polynomial Integer 



Type : Polynomial Integer 



zz:= inverse (xx)*yy 

+- al2 y2 + a22 yl+ 
1 

all a22 - al2 a21| 

I 

all y2 - a21 yl | 

1 

+all a22 - al2 a21+ 



Type: Matrix Fraction Polynomial Integer 



inverse (xx)**2 

[ 



2 
a22 + al2 a21 
[ , 

2 2 2 2 

all a22 - 2all al2 a21 a22 + al2 a21 

- al2 a22 - all al2 
] 

2 2 2 2 

all a22 - 2all al2 a21 a22 + al2 a21 



- a21 a22 - all a21 
[ , 

2 2 2 2 

all a22 - 2all al2 a21 a22 + al2 a21 

2 

al2 a21 + all 
] 

2 2 2 2 

all a22 - 2all al2 a21 a22 + al2 a21 
] 

Type: Matrix Fraction Polynomial Integer 



v := matrix([[2, -1, 1], [0, 1, 1], [-1, 1, 1]]) 

+2 - 1 1+ 

I I 

10 1 II 

I I 

+- 1 1 1+ 

Type: Matrix Integer 

eigenvectors v 

+0+ +1+ 

II II 

[[eigval= 2,eigmult= l,eigvec= [ | 1 | ] ] , [eigval= l,eigmult= 2,eigvec= [|1|]]] 

II II 

+1+ +0+ 

Type: List Union (Record (algrel: Fraction Polynomial Integer ,algmult : Integer ,algvec: List Matrix Fraction Po 

4.1.7 Graphics 

In Axiom 

For comparison with other CAS choose from: Derive (see chapter 4.2.7) Macsyma (see chapter 4.3.8) Maple (see 
chapter 4.4.7) Mathematica (see chapter 4.5.7) Reduce (see chapter 4.6.8) 

2D Graphics 

In Axiom 

For comparison with other CAS choose from: Derive (see chapter 4.2.7) Macsyma (see chapter 4.3.8) Maple (see 
chapter 4.4.7) Mathematica (see chapter 4.5.7) Reduce (see chapter 4.6.8) 



• function of one parameter 
draw(sin(%e~x), x = 0..%pi) 

• graph of several functions (Bessel functions J(n,x),n = 0,2,5) 
draw([besselj(0, x), besselj(2, x), besselj(5, x)], x = 0..10, _ 
title == "Bessel Functions BesselJ(n, x)") 

3D Graphics 

In Axiom 

For comparison with other CAS choose from: Derive (see chapter 4.2.7) Macsyma (see chapter 4.3.8) Maple (see 
chapter 4.4.7) Mathematica (see chapter 4.5.7) Reduce (see chapter 4.6.8) 

• graph of a function with 2 parameters 
draw(sin(%pi*sin(x+y)), x = -3..3, y = -3. .3) 

• graph of another function with 2 parameters 

draw(tan(x*y), x = -2/3*%pi..2/3*%pi, y = -2/3*%pi..2/3*%pi) 

Parametric plots 

In Axiom 

For comparison with other CAS choose from: Derive (see chapter 4.2.7) Macsyma (see chapter 4.3.8) Maple (see 
chapter 4.4.7) Mathematica (see chapter 4.5.7) Reduce (see chapter 4.6.8) 

• surface defined parametrically 
draw(surface(sin(v), sin(2*v)*sin(u), sin(2*v)*cos(u)), _ 
u = 0..2*%pi, v = -%pi/2..%pi/2) 

4.2 Derive 

• inputs 

outputs 

4.2.1 Number domains 

In Derive 

For comparison with other CAS choose from: Axiom (see chapter 4.1.1) Macsyma (see chapter 4.3.1) Maple (see 
chapter 4.4.1) Mathematica (see chapter 4.5.1) Reduce (see chapter 4.6.1) 

Big integers 

In Derive 

For comparison with other CAS choose from: Axiom (see chapter 4.1.1) Macsyma (see chapter 4.3.1) Maple (see 
chapter 4.4.1) Mathematica (see chapter 4.5.1) Reduce (see chapter 4.6.1) 

• integers of arbitrary size 
23~12 

12 
23 



21914624432020321 



60! 

60! 



8320987112741390144276341183223364380754172606361245952449277696409600 
000000000000 



23"4 37 59 101 



4 
23 37 59 101 



61700183203 



• factorization of integers 
factor 

4 
23 37 59 101 



Rational numbers 

In Derive 

For comparison with other CAS choose from: Axiom (see chapter 4.1.1) Macsyma (see chapter 4.3.1) Maple (see 
chapter 4.4.1) Mathematica (see chapter 4.5.1) Reduce (see chapter 4.6.1) 

• precise calculation with rational numbers 
1234567890/98765432 

1234567890 
98765432 



617283945 
49382716 



1/24-2/15-64/47 



1 2 64 

2 15 47 



1027 

1410 



Complex numbers 

In Derive 

For comparison with other CAS choose from: Axiom (see chapter 4.1.1) Macsyma (see chapter 4.3.1) Maple (see 
chapter 4.4.1) Mathematica (see chapter 4.5.1) Reduce (see chapter 4.6.1) 

• precise calculation with complex numbers 
(2 + 3 #i) (15 - 6 #i) + 2/(2 - 4 #i) 

2 
(2 + 3 i) (15 - 6 i) + 

2 - 4 i 



241 167 i 

+ 

5 5 



Radicals 

In Derive 

For comparison with other CAS choose from: Axiom (see chapter 4.1.1) Macsyma (see chapter 4.3.1) Maple (see 
chapter ??) Mathematica (see chapter ??) Reduce (see chapter 4.6.1) 

• computation with radicals 

(x-24-2 sqrt(2) x+2)/(x+sqrt(2)) 

2 
x +2 2 x + 2 



x + 2 



x + 2 



sqrt(x~2-2 sqrt(2) x y+2 y~2) 

2 2 

(x -2 2xy + 2y) 



2 2 2 x - 4 y 

4 



Big floating point numbers 

In Derive 

For comparison with other CAS choose from: Axiom (see chapter 4.1.1) Macsyma (see chapter 4.3.1) Maple (see 
chapter 4.4.1) Mathematica (see chapter 4.5.1) Reduce (see chapter 4.6.1) 

• computation with floating point numbers 
- 1027/1410 



1027 
1410 

-0.7283687943 

(167 #i + 241)/5 
241 167 i 

5 5 

48.2 + 33.4 i 

pi 

3.141592920 



cos(pi) 

COS ( ) 



sin(l) 

SIN (1) 



0.8414709137 



• computation with arbitrary number of digits 
options, precision, digits 50 
pi 



3 . 1415926535897932384626433832795028841971693993751 



should be cos(pi/6) = sqrt(3)/2 
cos(pi/6) 



cos — 

6 



. 86602540378443864676372317075293618347140262690519 



#2-2 



2 
. 86602540378443864676372317075293618347140262690519 



0.75 



• complex functions 

options, precision, digits 10 
tan(1.0 + 1.0 #i) 

TAN (1 + 1 i) 



0.2717525893 + 1.083923333 i 

log(1.0 + 1.0 #i) 

LOG (1 + 1 i) 

0.3465735900 + 0.7853981634 i 



4.2.2 Polynomials 

In Derive 

For comparison with other CAS choose from: Axiom (see chapter 4.1.2) Macsyma (see chapter 4.3.2) Maple (see 
chapter 4.4.2) Mathematica (see chapter 4.5.2) Reduce (see chapter 4.6.2) 

Basic operations 

In Derive 

For comparison with other CAS choose from: Axiom (see chapter 4.1.2) Macsyma (see chapter 4.3.2) Maple (see 
chapter 4.4.2) Mathematica (see chapter 4.5.2) Reduce (see chapter 4.6.2) 

• by default parantheses are not expanded 
p := (a+b+c)~4 

4 
P := (a + b + c) 



expand 



43 3 222 22 3 2 

a + 4 a b+4a c+6a b +12 a bc+6a c +4ab + 12 a b c 



2 343 22 34 

+ 12 a b c +4ac +b +4b c + 6b c +4bc +c 



• derivation 
d := dif(p,a) 

d 
D := — P 

da 



3 2 2 2 3 2 2 

4a +12 a (b + c) + 12 a (b +2bc + c)+4b +12b c+12bc 



3 
+ 4 c 



ddd43 3 22 2 22 3 

— — — a +4a b + 4a c + 6a b +12 a bc + 6a c +4ab 
db db da 



2 2 343 22 34 

+ 12 a b c + 12 a b c +4ac +b +4b c + 6b c +4bc +c 



24 a + 24 b + 24 c 

• integration 
i := int(d,a) 

I := D da 



43 22 2 32 2 

a +4a (b + c)+6a (b +2bc + c)+a(4b +12b c+12bc 



3 
+ 4 c ) 



• verification 
i - p 

I - P 



4 3 2 2 3 4 

(b +4b c + 6b c +4bc +c) 



• greatest common divisor of polynomials 
(a~2-b~2)/(a~2-2 a b+b~2) 

2 2 

a - b 



2 2 

a - 2 a b + b 



a + b 
a - b 



g := 34 x~19-91 x+70 x~7-25 x'16+20 x~3-86 

f := g (64 x'34-21 x~47-126 x~8-46 x~5-16 x'60-81) 

h := g (72 x'60-25 x'25-19 x'23-22 x~39-83 x~52+54 x~10+81) 

f/h insufficient memory 

Factorization 

In Derive 

For comparison with other CAS choose from: Axiom (see chapter 4.1.2) Macsyma (see chapter 4.3.2) Maple (see 
chapter 4.4.2) Mathematica (see chapter 4.5.2) Reduce (see chapter 4.6.2) 

• factorization is transformation of a polynomial into a product of polynomials 
a~2-b~2 

2 2 

a - b 



factor 

(a - b) (a + b) 

a~2+b~2 

2 2 

a + b 



factor, complex 

(a - i b) (a + i b) 

(x~2 z+y~4 z~2+5) (-x~3 y+z~2+3) 

3 2 2 4 2 

(- x y + z + 3) (x z + y z +5) 



expand 



5 352 3 23 2 44 42 2 

-x yz-x y z -5x y + x z +3x z + y z +3y z +5z + 

15 



factor 



2 4 2 3 2 

(x z+y z + 5) (x y - (z +3)) 



for bigger expressions was insufficient memory 

4.2.3 Rational functions 

In Derive 

For comparison with other CAS choose from: Axiom (see chapter 4.1.3) Macsyma (see chapter 4.3.3) Maple (see 
chapter 4.4.3) Mathematica (see chapter 4.5.3) Reduce (see chapter 4.6.3) 

• integration 

w := int((3 a b~2-5 a~2 b)/(a~4-2),a) 

2 2 

3 a b - 5 a b 

W := da 

4 
a - 2 



2 
3/4 2 2 a - 2 3/4 

3/4 a 2 3 2 b LN 5 b 2 LN 

5 b 2 ATAN 2 

2 2 a + 2 
+ + 

4 8 8 



1/4 



a + 2 

1/4 



• verification by derivation 
dif(w,a); 

d 

— W 
da 



2 a b (5 a - 3 b) 



2 1/4 3/4 

(2a +2) (a +2 ) ( 2 a - 2 ) 



by using build the denominator is chosen 

2 1/4 3/4 

(2a +2) (a +2 ) ( 2 a - 2 ) 



4 
2 a - 4 



final result 

a b (5 a - 3 b) 



4 
a - 2 



4.2.4 Solving equations 

In Derive 

For comparison with other CAS choose from: Axiom (see chapter 4.1.4) Macsyma (see chapter 4.3.4) Maple (see 
chapter 4.4.4) Mathematica (see chapter 4.5.4) Reduce (see chapter 4.6.4) 

Linear systems 

In Derive 

For comparison with other CAS choose from: Axiom (see chapter 4.1.4) Macsyma (see chapter 4.3.4) Maple (see 
chapter 4.4.4) Mathematica (see chapter 4.5.4) Reduce (see chapter 4.6.4) 

[2x + y + 3z-9, x-2y + z + 2, 3x + 2y + 2z-7] 

[2x+y+3z-9, x-2y+z+2, 3x+2y+2z-7] 
[x=-l, y=2, z=3] 



Nonlinear equations 

In Derive 

For comparison with other CAS choose from: Axiom (see chapter 4.1.4) Macsyma (see chapter 4.3.4) Maple (see 
chapter 4.4.4) Mathematica (see chapter 4.5.4) Reduce (see chapter 4.6.4) 

• solving polynomial equations 
x~3 + 5 x~2 - x + 2 

3 2 

x +5x -x + 2 



3777 349 1/3 349 3777 1/3 5 
18 54 54 18 3 



3777 349 1/3 349 3777 1/3 5 77897 

x = + + + i 

144 432 432 144 3 3456 



349 3777 1/6 349 3777 77897 1/6 
1152 1152 3456 



3777 349 1/3 349 3777 1/3 5 

x = + + + i 

144 432 432 144 3 



349 3777 77897 1/6 77897 349 3777 1/6 
+ - 

1152 3456 3456 1152 



• multiple use of inversion functions 
LOG(ACOS(ASIN(x~(2/3) - b) - 1)) 

2/3 
LOG (ACOS (ASIN (x - b) - 1)) 



1/3 
x = - (SIN (COS (1) + 1) + b) 



1/3 
x = (SIN (COS (1) + 1) + b) 



4.2.5 Analytical operations 

In Derive 

For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Macsyma (see chapter 4.3.5) Maple (see 
chapter 4.4.5) Mathematica (see chapter 4.5.5) Reduce (see chapter 4.6.5) 



Limits 

In Derive 

For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Macsyma (see chapter 4.3.5) Maple (see 
chapter 4.4.5) Mathematica (see chapter 4.5.5) Reduce (see chapter 4.6.5) 

lim(sin(x)/x,x,0) 

SIN (x) 

lim 

x->0 x 



lim((3 sin(pi x) - sin(3 pi x))/x~3,x,0) 

3 SIN ( x) - SIN (3 x) 

lim 

x->0 3 

x 



3 
4 



lim((2x+5)/(3x-2),x,inf) 

2 x + 5 

lim 

x-> 3 x - 2 



Taylor series 

In Derive 

For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Macsyma (see chapter 4.3.5) Maple (see 
chapter 4.4.5) Mathematica (see chapter 4.5.5) Reduce (see chapter 4.6.5) 

taylor (#e~x, x, 0, 4) 

x 
TAYLOR (e , x, 0, 4) 



4 3 2 

XXX 

+ + + x + I 

24 6 2 



taylor(taylor(#e~(x+y), x, 0, 2), y, 0, 2) 



x + y 
TAYLOR (TAYLOR (e , x, 0, 2) , y , 0, 2) 



2 2 

(y + 2 y + 2) (x + 2 x + 2) 



4 



Summation and Products 



In Derive 

For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Macsyma (see chapter 4.3.5) Maple (see 
chapter 4.4.5) Mathematica (see chapter 4.5.5) Reduce (see chapter 4.6.5) 

sum(m~2 x~m,m,l,n); 



n 2 m 
m x 
m=l 



2 2 2 2 n 

x (n x - x (2 n + 2 n - 1) + (n + 1) ) x x (x + 1) 



3 3 

(x - 1) (x - 1) 



sum(cos((2 m-1) pi/(2 n+l)),m,l,r) 

r (2 m - 1) 
COS 

m=l 2 n + 1 



2 r 
SIN 

2 n + 1 



2 SIN 

2 n + 1 



product (#e~(sin(m x)),m,l,n) 

n SIN (m x) 
e 
m=l 



COT (x) / 2 + 1 / (2 SIN (x)) - COS (x (n + 1/2)) / (2 SIN (x/2)) 
e 



Integration 

In Derive 

For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Macsyma (see chapter 4.3.5) Maple (see 
chapter 4.4.5) Mathematica (see chapter 4.5.5) Reduce (see chapter 4.6.5) 

int(x d~x sin(x),x) 

x 
x d SIN (x) dx 



x x 2 x 
x SIN (x) LN (d) d x COS (x) d SIN (x) LN (d) d 
+ 

2 2 2 2 

LN (d) +1 LN (d) + 1 (LN (d) + 1) 



x x 
2 COS (x) LN (d) d SIN (x) d 
+ 

2 2 2 2 

(LN (d) + 1) (LN (d) + 1) 



4.2.6 Matrices 

In Derive 

For comparison with other CAS choose from: Axiom (see chapter 4.1.6) Macsyma (see chapter 4.3.6) Maple (see 
chapter 4.4.6) Mathematica (see chapter 4.5.6) Reduce (see chapter 4.6.6) 

x := [[a,b],[c,d]] 

a b 
X : = 

c d 



y == [[u],[v]] 



u 
Y : = 

v 



det(x) 

DET (X) 

a d - b c 



= x-(-i) • y 

-l 

Z := X Y 



d u - b v 

a d - b c 

a v - c u 

a d - b c 



l/x-2 

1 



2 

X 



2 

b c + d b (a + d) 



22 22 22 22 

a d -2abcd + b c a d -2abcd + b c 

2 
c (a + d) a + b c 



22 2222 22 

a d -2abcd + b c a d -2abcd + b c 



eigenvalues([[2,-l,l], [0,1,1], [-1,1,1]]) 

2-11 

EIGENVALUES 1 1 , e 

-111 

[e = 1, e = 2] 



4.2.7 Graphics 

In Derive 

For comparison with other CAS choose from: Axiom (see chapter 4.1.7) Macsyma (see chapter 4.3.8) Maple (see 
chapter 4.4.7) Mathematica (see chapter 4.5.7) Reduce (see chapter 4.6.8) 

2D Graphics 

In Derive 

For comparison with other CAS choose from: Axiom (see chapter 4.1.7) Macsyma (see chapter 4.3.8) Maple (see 
chapter 4.4.7) Mathematica (see chapter 4.5.7) Reduce (see chapter 4.6.8) 

• function of one parameter 

sin(exp(x)) plot, setting scale and center, plot 



• graph of several functions (Bessel functions J{n,x),n = 0,2,5) 
transfer, load, derive, bessel 

BESSEL_J(0,z) 
BESSEL_J(2,z) 
BESSEL _J(5,z) plot, setting scale a center, plot 

3D Graphics 

In Derive 

For comparison with other CAS choose from: Axiom (see chapter 4.1.7) Macsyma (see chapter 4.3.8) Maple (see 
chapter 4.4.7) Mathematica (see chapter 4.5.7) Reduce (see chapter 4.6.8) 

• graph of a function with 2 parameters 
SIN(pi SIN(x)+y) 

• graph of another function with 2 parameters 
TAN(x y) 

Parametric plots 

In Derive 

For comparison with other CAS choose from: Axiom (see chapter 4.1.7) Macsyma (see chapter 4.3.8) Maple (see 
chapter 4.4.7) Mathematica (see chapter 4.5.7) Reduce (see chapter 4.6.8) 

• parametricaly given surface 
transfer, load, derive, graphics 

ISOMETRICS([SIN(v), SIN(2 v) SIN(u), SIN(2 v) COS(u)], u,0,2 pi,20,v,-pi/2,pi/2,10) 
ISOMETRICS([SIN(v), SIN(2 v) SIN(u), SIN(2 v) COS(u)], v,-pi/2,pi/2,10,u,0,2 pi,20) 

4.3 Macsyma 

• inputs 

outputs 

4.3.1 Number domains 

In Macsyma 

For comparison with other CAS choose from: Axiom (see chapter 4.1.1) Derive (see chapter 4.2.1) Maple (see 
chapter 4.4.1) Mathematica (see chapter 4.5.1) Reduce (see chapter 4.6.1) 

Big integers 

In Macsyma 

For comparison with other CAS choose from: Axiom (see chapter 4.1.1) Derive (see chapter 4.2.1) Maple (see 
chapter 4.4.1) Mathematica (see chapter 4.5.1) Reduce (see chapter 4.6.1) 

• integers of arbitrary size 
23-12; 

21914624432020321 



60!: 



8320987112741390144276341183223364380754172606361245952449277696409600 
000000000000 



bi: 23~4*37*59*101: 



61700183203 



• factorization of integers 
factor (bi); 



23 37 59 101 



bia: 23*11~6; 



40745903 



• integer greatest common divisor 
gcd(bi,bia); 



23 



Rational numbers 

In Macsyma 

For comparison with other CAS choose from: Axiom (see chapter 4.1.1) Derive (see chapter 4.2.1) Maple (see 
chapter 4.4.1) Mathematica (see chapter 4.5.1) Reduce (see chapter 4.6.1) 

• exact calculation with rational numbers 
1234567890/98765432; 

617283945 
49382716 



rn: 1/2+2/15-64/47; 



1027 
1410 



Complex numbers 

In Macsyma 

For comparison with other CAS choose from: Axiom (see chapter 4.1.1) Derive (see chapter 4.2.1) Maple (see 
chapter 4.4.1) Mathematica (see chapter 4.5.1) Reduce (see chapter 4.6.1) 

• exact calculation with complex numbers 
en: (2+3*%i)*(15-6*%i)+2/(2-4*%i); 

2 

(15 - 6 °/.i) (3 °/.i + 2) + 

2-4 °/.i 



en: rectform(cn); 



167 '/.i 241 

+ 

5 5 



Algebraic numbers 

In Macsyma 

For comparison with other CAS choose from: Axiom (see chapter 4.1.1) Derive (see chapter 4.2.1) Maple (see 
chapter ??) Mathematica (see chapter ??) Reduce (see chapter 4.6.1) 

• an algebraic number 
algebraic :true ; 



true 

tellrat(sqrt2**2-2); 

2 

[sqrt2 - 2] 

rat(l/(sqrt2+l)); 

/R/ sqrt2 - 1 

ogcd:gcd; 

spmod 
gcd: 'algebraic; 

algebraic 



rat((x**2+2*sqrt2*x+2)/(x+sqrt2)); 

/R/ x + sqrt2 

rat((x**3+(sqrt2-2)*x**2-(2*sqrt2+3)*x-3*sqrt2)/(x**2-2)); 



2 
x - 2 x - 3 
/R/ 
x - sqrt2 



gcd:ogcd; 

spmod 

radcan(sqrt(x**2-2*sqrt2*x*y+2*y**2)); 

2 2 

sqrt(2 y - 2 sqrt2 x y + x ) 

untellrat(sqrt2); 

[] 

• multiple algebraic extensions 

tellrat(sqrt5**2-5,cbrt3**3-3); 



3 2 

[cbrt3 - 3, sqrt5 - 5] 



rat(cbrt3**3); 

/R/ 3 
rat(sqrt5**2); 

/R/ 5 
rat(cbrt3); 

/R/ cbrt3 



rat(sqrt5); 

/R/ sqrt5 

radcan(sqrt(x**2+2*(sqrt5-cbrt3)*x+5-2*sqrt5*cbrt3+cbrt3**2)); 



2 2 

sqrt(x + (2 sqrt5 - 2 cbrt3) x - 2 cbrt3 sqrt5 + cbrt3 + 5) 



untellrat(sqrt5,cbrt3); 



[] 



Big floating point numbers 

In Macsyma 

For comparison with other CAS choose from: Axiom (see chapter 4.1.1) Derive (see chapter 4.2.1) Maple (see 
chapter 4.4.1) Mathematica (see chapter 4.5.1) Reduce (see chapter 4.6.1) 

rn: - 1027/1410; 

1027 
1410 

en: (167*i + 241)/5; 

167 °/.i + 241 



• computation with floating point numbers 
bfloat(rn); 

- 7. 28368794326241 13475b- 1 

bfloat(cn); 

3.34bl '/,i + 4.82bl 

bfloat(%pi); 

3. 1415926535897932385b0 

cos(bfloat(%pi)); 



- l.ObO 

sin(l.ObO); 

8 .4147098480789650665b- 1 



• computation with an arbitrary number of digits 
fpprec: 50$ 
bfloat(%pi); 

3 . 1415926535897932384626433832795028841971693993751b0 



cos(bfloat(%pi)); 

- l.ObO 



should be cos(pi/6) = sqrt(3)/2 
cos(bfloat(%pi/6)); 

8. 6602540378443864676372317075293618347140262690519b-! 



%"2; 

7.5b-l 

fpprec: 20$ 

• with normal defaults, underflows are converted to 
exp(-100000.1**2); 

0.0 



• complex functions 

bfloat(rectform(tan(1.0b0 + 1.0b0*%i))); 

1.0839233273386945435b0 '/.i + 2. 7175258531951171653b-! 



bfloat(log(1.0b0 + 1.0b0*%i)); 

7. 8539816339744830962b-! */,i + 3. 4657359027997265471b-! 



4.3.2 Polynomials 

In Macsyma 

For comparison with other CAS choose from: Axiom (see chapter 4.1.2) Derive (see chapter 4.2.2) Maple (see 
chapter 4.4.2) Mathematica (see chapter 4.5.2) Reduce (see chapter 4.6.2) 

Basic operations 

In Macsyma 

For comparison with other CAS choose from: Axiom (see chapter 4.1.2) Derive (see chapter 4.2.2) Maple (see 
chapter 4.4.2) Mathematica (see chapter 4.5.2) Reduce (see chapter 4.6.2) 

• by default, parentheses are not expanded 
pol: (a+b+c)~4; 

4 
(c + b + a) 



pol: expand(pol); 

4 3 322 2223 2 

c +4bc +4ac +6b c + 12 a b c +6a c +4b c+12ab c 

2 34 32234 

+ 12 a bc + 4a c + b +4ab +6a b +4a b + a 



• differentiation 

dpol: diff(pol, a); 

3 2 2 2 2 3 2 

4c +12 be +12ac +12b c + 24 a b c + 12 a c+4b + 12 a b 

2 3 

+ 12 a b + 4 a 



diff(pol, a, 1, b, 2); 

24 c + 24 b + 24 a 

• integration 

integrate (dpol, a); 

3 222 2 2 3 3 

4ac + 12 a b c +6a c + 12 a b c + 12 a bc + 4a c + 4ab 

2 2 3 4 
+ 6a b +4a b + a 



• verification 
%-pol; 



4 3 2 2 3 4 

c -4bc -6b c -4b c-b 



• polynomial greatest common divisor 
(a~2-b~2)/(a~2-2*a*b+b~2); 



2 2 

a - b 



2 2 

b - 2 a b + a 



ratsimp(%); 



b + a 
b - a 



g : 34*x~19-91*x+70*x~7-25*x~16-|-20*x~3-86; 



19 16 7 3 

34 x - 25 x + 70 x + 20 x - 91 x - 86 



fl : g * (64*x~34-21*x~47-126*x~8-46*x~5-16*x~60-81); 

19 16 7 3 

(34 x - 25 x + 70 x + 20 x - 91 x - 86) 



60 47 34 8 5 

(- 16 x - 21 x + 64 x - 126 x - 46 x - 81) 



f2 : g * (72*x~60-25*x~25-19*x~23-22*x~39-83*x~52+54*x~10-|-81); 

19 16 7 3 

(34 x - 25 x + 70 x + 20 x - 91 x - 86) 

60 52 39 25 23 10 
(72 x - 83 x - 22 x - 25 x - 19 x + 54 x +81) 



gcd(fl,f2); 



19 16 7 3 
34 x - 25 x + 70 x + 20 x - 91 x - 86 



Factorization 

In Macsyma 

For comparison with other CAS choose from: Axiom (see chapter 4.1.2) Derive (see chapter 4.2.2) Maple (see 
chapter 4.4.2) Mathematica (see chapter 4.5.2) Reduce (see chapter 4.6.2) 

• factorization is the transformation of a polynomial into a product of polynomials 
factor(a~2-b~2); 

- (b - a) (b + a) 



gfactor(a~2+b~2); 

(b - Xi a) (b + Xi a) 

fa: expand( (x~2*z+y~4*z~2+5)* (x*y~3+z~2)* (-x~3*y+z~2+3)* (x~3*y~4+z~2) ); 

48 2 7 386 7 6 356 46 6 545 

y z +x z +x y z +xy z -x y z +3y z +5z +x y z 

335 5 5 25 4 11 4 694 484 384 

+ xyz-xyz+3xz+xy z-xyz-xyz+3xyz 

74 344 34 3 4 4 673 853 

+ 3xy z +5x y z +5xy z -5x yz +15z +x y z -x y z 

643 543 333 7 12 2 4 11 2 472 

-xyz+3xyz+3xyz-xy z+3xy z+5xyz 

652 442 342 32 98 67 

-5x y z -5x y z + 15 x y z + 15 x y z -x y z + 3x y z 

7 8 4 7 

-5x y + 15 x y 



fact or (fa); 



23 2 3234422 

(z -x y + 3)(z +xy)(z +x y)(y z +x z + 5) 



Decomposition 

In Macsyma 

For comparison with other CAS choose from: Axiom (see chapter ??) Derive (see chapter ??) Maple (see chapter 
4.4.2) Mathematica (see chapter 4.5.2) Reduce (see chapter 4.6.2) 

• decomposing a polynomial into simpler polynomials which when composed, produce the original polynomial 
polydecomp(x~6+9*x~5+52*x~4+177*x~3+435*x~2+630*x+593, x); 



2 
3 2 x - 9 

[x + 25 x + 210 x + 593, , 2 x + 3] 

4 



polydecomp(x~44-2*x~3*y + 3*x~2*y~2 + 2*x*y~3 + y~4 + 2*x~2*y 4-2*x*y~2 + 2*y~3 + 5* x ~2 
+ 5*x*y + 6*y~2 + 5*y + 9, x); 

2 2 2 
x +11 3y +4y+x +10 
[ s s y + 2 x] 

4 2 



Grobner bases 

In Macsyma 

For comparison with other CAS choose from: Axiom (see chapter 4.1.2) Derive (see chapter ??) Maple (see 
chapter 4.4.2) Mathematica (see chapter 4.5.2) Reduce (see chapter 4.6.2) 

• system of polynomials 

polys : [45*p + 35*s - 165*b - 36, 35*p + 40*z 4- 25*t - 27*s, 15*w + 25*p*s 4- 30*z - 18*t 
- 165*b~2, - 9*w + 15*p*t + 20*z*s, w*p + 2*z*t - ll*b~3, 99*w - ll*s*b + 3*b~2, b~2 4- 
33/50*b + 2673/10000]; 

[35 s + 45 p - 165 b - 36, 40 z + 25 t - 27 s + 35 p, 

2 

30 z + 15 w - 18 t + 25 p s - 165 b,20sz-9w+15pt, 

3 2 2 33 b 2673 

2tz + pw-llb,99w-llbs + 3b,b + + ] 

50 10000 



vars : [w, p, z, t, s, b]; 

[w, p, z, t, s, b] 
total degree ordering grobner _tot_order:true; 

true 
gpolys:grobner (polys, vars); 

/R/ [- 500 b + 200 s - 9, 1850 b - 750 t - 81, - 9500 b - 60000 w - 3969, 

2 
10000 b + 6600 b + 2673 
- 24500 b - 18000 z - 10287, - 3100 b + 1800 p - 1377, ] 

10000 



solve the Grobner basis solve(gpolys, vars); 



190 sqrt(ll) '/,i + 139 62 sqrt(ll) '/,i + 59 

[[w = , p = 

10000 300 

490 sqrt(ll) '/,i + 367 148 sqrt(ll) °/,i - 461 

z = , t = 

3000 500 



15 sqrt(ll) '/,i - 39 12 sqrt(ll) '/.i - 33 

s = , b = ], 

50 100 

190 sqrt(ll) '/,i - 139 62 sqrt(ll) '/,i - 59 

[w = , p = , 

10000 300 

490 sqrt(ll) '/,i - 367 148 sqrt(ll) '/,i + 461 

z = , t = 

3000 500 

15 sqrt(ll) '/,i + 39 12 sqrt(ll) '/,i + 33 
s = s b = ]] 

50 100 



4.3.3 Rational functions 

In Macsyma 

For comparison with other CAS choose from: Axiom (see chapter 4.1.3) Derive (see chapter 4.2.3) Maple (see 
chapter 4.4.3) Mathematica (see chapter 4.5.3) Reduce (see chapter 4.6.3) 

rf: (3*a*b~2-5*a~2*b)/(a~4-2); 

2 2 

3 a b - 5 a b 



4 
a - 2 



• integration 

integrate (rf, a); 



1/4 

a - 2 

5 log( ) 

2 2 1/4 

3 log(a + sqrt(2)) b 3 log(a - sqrt(2)) b a + 2 
b ( + 

4 sqrt(2) 4 sqrt(2) 1/4 

4 2 

a 

5 atan( ) 

1/4 
2 



) 

1/4 
2 2 



• partial fraction decomposition 

partfrac((10*x~2-ll*x-6)/(x~3-x~2-2*x), x); 



5 3 2 

x + 1 x x - 2 



4.3.4 Solving equations 

In Macsyma 

For comparison with other CAS choose from: Axiom (see chapter 4.1.4) Derive (see chapter 4.2.4) Maple (see 
chapter 4.4.4) Mathematica (see chapter 4.5.4) Reduce (see chapter 4.6.4) 

Linear systems 

In Macsyma 

For comparison with other CAS choose from: Axiom (see chapter 4.1.4) Derive (see chapter 4.2.4) Maple (see 
chapter 4.4.4) Mathematica (see chapter 4.5.4) Reduce (see chapter 4.6.4) 

solve([2*xl+x2+3*x3-9, xl-2*x2+x3+2, 3*xl+2*x2+2*x3-7], [xl, x2, x3]); 

[[xl = - 1, x2 = 2, x3 = 3]] 



Nonlinear equations 

In Macsyma 

For comparison with other CAS choose from: Axiom (see chapter 4.1.4) Derive (see chapter 4.2.4) Maple (see 
chapter 4.4.4) Mathematica (see chapter 4.5.4) Reduce (see chapter 4.6.4) 

• using decomposition to solve high degree polynomials 

solve(x~8-8*x~7+34*x~6-92*x~5+175*x~4-236*x~3+226*x~2-140*x+46, x); 

sqrt(- sqrt(4 sqrt(3) + 3) '/.i - 3) - sqrt(2) 
[x = 

sqrt(2) 

sqrt(- sqrt(4 sqrt(3) + 3) '/.i - 3) + sqrt(2) 

x = , 

sqrt(2) 

sqrt(sqrt(4 sqrt(3) + 3) '/.i - 3) - sqrt(2) 

x = - , 

sqrt(2) 

sqrt(sqrt(4 sqrt(3) + 3) '/.i - 3) + sqrt(2) 

x = , 

sqrt(2) 

sqrt(sqrt(4 sqrt(3) - 3) + 3) '/.i - sqrt(2) 
x = , 



sqrt(2) 

sqrt(sqrt(4 sqrt(3) - 3) + 3) '/.i + sqrt(2) 

x = , 

sqrt(2) 

sqrt(3 - sqrt(4 sqrt(3) - 3)) '/.i - sqrt(2) 

x = - , 

sqrt(2) 

sqrt(3 - sqrt(4 sqrt(3) - 3)) '/.i + sqrt(2) 
x = ] 

sqrt(2) 

• multiple use of inversion functions 
assume (b l 0); 

[b > 0] 

solve(log(acos(asin(x~(2/3)-b)-l))+2, x); 

- 2 3/2 1/3 - 2 

[x = (b + sin(cos(°/,e ) + 1)) , x = - sqrt(b + sin(cos(7,e ) + 1))] 



forget (b l 

4.3.5 Analytical operations 

In Macsyma 

For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Derive (see chapter 4.2.5) Maple (see 
chapter 4.4.5) Mathematica (see chapter 4.5.5) Reduce (see chapter 4.6.5) 

Limits 

In Macsyma 

For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Derive (see chapter 4.2.5) Maple (see 
chapter 4.4.5) Mathematica (see chapter 4.5.5) Reduce (see chapter 4.6.5) 

limit(sin(x)/x, x, 0); 



limit((3*sin(%pi*x) - sin(3*%pi*x))/x~3, x, 0); 

3 
4 '/.pi 

limit((2*x+5)/(3*x-2), x, %inf); 

2 °/„inf + 5 

3 y.inf - 2 



Taylor series 

In Macsyma 

For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Derive (see chapter 4.2.5) Maple (see 
chapter 4.4.5) Mathematica (see chapter 4.5.5) Reduce (see chapter 4.6.5) 

taylor(%e~x, x, 0, 4); 

2 3 4 

XXX 

111 1 + x + — + — + — + . . . 

2 6 24 



taylor(%e~(x+y), x, 0, 2, y, 0, 2); 

2 2 2 

y y y y i 2 

111 1 + y + — + . . . + ( — + y + 1 + . . .) x + ( — +-+-+. . .) x 
2 2 4 2 2 



%"2: 



2 2 

111 l+2y+2y + . . . + (4 y +4y+2+...)x 



2 2 

+ (4y +4y + 2+...)x +... 



Summation and Products 

In Macsyma 

For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Derive (see chapter 4.2.5) Maple (see 
chapter 4.4.5) Mathematica (see chapter 4.5.5) Reduce (see chapter 4.6.5) 

load(nusuml)$ 

closedform(sum(m~2*x~m, m, 1, n)); 

n + 1 2 2 2 2 

x (n x +(-2n -2n+l)x + n +2n+l) x(x+l) 



3 3 

(x - 1) (x - 1) 



closedform(sum(cos((2*m-l)*%pi/(2*n+l)), m, 1, r)); 



\ °/.pi (2 m - 1) 
> cos( ) 

/ 2 n + 1 

m = 1 



closedform(product(%e~(sin(m*x)), m, 1, n)); 

7,i n x + °/,i x 
°/.i °/.i °/.e 

exptC/.e, 

°/,i n x + °/,i x °/,i n x °/,i x 

2 '/,e - 2 7,e 2 '/,e - 2 

Xi x 

°/.i °/.e yi 
+ + ) 

y.i x y,i x 

2 y,e - 2 2 y.e - 2 



closedform(sum(m*2~m/(m+2)!, m, 1, n)); 

n 
2 2 

1 

(n + 2)! 



Integration 

In Macsyma 

For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Derive (see chapter 4.2.5) Maple (see 
chapter 4.4.5) Mathematica (see chapter 4.5.5) Reduce (see chapter 4.6.5) 

integrate(x~2*(a+b*x)~p, x); 

32 322 22 3p log(b x + a) 

(b (p + 3 p + 2) x + a b (p + p) x - 2 a b p x + 2 a ) °/,e 

3 3 2 

b (p + 6 p + 11 p + 6) 



integrate(x"2*log(x"2+a"2), x); 

3 2 

3 x x - 3 a x 

3 2 2 2 (a atan(-) + ) 

x log(x + a ) a 3 



integrate(x*d~x*sin(x), x); 

3 2 log(d) x 

(((log (d) + log(d)) x - log (d) + 1) y.e sin(x) 

2 log(d) x 4 2 

+ ((- log (d) - 1) x + 2 log(d)) y.e cos(x))/(log (d) + 2 log (d) + 1) 



integrate (x*sqrt(a+b*x)~p, x); 

p log(b x + a) 

2 2 2 2 

(b (2 p + 4) x + 2 a b p x - 4 a ) '/.e 

2 2 
b (p + 6 p + 8) 



radcan(%); 



p/2 2 2 2 2 

(bx + a) ((2 b p + 4b)x +2abpx-4a) 

2 2 2 2 

b p +6b p + 8b 



integrate(2*x*%e~(x~2)*log(x)+%e~(x~2)/x+(log(x)-2)/(log(x)-2+x)-2+((2/x)*log(x) + (l/x)4-l)/(log(x)-2+x), 

x); 

2 
2 log(x) x 

log(log (x) + x) + °/,e log(x) 

2 
log (x) + X 



Ordinary differential equations 

In Macsyma 

For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Derive (see chapter ??) Maple (see 
chapter 4.4.5) Mathematica (see chapter 4.5.5) Reduce (see chapter 4.6.5) 

ode('diff(y, x) + y * sin(x)/cos(x) - l/cos(x), y, x); 

y = cos(x) (tan(x) + °/c) 



• Bernoulli equation 

ode(x*(l-x~2)*'diff(y, x) + (2*x~2 -l)*y - x~3*y~3, y , x ); 

log(x + 1) log(x - 1) 

+ 

2 2 

x y. e 

y = 

3 5 3 
x +3x 3x +5x + 15 x 
sqrtC/.c - 2 ( )) 

3 15 



radcan(%); 



sqrt(5) sqrt(x - 1) x sqrt(x + 1) 

y = 

5 
sqrt(2 x + 5 °/,c) 



ode('diff(y, x, 2)+4*'diff(y, x)+4*y-x*exp(x), y, x); 



x 
(3 x - 2) y.e - 2 x 

y = + c/ (k2 x + Xkl) °/.e 

27 



Substitutions - pattern matching 

In Macsyma 

For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Derive (see chapter ??) Maple (see 
chapter ??) Mathematica (see chapter ??) Reduce (see chapter 4.6.5) 

matchdeclare([x, y], true)$ 

letsimp((cos(x)*cos(y), (cos(x+y) + cos(x-y))/2)); 

cos(y + x) + cos(y - x) 
2 

letsimp((cos(x)*sin(y), (sin(x+y) - sin(x-y))/2)); 

sin(y + x) + sin(y - x) 
2 

letsimp((sin(x)*sin(y), (cos(x-y) - cos(x+y))/2)); 

cos(y - x) - cos(y + x) 
2 

letsimp((cos(x)"2, (1 + cos(2*x))/2)); 

cos (2 x) + 1 

2 

letsimp((sin(x)~2, (1 - cos(2*x))/2)); 

1 - cos (2 x) 
2 



letsimp (expand ((al*cos(wt) + a3*cos(3*wt) + bl*sin(wt) + b3*sin(3*wt))~3)); 



3 3 2 2 

b3 sin (3 wt) + 3 a3 b3 cos (3 wt) sin (3 wt) 

2 2 2 2 

+ 3 bl b3 sin(wt) sin (3 wt) + 3 al b3 cos(wt) sin (3 wt) 

2 2 
+ 3 a3 b3 cos (3 wt) sin(3 wt) + 6 a3 bl b3 sin(wt) cos(3 wt) sin(3 wt) 

2 2 
+ 6 al a3 b3 cos(wt) cos(3 wt) sin(3 wt) + 3 bl b3 sin (wt) sin(3 wt) 

2 2 
+ 6 al bl b3 cos(wt) sin(wt) sin(3 wt) + 3 al b3 cos (wt) sin(3 wt) 

3 3 2 2 2 2 

+ a3 cos (3 wt) + 3 a3 bl sin(wt) cos (3 wt) + 3 al a3 cos(wt) cos (3 wt) 

2 2 
+ 3 a3 bl sin (wt) cos (3 wt) + 6 al a3 bl cos(wt) sin(wt) cos (3 wt) 

2 2 3 3 2 2 

+ 3 al a3 cos (wt) cos (3 wt) + bl sin (wt) + 3 al bl cos(wt) sin (wt) 

2 2 3 3 

+ 3 al bl cos (wt) sin(wt) + al cos (wt) 



declare(int, linear)$ 
matchdeclare(p, is(diff(p, x) = 0))$ 
tellsimp(int(x~p, x), x~(p+l)/(p+l)); 

[intrulel, simpargsl] 



tellsimp(int(l, x), x); 

[intrule2, intrulel, simpargsl] 



int(a~2*b+a~b+3*a-5, a); 



b + 1 3 2 
a a b 3 a 
+ + 5 a 

b + 1 3 2 



int(a~(a+l), a); 

a + 2 



a 



a + 2 



4.3.6 Matrices 

In Macsyma 

For comparison with other CAS choose from: Axiom (see chapter 4.1.6) Derive (see chapter 4.2.6) Maple (see 
chapter 4.4.6) Mathematica (see chapter 4.5.6) Reduce (see chapter 4.6.6) 

xx: matrix([all, al2], [a21, a22]); 

[ all al2 ] 
[ ] 

[ a21 a22 ] 



yy: matrix([yl], [y2]); 



determinant (xx) ; 



zz: xx~~(-l).yy; 



[ yi 1 
[ ] 
[ y2 ] 



all a22 - al2 a21 



a22 yl al2 y2 

all a22 - al2 a21 all a22 - al2 a21 

all y2 a21 yl 

all a22 - al2 a21 all a22 - al2 a21 



•(-2); 



Col 1 



a22 



al2 a21 



2 2 

(all a22 - al2 a21) (all a22 - al2 a21) 



a21 a22 



all a21 



2 2 

(all a22 - al2 a21) (all a22 - al2 a21) 



al2 a22 



all al2 



2 2 

(all a22 - al2 a21) (all a22 - al2 a21) 



Col 2 



al2 a21 



all 



+ 



[ (all a22 - al2 a21) (all a22 - al2 a21) ] 



factor (%); 



a22 + al2 a21 



al2 (a22 + all) 



2 2 

(all a22 - al2 a21) (all a22 - al2 a21) 



a21 (a22 + all) 



al2 a21 + all 



2 2 

(all a22 - al2 a21) (all a22 - al2 a21) 



v : matrix([2, -1, 1], [0, 1, 1], [-1, 1, 1]); 

[2 -11] 

[ ] 

[0 11] 

[ ] 

[ - 1 1 1 ] 

eigenvectors (v); 

[[2, 1, [[0, 1, 1]]], [1, 2, [[1, 1, 0]]]] 



4.3.7 Code generation 

In Macsyma 

For comparison with other CAS choose from: Axiom (see chapter ??) Derive (see chapter ??) Maple (see chapter 
??) Mathematica (see chapter ??) Reduce (see chapter 4.6.7) 

• package gentran 

gentranout ("gentst_mac.f"); 



gentst_mac.f 



m : matrix( [ 18*COS(Q3)*COS(Q2)*M30*P~2 - 9*SIN(Q3)~2*P~2*M30 - SIN(Q3)~2*J30Y + 
SIN(Q3)~2*J30Z + P~2*M10 + 18*P~2*M30 + J10Y + J30Y, 9*COS(Q3)*COS(Q2)*M30*P~2 
- SIN(Q3)~2*J30Y + SIN(Q3)~2*J30Z - 9*SIN(Q3)~2*M30*P~2 + J30Y + 9*M30*P~2, -9*SIN(Q3)*SIN(Q2' 

]. 

[ 9*COS(Q3)*COS(Q2)*M30*P~2 - SIN(Q3)~2*J30Y + SIN(Q3)~2*J30Z - 9*SIN(Q3)~2*M30*P~2 
+ J30Y + 9*M30*P~2, -SIN(Q3)~2*J30Y + SIN(Q3)"2*J30Z - 9*SIN(Q3)~2*M30*P~2 + J30Y 
+ 9*M30*P~2, ], 

[ -9*SIN(Q3)*SIN(Q2)*M30*P~2, 0, 9*M30*P~2 + J30X ] ); 



2 2 2 2 

matrix([- 9 m30 p sin (q3) + j30z sin (q3) - j30y sin (q3) 



2 2 2 

+ 18 m30 p cos(q2) cos(q3) + 18 m30 p + mlO p + j30y + jlOy, 

2 2 2 2 2 

- 9 m30 p sin (q3) + j30z sin (q3) - j30y sin (q3) + 9 m30 p cos(q2) cos(q3) 

2 2 

+ 9 m30 p + j30y, - 9 m30 p sin(q2) sin(q3)] , 

2 2 2 2 2 
[- 9 m30 p sin (q3) + j30z sin (q3) - j30y sin (q3) + 9 m30 p cos(q2) cos(q3) 

2 2 2 2 2 

+ 9 m30 p + j30y, - 9 m30 p sin (q3) + j30z sin (q3) - j30y sin (q3) 

2 2 2 

+ 9 m30 p + j30y, 0], [-9 m30 p sin(q2) sin(q3) , 0, 9 m30 p + j30x]) 

• we know that matrix m is symmetric. We wish to generate FORTRAN code to compute numerical values for 
matrix m, and its inverse matrix, minv. 

mm : copymatrix(m)$ 



gentran( literal( "C", cr, "C — Calculate Matrix Values — ", cr, "C", cr ) ); 



lent st mac.f 



for i:l thru 3 do for j:i thru 3 do gentran( m[eval(i),eval(j)] : eval(m[ij]) 



gentran( literal( "C", cr, "C — Assign Non-Zero Matrix Values to Temporary ", "Variables 

— ", cr, "C", cr ) ); 

gentst_mac.f 

for i:l thru 3 do for j:i thru 3 do if m[ij]#0 then ( var : tempvar (false), markvar(var), m[ij] : 
var, mp,i] : var, gentran( eval(var) : m[eval(i),eval(j)] ) )$ 



• matrix m contains m; 

[ tO tl t2 ] 
[ ] 

[ tl t3 ] 
[ ] 

[ t2 t4 ] 



minv : m 



'(-I) 



gentran( literal( "C", cr, "C — Calculate Inverse Matrix Values — ", cr, "C", cr ) ); 

gentst_mac.f 
for i:l thru 3 do for j:i thru 3 do gentran( minv[eval(i),eval(j)] : eval(minv[i j]) )$ 



gentran( literal( "C", cr, "C — Copy Entries Across Main Diagonals — ", cr, "C", cr ), for i:l 
thru 3 do for j:i+l thru 3 do ( m[j,i] : m[ij], minv[j,i] : minv[ij] ) ); 



gentst_mac.f 
gentranshut ("gentst_mac.f ' ); 



true 



Generated FORTRAN program 

Generated program 



C 

C Calculate Matrix Values 

C 

M(l,l)=-(9*M30*P**2*SIN(Q3)**2)+J30Z*SIN(Q3)**2-(J30Y*SIN(q3)**2)+ 
. 18*M30*P**2*COS(Q2)*COS(q3)+18*M30*P**2+M10*P**2+J30Y+J10Y 

M(l,2)=-(9*M30*P**2*SIN(Q3)**2)+J30Z*SIN(Q3)**2-(J30Y*SIN(q3)**2)+ 
. 9*M30*P**2*COS(Q2)*COS(q3)+9*M30*P**2+J30Y 

M(1,3)=-(9*M30*P**2*SIN(Q2)*SIN(Q3)) 

M(2,2)=-(9*M30*P**2*SIN(Q3)**2)+J30Z*SIN(q3)**2-(J30Y*SIN(q3)**2)+ 
. 9*M30*P**2+J30Y 

M(2,3)=0 

M(3,3)=9*M30*P**2+J30X 
C 

C Assign Non-Zero Matrix Values to Temporary Variables 

C 

T0=M(1,1) 

T1=M(1,2) 

T2=M(1,3) 

T3=M(2,2) 

T4=M(3,3) 
C 

C Calculate Inverse Matrix Values 

C 

MINV (1,1)= (T3*T4) / ( (T0*T3-T1**2) *T4- (T2**2*T3) ) 



MINV ( 1 , 2) =- ( (T1*T4) / ( (T0*T3-T1**2) *T4- (T2**2*T3) ) ) 
MINV ( 1 , 3) =- ( (T2*T3) / ( (T0*T3-T1**2) *T4- (T2**2*T3) ) ) 
MINV (2 , 2) = (T0*T4-T2**2) / ( (T0*T3-T1**2) *T4- (T2**2*T3) ) 
MINV(2 , 3) = (T1*T2) / ( (T0*T3-T1**2) *T4- (T2**2*T3) ) 
MINV(3,3)=(T0*T3-T1**2)/((T0*T3-T1**2)*T4-(T2**2*T3)) 
C 

C Copy Entries Across Main Diagonals 

C 

DO 25001 1=1,3 

DO 25002 J=I+1,3 
M(J,I)=M(I,J) 
MINV(J,I)=MINV(I,J) 
25002 CONTINUE 

25001 CONTINUE 

4.3.8 Graphics 

In Macsyma 

For comparison with other CAS choose from: Axiom (see chapter 4.1.7) Derive (see chapter 4.2.7) Maple (see 
chapter 4.4.7) Mathematica (see chapter 4.5.7) Reduce (see chapter 4.6.8) 

2D Graphics 

In Macsyma 

For comparison with other CAS choose from: Axiom (see chapter 4.1.7) Derive (see chapter 4.2.7) Maple (see 
chapter 4.4.7) Mathematica (see chapter 4.5.7) Reduce (see chapter 4.6.8) 

• function of one parameter 

plot(sin(%e-x), x, 0, %pi, "Time", "Signal", false); 

• graph of several functions (Bessel functions J(n, x), n = 0, 2, 5) 

plot([bessel_j[0](x), bessel_j[2](x), bessel_j[5] (x)], x, 0, 10, false, "Value", "Bessel Functions 
BesselJ(n, x)"); 

3D Graphics 

In Macsyma 

For comparison with other CAS choose from: Axiom (see chapter 4.1.7) Derive (see chapter 4.2.7) Maple (see 
chapter 4.4.7) Mathematica (see chapter 4.5.7) Reduce (see chapter 4.6.8) 

• graph of a function with 2 parameters 
plot3d(sin(%pi*sin(x+y)), x, -3, 3, y, -3, 3); 

• graph of another function with 2 parameters 
plot3d(tan(x*y), x, -2/3*%pi, 2/3*%pi, y, -2/3*%pi, 2/3*%pi); 

Parametric plots 

In Macsyma 

For comparison with other CAS choose from: Axiom (see chapter 4.1.7) Derive (see chapter 4.2.7) Maple (see 
chapter 4.4.7) Mathematica (see chapter 4.5.7) Reduce (see chapter 4.6.8) 

• surface defined parametrically 

plotsurf([[sin(v), sin(2*v)*sin(u), sin(2*v)*cos(u)]], u, 0, 2*%pi, v, -%pi/2, %pi/2); 



Contour maps 

In Macsyma 

For comparison with other CAS choose from: Axiom (see chapter ??) Derive (see chapter ??) Maple (see chapter 

4.4.7) Mathematica (see chapter 4.5.7) Reduce (see chapter 4.6.8) 

• contour map of a function with 2 parameters 
contourplot(sin(x + cos(y)), x, 0, 3/2*%pi, y, 0, 3/2*%pi); 

• contour map of another function with 2 parameters 
contourplot(tan(x*y), x, -%pi, %pi, y, -%pi, %pi); 

• contour map of yet another function with 2 parameters 
contourplot(%e~(-sqrt(x~2 + y~2))*cos(atan(x/y)), x, -1, 1, y, -1, 1); 

4.3.9 Graphical presentation of formulas 

In Macsyma 

For comparison with other CAS choose from: Axiom (see chapter ??) Derive (see chapter ??) Maple (see chapter 

4.4.8) Mathematica (see chapter 4.5.8) Reduce (see chapter 4.6.9) 

4.4 Maple 

• inputs 

outputs 

4.4.1 Number domains 

In Maple 

For comparison with other CAS choose from: Axiom (see chapter 4.1.1) Derive (see chapter 4.2.1) Macsyma (see 
chapter 4.3.1) Mathematica (see chapter 4.5.1) Reduce (see chapter 4.6.1) 

Big integers 

In Maple 

For comparison with other CAS choose from: Axiom (see chapter 4.1.1) Derive (see chapter 4.2.1) Macsyma (see 
chapter 4.3.1) Mathematica (see chapter 4.5.1) Reduce (see chapter 4.6.1) 

• integers of arbitrary size 
23-12; 

21914624432020321 



60!; 



832098711274139014427634118322336438075417260636124595244927769640960\ 
0000000000000 



bi := 23-4*37*59*101; 

bi := 61700183203 



• factorization of integers 
ifactor( bi ); 



(23) (37) (59) (101) 



bia := 23*11~6; 



bia := 40745903 



• greatest common divisor of integers 
gcd( bi, bia ); 



23 



Rational numbers 

In Maple 

For comparison with other CAS choose from: Axiom (see chapter 4.1.1) Derive (see chapter 4.2.1) Macsyma (see 
chapter 4.3.1) Mathematica (see chapter 4.5.1) Reduce (see chapter 4.6.1) 

• precise calculation with rational numbers 
1234567890 / 98765432; 

617283945 

49382716 



rn := 1/24-2/15-64/47; 



1027 



1410 



Complex numbers 

In Maple 

For comparison with other CAS choose from: Axiom (see chapter 4.1.1) Derive (see chapter 4.2.1) Macsyma (see 
chapter 4.3.1) Mathematica (see chapter 4.5.1) Reduce (see chapter 4.6.1) 

• precise calculation with complex numbers 
en := (2+3*I)*(15-6*I)+2/(2-4*I); 

en := 241/5 + 167/5 I 



Big floating point numbers 

In Maple 

For comparison with other CAS choose from: Axiom (see chapter 4.1.1) Derive (see chapter 4.2.1) Macsyma (see 
chapter 4.3.1) Mathematica (see chapter 4.5.1) Reduce (see chapter 4.6.1) 

rn := - 1027/1410; 

1027 

rn := - 

1410 

en := (167*1 + 241)/5; 

en := 241/5 + 167/5 I 



• computation with floating point numbers 
evalf( rn ); 

-.7283687943 



evalf( en ); 



48.20000000 + 33.40000000 I 



Pi; 



Pi 



evalf( Pi ); 



3.141592654 



evalf( cos( Pi ) ); 



i(i- ); 



-l. 



• computation with arbitrary number of digits 
evalf( rn, 20 ); 



. 728368794326241 13475 



Digits := 50; 



Digits := 50 



evalf( Pi ); 



3 . 1415926535897932384626433832795028841971693993751 



should be cos(pi/6) = sqrt(3)/2 
evalf( cos( Pi/ 6), 50); 



. 86602540378443864676372317075293618347140262690520 



"~2; 



0.7500000000 



Digits := 10; 



Digits := 10 



• underflow gives an error 
exp(-100000.1**2); 

Error, (in evalf /exp/general) argument too large 



• complex functions 
tan(1.0 + I ); 

0.2717525853 + 1.083923327 I 



log( 1.0 + I ); 



0.3465735903 + .7853981634 I 



4.4.2 Polynomials 

In Maple 

For comparison with other CAS choose from: Axiom (see chapter 4.1.2) Derive (see chapter 4.2.2) Macsyma (see 
chapter 4.3.2) Mathematica (see chapter 4.5.2) Reduce (see chapter 4.6.2) 



Basic operations 

In Maple 

For comparison with other CAS choose from: Axiom (see chapter 4.1.2) Derive (see chapter 4.2.2) Macsyma (see 
chapter 4.3.2) Mathematica (see chapter 4.5.2) Reduce (see chapter 4.6.2) 

• by default parantheses are not expanded 
pol := (a+b+c)~4; 

4 
pol := (a + b + c) 

expand( pol ); 

2 2 24443 3 

12 a bc + 12ab c+12abc +a +b +c +4a b + 4a c 

2222 3 33 22 3 

+ 6a b +6a c +4ab +4ac +4b c + 6b c +4bc 



• derivation 

dpol := diff( expand( pol ), a ); 

2 2 3 2 2 

dpol :=24abc+12b c+12bc +4a +12 a b + 12 a c 

2 2 3 3 

+ 12ab +12ac +4b +4c 



diff( pol, a, b, b ); 

24 a + 24 b + 24 c 

• integration 
int( dpol, a ); 

2 2 243 322 

12 a bc+12ab c+12abc +a +4a b + 4a c + 6a b 

2 2 3 3 

+ 6a c +4ab +4ac 



• verification 

simplify ( " - pol ); 

4 4 3 2 2 3 

-b -c -4b c-6b c -4bc 



• greatest common divisor of polynomials 
(a~2-b~2)/(a~2-2*a*b + b~2); 



2 2 

a - b 



2 2 

a - 2 a b + b 



simplify ( " ); 

a + b 
a - b 

g := 34*x~19 - 25*x~16 + 70*x~7 + 20*x~3 - 91*x - 86; 

19 16 7 3 
g := 34 x - 25 x + 70 x + 20 x - 91 x - 86 



fl := expand( g*(64*x~34 - 21*x~47 - 126*x~8 - 46*x~5 - 16*x~60 - 81) ); 

7 16 19 3 

fl := 7371 x - 5670 x + 2025 x - 2754 x - 1620 x + 6966 

34 47 8 5 60 53 

- 5504 x + 1806 x + 9916 x + 3956 x + 1376 x + 2176 x 

66 27 24 79 50 63 

- 714 x - 4284 x + 1586 x - 544 x - 2020 x + 205 x 

21 76 41 54 15 12 
+ 1150 x + 400 x + 4480 x - 1470 x - 8820 x - 3220 x 

67 37 11 35 48 9 

- 1120 x + 1280 x - 2520 x - 5824 x + 1911 x + 11466 x 

6 61 
+ 4186 x + 1456 x 



f2 := expand( g*(72*x~60 - 25*x~25 - 19*x~23 - 22*x~39 - 83*x~52 + 54*x~10 + 81) ); 

7 16 19 3 60 

±1 := - 7371 x + 5670 x - 2025 x + 2754 x + 1620 x - 6192 x 

53 24 79 63 76 41 

+ 7553 x + 1729 x + 2448 x + 1440 x - 1800 x + 625 x 

67 11 61 25 23 39 

+ 5040 x - 4914 x - 6552 x + 2150 x + 1634 x + 2367 x 



52 10 44 42 58 71 

+ 7138 x - 4644 x - 850 x - 1086 x - 748 x - 2822 x 

29 55 68 26 32 30 

+ 1836 x - 1110 x + 2075 x + 545 x - 1750 x - 1330 x 

46 59 17 28 13 40 

- 1540 x - 5810 x + 3780 x - 500 x + 1080 x + 2002 x 

- 6966 



gcd( fl, £2 ); 



19 16 7 3 
34 x - 25 x + 70 x + 20 x - 91 x - 86 



Factorization 

In Maple 

For comparison with other CAS choose from: Axiom (see chapter 4.1.2) Derive (see chapter 4.2.2) Macsyma (see 
chapter 4.3.2) Mathematica (see chapter 4.5.2) Reduce (see chapter 4.6.2) 

• factorization is transformation of a polynomial into a product of polynomials 
factor( a~2 - b~2 ); 

(a - b) (a + b) 



factor( a~2 + b~2 ); 

2 2 
a + b 



factor( a~2 + b~2, I ); 



(a - I b) (a + I b) 



fa := expand( (x~2* z + y~4* z"2 + 5) * (x* y~3 + z~2) 
• (-x~3* y + z~2 + 3) * ( x~3* y~4 + z~2) ); 

333 48 76 11 44 844 12 27 

fa:=3xzy+yz+yzx + y zx-yzx-y zx 

534 25 554 55 835 

+ 3x z y +3x z +x z y -x z y - x z y 

6 7 353 637 634 9 8 

+ 3x zy +x z y +x z y -x z y -x zy 



32 7 4 563 946 11 24 

+ 15 x y z +3y z x - y z x -y z x +3y z x 

47 78 234 434 43 

+ 15 x y -5x y + 15 z x y +5z x y -5z x y 

265 34 472 442 843 

-5z x y +5xy z +5x y z -5x y z +3y z x 

863 6 427 46 

+ y z x +5z + 15 z +x z +3y z 



fact or ( fa ); 



2 42 342 32 2 3 

(x z + y z +5)(x y +z)(xy +z)(-z -3 + x y) 



Decomposition 

In Maple 

For comparison with other CAS choose from: Axiom (see chapter ??) Derive (see chapter ??) Macsyma (see 
chapter 4.3.2) Mathematica (see chapter 4.5.2) Reduce (see chapter 4.6.2) 

• finding of polynomial substitution which transforms a polynomial into another polynomial in further polyno- 
mials 

readlib( compoly ); 

compoly( x~6 + 9*x~5 + 52*x~4 + 177*x~3 + 435*x~2 + 630*x + 593, x ); 

3 11 2 

x + + 5/3 x, x = 25/3 + 3 x + x 

27 



compoly( x~4 + 2*x~3* y + 3*x~2* y~2 + 2*x* y~3 + y~4 + 2*x~2* y + 2*x* y~2 + 2*y~3 + 
5*x~2 + 5*x*y + 6*y"2 + 5*y + 9, x ); 

2 2 2 

x + 11/4, x = y + y + 5/2 + y x + x 



Grobner bases 

In Maple 

For comparison with other CAS choose from: Axiom (see chapter 4.1.2) Derive (see chapter ??) Macsyma (see 
chapter 4.3.2) Mathematica (see chapter 4.5.2) Reduce (see chapter 4.6.2) 

• moduls groebner and gbasis dealing with Grobner bases 
with( grobner, gbasis ); 

[gbasis] 



polys := [ 45*p + 35*s - 165*b - 36, 35*p + 40*z + 25*t - 27*s, 15*w + 25*p*s + 30*z - 18*t 
- 165*b~2, -9*w + 15*p*t + 20*z*s, w*p + 2*z*t - ll*b~3, 99*w - ll*s*b + 3*b~2, b~2 + 
33/50*b + 2673/10000 ] ; 

polys := [45 p + 35 s - 165 b - 36, 35 p + 40 z + 25 t - 27 s , 

2 

15 w + 25 p s + 30 z - 18 t - 165 b , -9w+15pt+20zs, 

3 2 2 33 2673 

wp+2zt-llb,99w-llsb+3b,b + b + ] 

50 10000 



vars := [ w, p, z, t, s, b ]; 

vars := [w, p, z, t, s, b] 

gbasis( polys, vars, tdeg ); 

[1800 p - 3100 b - 1377, 18000 z + 24500 b + 10287, 

2 
60000 w + 9500 b + 3969, 10000 b + 6600 b + 2673, 

750 t - 1850 b + 81, - 500 b - 9 + 200 s] 

4.4.3 Rational functions 

In Maple 

For comparison with other CAS choose from: Axiom (see chapter 4.1.3) Derive (see chapter 4.2.3) Macsyma (see 
chapter 4.3.3) Mathematica (see chapter 4.5.3) Reduce (see chapter 4.6.3) 

• integration 

int( (3*a*b~2 - 5*a~2* b) / (a~4 - 2), a ); 

2 1/2 
2 1/2 - 2 + a 2 3/4 3/4 

3/8 b 2 ln( ) - 5/4 b 2 arctan(l/2 a 2 ) 

2 1/2 
- 2 - a 2 

1/4 
3/4 a + 2 

+ 5/8 b 2 ln( ) 

1/4 
a - 2 



• verification by derivation 
diff(",a); 



/ 1/2 2 1/2 l/2\ 

2 1/2 1 a 2 (- 2 + a 2 ) a 2 | 2 1/2 

b 2 |2 + 2 1 (- 2 - a 2 ) 

12 1/2 2 1/2 2 | 

\ - 2 - a 2 (- 2 - a 2 ) / 

2 1/2 
- 2 + a 2 

1/2 
b 2 
- 5/4 

2 1/2 
1 + 1/2 a 2 

/ 1/4 \ 

3/4|l a + 2 | 1/4 

b 2 | | (a - 2 ) 

I 1/4 1/4 2| 

\a - 2 (a - 2 ) / 

+ 5/8 

1/4 
a + 2 



simplify ( " ); 



1/2 
(5 a - 3 b) b 2 a 



2 1/2 1/4 1/4 

(2 + a 2 ) (a - 2 ) (a + 2 ) 



simplify (numer (" ) /expand(denom(" )) ) ; 

(5 a - 3 b) b a 



4 
a - 2 



• partial fraction decomposition 

convert((10*x~2- ll*x- 6) / (x~3- x~2- 2*x),parfrac,x); 



5 2 

3/x + + 

x + 1 x - 2 



4.4.4 Solving equations 

In Maple 

For comparison with other CAS choose from: Axiom (see chapter 4.1.4) Derive (see chapter 4.2.4) Macsyma (see 
chapter 4.3.4) Mathematica (see chapter 4.5.4) Reduce (see chapter 4.6.4) 



Linear systems 

In Maple 

For comparison with other CAS choose from: Axiom (see chapter 4.1.4) Derive (see chapter 4.2.4) Macsyma (see 
chapter 4.3.4) Mathematica (see chapter 4.5.4) Reduce (see chapter 4.6.4) 

solve( 2*xl + x2 + 3*x3 - 9 = 0, xl - 2*x2 + x3 + 2 = 0, 3*xl + 2*x2 + 2*x3 -7 = 0, xl, x2, 
x3); 

{xl = -1, x2 = 2, x3 = 3} 



Nonlinear equations 

In Maple 

For comparison with other CAS choose from: Axiom (see chapter 4.1.4) Derive (see chapter 4.2.4) Macsyma (see 
chapter 4.3.4) Mathematica (see chapter 4.5.4) Reduce (see chapter 4.6.4) 

• using decomposition! for solving high degree polynomials 

solve( x~8 - 8*x~7 + 34*x~6 - 92*x~5 + 175*x~4 - 236*x~3 + 226*x~2 - 140*x + 46 ); 

1/2 1/2 1/2 1/2 

1+1/2 2 (- 3 + 7.2) ,1-1/22 (- 3 + 7.2) 

1/2 1/2 1/2 1/2 

1+1/2 2 (- 3 - 7.2) ,1-1/22 (- 3 - 7.2) 

1/2 1/2 1/2 1/2 

1+1/2 2 (- 3 + 7.D ,1-1/22 (- 3 + 7.D 

1/2 1/2 1/2 1/2 

1+1/2 2 (- 3 - 7.D ,1-1/22 (- 3 - 7.D 

1/2 1/2 

7.1 : = (-3-43) 

1/2 1/2 

7.2 : = (-3 + 43) 



• multiple use of inversion functions 

solve( log( arccos( arcsin(x~(2/3)-b)- 1)) + 2 ); 

3 2 

{x = \_XX1 , b = - sin(cos(exp(-2)) + 1) + \_XX1 } 



Nonlinear systems 

In Maple 

For comparison with other CAS choose from: Axiom (see chapter 4.1.4) Derive (see chapter ??) Macsyma (see 
chapter ??) Mathematica (see chapter 4.5.4) Reduce (see chapter 4.6.4) 

• polynomial systems by the use of (see chapter2.7) 

solve( alpha* cl - beta* cl~2 - gamma* cl*c2 + epsilon* c3 = 0, -gamma* cl*c2 + (ep- 
silon+theta)* c3 - eta* c2 = 0, gamma* cl*c2 + eta* c2 - (epsilon+theta)* c3 = ,c3,c2,cl); 



{cl = cl, 
c2 = 

cl (- epsilon alpha + cl epsilon beta - theta alpha + cl theta beta) 
- eta epsilon + cl theta gamma 



c3 



2 
cl (- alpha eta + beta cl eta - gamma cl alpha + gamma cl beta) 

- eta epsilon + cl theta gamma 



4.4.5 Analytical operations 

In Maple 

For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Derive (see chapter 4.2.5) Macsyma (see 
chapter 4.3.5) Mathematica (see chapter 4.5.5) Reduce (see chapter 4.6.5) 

Limits 

In Maple 

For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Derive (see chapter 4.2.5) Macsyma (see 
chapter 4.3.5) Mathematica (see chapter 4.5.5) Reduce (see chapter 4.6.5) 

limit ( sin(x) / x, x = ); 



limit( (3*sin(Pi*x)- sin(3*Pi*x) ) / x~3, x = ); 

3 
4 Pi 

limit ( (2*x + 5) / (3*x - 2), x = infinity ); 

2/3 



Taylor series 

In Maple 

For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Derive (see chapter 4.2.5) Macsyma (see 
chapter 4.3.5) Mathematica (see chapter 4.5.5) Reduce (see chapter 4.6.5) 

series( E~x, x=0, 5 ); 

2 3 4 5 

1 + x + 1/2 x + 1/6 x + 1/24 x + 0(x ) 



series( E~(x+y), x=0, 3); 



2 3 

exp(y) + exp(y) x + 1/2 exp(y) x + 0(x ) 



series( ", y=0, 3 ); 

0(1) 



Summation and Products 

In Maple 

For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Derive (see chapter 4.2.5) Macsyma (see 
chapter 4.3.5) Mathematica (see chapter 4.5.5) Reduce (see chapter 4.6.5) 

sum( m~2 * x~m, m=l..n ); 

(n + 1) 2 2 2 2 

x (- 2 (n + 1) x + 2 (n + 1) x + x + x + (n + 1) x 

2 2 / 3 x (x + 1) 

- 2 (n + 1) x + (n + 1) ) / (x - 1) 

/ 3 

(x - 1) 



sum( cos( (2*m-l)* Pi / (2*n+l) ), m=l..r ); 

Pi (r +1) Pi 2 Pi (r + 1) 

- cos( ) (- 2 cos( ) sin( ) 

2n+l 2n+l 2n+l 

Pi Pi (r + 1) Pi Pi (r + 1) 

+ 2 cos( ) cos( ) sin( ) + sin( )) 

2n+l 2n+l 2n+l 2n+l 

/ Pi Pi 

/ sin( ) + cos( ) 

/ 2n+l 2n+l 



simplify ( " ); 



Pi (r +1) Pi 2 Pi (r + 1) 

(2 cos( ) cos( ) sin( ) 

2n+l 2n+l 2n+l 

Pi Pi (r + 1) 2 Pi 

- 2 cos( ) cos( ) sin( ) 

2n+l 2n+l 2n+l 

Pi (r + 1) Pi (r + 1) Pi Pi 

- cos( ) sin( ) + cos( ) sin( )) 

2n+l 2n+l 2n+l 2n+l 



/ Pi 
/ sin( ) 

/ 2 n + 1 



product( exp( sin(m*x) ), m=l..n ); 

(exp( 

- sin((n + 1) x) cos(x) + sin((n + 1) x) + cos((n + 1) x) sin(x) 
1/2 

cos(x) - 1 

/ sin(x) 
)) / exp(l/2 ) 

/ cos(x) - 1 



Integration 

In Maple 

For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Derive (see chapter 4.2.5) Macsyma (see 
chapter 4.3.5) Mathematica (see chapter 4.5.5) Reduce (see chapter 4.6.5) 

int( x~2 * (a+b*x)~p, x); 

3 p 2 p 
a (a + b x) a p x (a + b x) 
2 2 

3 3 2 2 3 2 

b (p + 6 p + 11 p + 6) b (p + 6 p + 11 p + 6) 

2 p 3 p 

pax (a + b x) x (a + b x) 

+ + 

2 p + 3 

b (p + 5 p + 6) 



simplify ( " ); 

p222332 22 2 33 

(a + bx) (p ax b +x b p +pax b -2a pxb + 3x b p 

3 3 3/ 3 

+ 2 a + 2 x b ) / ((p + 3) (p + 2) (p + 1) b ) 

/ 



int( x~2* log(x~2 + a~2), x); 

3 2 2 2 3 3 
1/3 x ln(x + a ) + 2/3 ax- 2/9 x - 2/3 a arctan(x/a) 



int( x* d~x * sin(x), x ); 



x 2 2 x 
x d tan(l/2 x) (- 2 ln(d) + 2) d tan(l/2 x) 
( + 

2 2 2 

ln(d) + 1 (ln(d) + 1) 

x x 2 x 

ln(d) x d tan(l/2 x) ln(d) d tan(l/2 x) x d 

+ 2 2 

2 2 2 2 

ln(d) + 1 (ln(d) + 1) ln(d) + 1 

x 
ln(d) d / 2 

+ 2 ) / (1 + tan(l/2 x) ) 

2 2 / 

(ln(d) + 1) 



simplify ( " ); 



2 
(- 2 ln(d) cos(x) + x cos(x) ln(d) + x cos(x) - sin(x) 

3 2 x 

- ln(d) x sin(x) - ln(d) x sin(x) + sin(x) ln(d) ) d 

/ 4 2 

/ (ln(d) + 2 ln(d) + 1) 
/ 



int( x* sqrt(a + b*x)~p, x ); 

/ 



1/2 p 
x ((a + b x) ) dx 



int( 2*x* exp(x~2)* log(x)+ exp(x~2)/x + (log(x)-2)/(log(x)~2+x)~2 + ((2/x)* log(x)+ (l/x) + l)/ 
(log(x)~2+x), x ); 

2 ln(x) 2 

exp(x ) ln(x) + ln(ln(x) + x) 

2 
ln(x) + x 



Ordinary differential equations 

In Maple 

For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Derive (see chapter ??) Macsyma (see 
chapter 4.3.5) Mathematica (see chapter 4.5.5) Reduce (see chapter 4.6.5) 

dsolve( diff(y(x),x) + y(x)*sin(x)/cos(x) - l/cos(x) = 0, y(x) ); 



y(x) = sin(x) + cos(x) \_C1 



• Bernoulli equation 

dsolve( x*(l-x~2)* diff(y(x),x) + (2*x~2-l)* y(x) - x~3* y~3 = 0, y(x)); 

5 
1 2 x + 5 \_C1 
= 1/5 

2 2 2 

y(x) x (- 1 + x ) 



solve( ", y(x) ); 



1/2 5 7 2 1/2 

5 x (- 2 x +2x -5 \_C1 + 5 \_C1 x ) 

» 

5 
- 2 x - 5 \_C1 

1/2 5 7 2 1/2 

5 x(-2x+2x-5 \_C1 + 5 \_C1 x ) 

5 
- 2 x - 5 \_C1 



dsolve( diff(y(x),x,x) + 4* diff(y(x),x) + 4*y(x) - x* exp(x) = 0,y(x)); 

y(x) = 1/9 x exp(x) - 2/27 exp(x) + \_C1 exp(- 2 x) + \_C2 exp(- 2 x) x 

4.4.6 Matrices 

In Maple 

For comparison with other CAS choose from: Axiom (see chapter 4.1.6) Derive (see chapter 4.2.6) Macsyma (see 
chapter 4.3.6) Mathematica (see chapter 4.5.6) Reduce (see chapter 4.6.6) 

xx := array( [[all,al2], [a21,a22]] ); 

[ all al2 ] 

xx := [ ] 

[ a21 a22 ] 



yy := array( [yl,y2] ); 



yy := [ yl, y2 ] 



evalm( xx~(-l) * yy); 



a22 yl - al2 y2 a21 yl - all y2 
[ , ] 

all a22 - al2 a21 all a22 - al2 a21 



with(linalg); 

det( xx ); 



all a22 - al2 a21 



evalm( xx~(-2) ); 



a22 + al2 a21 



al2 (a22 + all) 



2 2 

(all a22 - al2 a21) (all a22 - al2 a21) 



a21 (a22 + all) 



al2 a21 + all 



2 2 

(all a22 - al2 a21) (all a22 - al2 a21) 



eigenvals( [[2,-1,1], [0,1,1], [-1,1,1]] ); 



2, 1, 1 



4.4.7 Graphics 

In Maple 

For comparison with other CAS choose from: Axiom (see chapter 4.1.7) Derive (see chapter 4.2.7) Macsyma (see 
chapter 4.3.8) Mathematica (see chapter 4.5.7) Reduce (see chapter 4.6.8) 

2D Graphics 

In Maple 

For comparison with other CAS choose from: Axiom (see chapter 4.1.7) Derive (see chapter 4.2.7) Macsyma (see 
chapter 4.3.8) Mathematica (see chapter 4.5.7) Reduce (see chapter 4.6.8) 

• function of one parameter 

• graph of several functions (Bessel functions J(n,x),n = 0,2,5) 



3D Graphics 

In Maple 

For comparison with other CAS choose from: Axiom (see chapter 4.1.7) Derive (see chapter 4.2.7) Macsyma (see 
chapter 4.3.8) Mathematica (see chapter 4.5.7) Reduce (see chapter 4.6.8) 

• graph of a function with 2 parameters 

• graph of another function with 2 parameters 



Parametric plots 

In Maple 

For comparison with other CAS choose from: Axiom (see chapter 4.1.7) Derive (see chapter 4.2.7) Macsyma (see 
chapter 4.3.8) Mathematica (see chapter 4.5.7) Reduce (see chapter 4.6.8) 

• parametricaly given surface 

Contour maps 

In Maple 

For comparison with other CAS choose from: Axiom (see chapter ??) Derive (see chapter ??) Macsyma (see 
chapter 4.3.8) Mathematica (see chapter 4.5.7) Reduce (see chapter 4.6.8) 

• map of a function with 2 parameters 

• map of another function with 2 parameters 

• map of further function with 2 parameters 

• map of transformation of cartesian coordinate grid by a complex function 

Polytopes 

In Maple 

For comparison with other CAS choose from: Axiom (see chapter ??) Derive (see chapter ??) Macsyma (see 
chapter ??) Mathematica (see chapter 4.5.7) Reduce (see chapter ??) 

• graphical presentation of polyhedrons 

4.4.8 Graphical presentation of formulas 

In Maple 

For comparison with other CAS choose from: Axiom (see chapter ??) Derive (see chapter ??) Macsyma (see 
chapter 4.3.9) Mathematica (see chapter 4.5.8) Reduce (see chapter 4.6.9) 

4.5 Mathematica 

• inputs 

outputs 

4.5.1 Number domains 

In Mathematica 

For comparison with other CAS choose from: Axiom (see chapter 4.1.1) Derive (see chapter 4.2.1) Macsyma (see 
chapter 4.3.1) Maple (see chapter 4.4.1) Reduce (see chapter 4.6.1) 

Big integers 

In Mathematica 

For comparison with other CAS choose from: Axiom (see chapter 4.1.1) Derive (see chapter 4.2.1) Macsyma (see 
chapter 4.3.1) Maple (see chapter 4.4.1) Reduce (see chapter 4.6.1) 

• integers of arbitrary size 
23~12 

21914624432020321 



60! 



8320987112741390144276341183223364380754172606361245952\ 
449277696409600000000000000 



Factorial[ 60 ] 

8320987112741390144276341183223364380754172606361245952\ 
449277696409600000000000000 



bi = 23"4*37*59*101 

61700183203 



• factorization of integers 
FactorInteger[ bi ] 

{{23, 4}, {37, 1}, {59, 1}, {101, 1}} 



bia = 23*11-6 

40745903 

• greatest common divisor of integers 
GCD[ bi, bia ] 

23 

Rational numbers 

In Mathematica 

For comparison with other CAS choose from: Axiom (see chapter 4.1.1) Derive (see chapter 4.2.1) Macsyma (see 
chapter 4.3.1) Maple (see chapter 4.4.1) Reduce (see chapter 4.6.1) 

• precise calculation with rational numbers 
1234567890/98765432 

617283945 
49382716 



rn = 1/24-2/15-64/47 

1027 
-( ) 

1410 



Complex numbers 

In Mathematica 

For comparison with other CAS choose from: Axiom (see chapter 4.1.1) Derive (see chapter 4.2.1) Macsyma (see 
chapter 4.3.1) Maple (see chapter 4.4.1) Reduce (see chapter 4.6.1) 

• precise calculation with complex numbers 
en = (2+3I)*(15-6I)+2/(2-4I) 

241 167 I 

5 5 



Big floating point numbers 

In Mathematica 

For comparison with other CAS choose from: Axiom (see chapter 4.1.1) Derive (see chapter 4.2.1) Macsyma (see 
chapter 4.3.1) Maple (see chapter 4.4.1) Reduce (see chapter 4.6.1) 

rn = - 1027/1410 

1027 
-( ) 

1410 



en = (167*1 + 241)/5 

241 167 I 

5 5 

• computation with floating point numbers with arbitrary number of digits 
N[ rn, 20 ] 

-0 . 728368794326241 13475 

N[ en, 10 ] 

48.2 + 33.4 I 



N[ Pi, 40 ] 

3 . 1415926535897932384626433832795028841972 

Cos[Pi] 

-1 

N[ Sin[l], 20 ] 

. 84147098480789650665 



• should be cos(pi/6) = sqrt(3)/2 
N[ Cos[Pi/6], 50 ] 

. 86602540378443864676372317075293618347140262690519 

%"2 

0.75 

• complex functions 
Tan[1.0 + 1.01] 

0.271753 + 1.08392 I 



Log[N[l,30]+I] 

. 346573590279972654708616060729 + 
. 78539816339744830961566084582 I 



Precision[%] 

30 



N[ Log[N[l,40] + I], 60 



. 3465735902799726547086160607290882840378 + 
. 7853981633974483096156608458198757210493 I 



Precision[%] 

40 



4.5.2 Polynomials 

In Mathematica 

For comparison with other CAS choose from: Axiom (see chapter 4.1.2) Derive (see chapter 4.2.2) Macsyma (see 
chapter 4.3.2) Maple (see chapter 4.4.2) Reduce (see chapter 4.6.2) 

Basic operations 

In Mathematica 

For comparison with other CAS choose from: Axiom (see chapter 4.1.2) Derive (see chapter 4.2.2) Macsyma (see 
chapter 4.3.2) Maple (see chapter 4.4.2) Reduce (see chapter 4.6.2) 

• by default parantheses are not expanded 
pol := (a+b+c)~4 
Expand [ pol ] 



4 3 2 2 3 4 3 

a +4a b + 6a b +4ab +b +4a c + 

2 2 3 2 2 2 

12 a bc + 12ab c + 4b c + 6a c +12abc + 

2 2 3 3 4 

6b c +4ac +4bc +c 



• derivation 

dpol = D[ Expand[pol], a ] 

3 2 2 3 2 

4a +12 a b + 12 a b +4b +12 a c+24abc+ 

2 2 2 3 

12 b c + 12 a c +12 be +4c 

D[ pol, a, b, 2 ] 

24 (a + b + c) 

• integration 
Integrate [ dpol, a ] 

4 3 2 2 3 

a + 4 a (b + c) + 6 a (b + c) + 4 a (b + c) 



• verification 
% - pol 

4 3 2 2 3 

a + 4 a (b + c) + 6 a (b + c) + 4 a (b + c) 

4 
(a + b + c) 



Expand [%] 



4 3 2 2 3 4 

-b -4b c-6b c -4bc -c 



• greatest common divisor of polynomials 
(a~2-b~2)/(a~2-2 a b+b~2) 

2 2 

a - b 



2 2 

a - 2 a b + b 



Simplify[%] 

a + b 
a - b 

g := 34 x~19 - 25 x~16 + 70 x~7 + 20 x~3 - 91 x - 86 

fl := g (64 x~34 - 21 x~47 - 126 x~8 - 46 x~5 - 16 x~60 - 81) 

fl 

3 7 16 19 
(-86 - 91 x + 20 x + 70 x - 25 x + 34 x ) 

5 8 34 47 60 
(-81 - 46 x - 126 x + 64 x - 21 x - 16 x ) 



fl = Expand[%] 

3 5 6 7 
6966 + 7371 x - 1620 x + 3956 x + 4186 x - 5670 x + 

8 9 11 12 15 
9916 x + 11466 x - 2520 x - 3220 x - 8820 x + 

16 19 21 24 27 
2025 x - 2754 x + 1150 x + 1586 x - 4284 x 

34 35 37 41 47 
5504 x - 5824 x + 1280 x + 4480 x + 1806 x + 

48 50 53 54 60 
1911 x - 2020 x + 2176 x - 1470 x + 1376 x + 

61 63 66 67 76 
1456 x + 205 x - 714 x - 1120 x + 400 x 

79 
544 x 



f2 := g (72 x"60 - 25 x~25 - 19 x~23 -22 x~39 - 83 x~52 + 54 x~10 + 81) 
f2 = Expand [f2] 

3 7 10 
-6966 - 7371 x + 1620 x + 5670 x - 4644 x 

11 13 16 17 19 
4914 x + 1080 x - 2025 x + 3780 x + 2754 x + 

23 24 25 26 28 
1634 x + 1729 x + 2150 x + 545 x - 500 x + 

29 30 32 39 40 
1836 x - 1330 x - 1750 x + 2367 x + 2002 x + 



41 42 44 46 52 
625 x - 1086 x - 850 x - 1540 x + 7138 x + 

53 55 58 59 60 
7553 x - 1110 x - 748 x - 5810 x - 6192 x 

61 63 67 68 71 
6552 x + 1440 x + 5040 x + 2075 x - 2822 x 

76 79 
1800 x + 2448 x 



PolynomialGCD[ fl, f2 ] 

3 7 16 19 

-86 - 91 x + 20 x + 70 x - 25 x + 34 x 



Factorization 

In Mathematica 

For comparison with other CAS choose from: Axiom (see chapter 4.1.2) Derive (see chapter 4.2.2) Macsyma (see 
chapter 4.3.2) Maple (see chapter 4.4.2) Reduce (see chapter 4.6.2) 

• factorization is transformation of a polynomial into a product of polynomials 
Factor[ a~2 - b~2 ] 

(a - b) (a + b) 



Factor[ a~2 + b~2 ] 



2 2 

a + b 



Factor [a" 2 + b~2, Gaussianlntegers -£ True ] 

(-1 a + b) (I a + b) 

fa := ( x~2 z + y~4 z~2 + 5) * ( x y~3 + z~2) * (-x~3 y + z~2 + 3) * ( x~3 y~4 + z~2 ) 
fa = Expand [fa] 

47 78 67 98 32 

15 x y -5x y +3x y z-x y z+15xy z + 

342 442 652 472 

15xyz -5xyz-5xyz+5xyz + 

4 11 2 7 12 2 333 543 

3xy z-xy z+3xyz+3xyz- 

643 853 673 4 3 4 

xyz -xyz +xyz+15z -5xyz + 



34 344 74 384 

5xy z +5x y z +3xy z +3x y z 

484 694 4 11 4 2 5 5 5 

xyz-xyz+xy z+3xz-xyz + 

335 545 6 46 356 

xyz+xyz+5z+3yz-xyz + 

76 386 27 48 

xy z +x y z +x z +y z 



Fact or [ fa 



3 2 3 2 3 4 2 

(3-x y + z)(xy +z)(x y +z) 

2 4 2 

(5 + x z + y z) 



Decomposition 

In Mathematica 

For comparison with other CAS choose from: Axiom (see chapter ??) Derive (see chapter ??) Macsyma (see 
chapter 4.3.2) Maple (see chapter 4.4.2) Reduce (see chapter 4.6.2) 

• finding of polynomial substitution which transforms a polynomial into another polynomial in further polyno- 
mials 



Decompose[ x ~ 6 + 9x ~ 5 + 52x~4 + 177x~3 + 435x~2 + 630x + 593 ] 

2 3 
{593 + 210 x + 25 x + x , x (3 + x)} 



Decompose[ x ~ 4 + 2x ~ 3 Y + 3x ~ 2 y~ 2 + 2x Y~ 3 + y~ 4 + 2x ~ 2 Y + 2x y~ 2 + 2 y~ 3 + 5x ~ 2 + 5 x 
y + 6y~2 + 5y + 9, x ] 

2 2 2 3 

{9 + 5 x + x +5y+2xy+6y +2xy +2y + 

4 
y , x (x + y)} 



Grobner bases 

In Mathematica 

For comparison with other CAS choose from: Axiom (see chapter 4.1.2) Derive (see chapter ??) Macsyma (see 
chapter 4.3.2) Maple (see chapter 4.4.2) Reduce (see chapter 4.6.2) 

polys := 45p + 35s - 165b - 36, 35p + 40z + 25t - 27s, 15w + 25p s + 30z - 18t - 165b~2, -9w + 
15p t + 20z s, w p + 2z t - llb~3, 99w - lis b + 3b~2, b~2 + 33/50b + 2673/10000 

vars := w,p,z,t,s,b 

GroebnerBasis[ polys, vars ] 



2 
{2673 + 6600 b + 10000 b , -9 - 500 b + 200 s, 

81 - 1850 b + 750 t, 10287 + 24500 b + 18000 z, 

-1377 - 3100 b + 1800 p, 3969 + 9500 b + 60000 w} 



4.5.3 Rational functions 

In Mathematica 

For comparison with other CAS choose from: Axiom (see chapter 4.1.3) Derive (see chapter 4.2.3) Macsyma (see 
chapter 4.3.3) Maple (see chapter 4.4.3) Reduce (see chapter 4.6.3) 

• integration 

Integrate[ ( 3a b ~ 2 " 5a ~ 2 b )/( a ~ 4 " 2 )> a ] 

a 

-5 b ArcTan[ ] 

1/4 
2 

1/4 
2 2 

3/4 2 3/4 

(10 b - 3 2 b ) Log [-2 +2 a] 

1/4 
8 2 

3/4 2 3/4 

(-10 b - 3 2 b ) Log [2 +2 a] 

1/4 
8 2 

2 2 

3 b Log [2 + Sqrt[2] a ] 

4 Sqrt[2] 



• verification by derivation 
D[%,a] 

2 
-5b 3 a b 

2 2 

a 2 (2 + Sqrt[2] a ) 

2 Sqrt[2] (1 + ) 

Sqrt [2] 

3/4 2 3/4 2 



-10b-32 b 10b-32 b 

3/4 3/4 

4 Sqrt[2] (2+2 a) 4 Sqrt [2] (-2 +2 a) 

Simplify[%] 

a b (-5 a + 3 b) 

4 
-2 + a 



• partial fraction decomposition 

Apart[ (10x~2- llx- 6)/(x~3- x~2- 2x)] 

2 3 5 

-2 + x x 1 + x 



4.5.4 Solving equations 

In Mathematica 

For comparison with other CAS choose from: Axiom (see chapter 4.1.4) Derive (see chapter 4.2.4) Macsyma (see 
chapter 4.3.4) Maple (see chapter 4.4.4) Reduce (see chapter 4.6.4) 

Linear systems 

In Mathematica 

For comparison with other CAS choose from: Axiom (see chapter 4.1.4) Derive (see chapter 4.2.4) Macsyma (see 
chapter 4.3.4) Maple (see chapter 4.4.4) Reduce (see chapter 4.6.4) 

Solve[2 xl + x2 + 3 x3 - 9 == 0, xl - 2 x2 + x3 + 2 == 0, 3 xl 4- 2 x2 4- 2 x3 - 7 == 0, xl, x2, 
x3 ] 

{{xl -> -1, x2 -> 2, x3 -> 3}} 



Nonlinear equations 

In Mathematica 

For comparison with other CAS choose from: Axiom (see chapter 4.1.4) Derive (see chapter 4.2.4) Macsyma (see 
chapter 4.3.4) Maple (see chapter 4.4.4) Reduce (see chapter 4.6.4) 

• using decomposition! for solving high degree polynomials 

Solve[ x~8 - 8x~7 + 34x~6 - 92x~5 + 175x~4 - 236x~3 4- 226x~2 - 140x 4- 46 == ] 

4 - Sqrt [16 - 8 (5 - Sqrt [-3 - 4 Sqrt [3]])] 
{{x -> }, 

4 

4 + Sqrt [16 - 8 (5 - Sqrt [-3 - 4 Sqrt [3]])] 
{x _> }) 

4 
4 - Sqrt [16 - 8 (5 + Sqrt [-3 - 4 Sqrt [3]])] 



{x _> }f 

4 

4 + Sqrt[16 - 8 (5 + Sqrt [-3 - 4 Sqrt[3]])] 
{x _> }) 

4 

4 - Sqrt [16 - 8 (5 - Sqrt [-3 + 4 Sqrt [3]])] 
{x _> }) 

4 

4 + Sqrt [16 - 8 (5 - Sqrt [-3 + 4 Sqrt [3]])] 
{x _> }) 

4 

4 - Sqrt [16 - 8 (5 + Sqrt [-3 + 4 Sqrt [3]])] 
{x _> }) 

4 

4 + Sqrt [16 - 8 (5 + Sqrt [-3 + 4 Sqrt [3]])] 
{x _> }} 

4 



• multiple use of inversion functions 

Solve[ Log[ ArcCos[ ArcSin[x~(2/3)-b]- 1]]+ 2 == 0, x ] 

3 2 -2 

{{x -> -Sqrt[b + 3 b Sin[l + Cos [E ]] + 

-2 2 -2 3 

3 b Sin[l + Cos[E ]] + Sin[l + Cos [E ]] ]}, 

3 2 -2 

{x -> Sqrt[b + 3 b Sin[l + Cos[E ]] + 

-2 2 -2 3 

3 b Sin[l + Cos[E ]] + Sin[l + Cos [E ]] ]}} 



Simplify[%] 



12 3 

{{x -> -Sqrt[(b + Sin[2 Cos [ ] ]) ]}, 

2 
2 E 

12 3 

{x -> Sqrt[(b + Sin[2 Cos [ ] ]) ]}} 

2 
2 E 



Nonlinear systems 
In Mathematica 



For comparison with other CAS choose from: Axiom (see chapter 4.1.4) Derive (see chapter ??) Macsyma (see 
chapter ??) Maple (see chapter 4.4.4) Reduce (see chapter 4.6.4) 

• polynomial systems by the use of (see chapter2.7) 

Solve[ alpha cl -beta cl~2 -gamma cl c2 +epsilon c3 == 0, -gamma cl c2 + (epsilon+theta) c3 
- eta c2 == 0, gamma cl c2 + eta c2 - (epsilon+theta) c3 == , c3,c2,cl ] 

{{c3 -> (cl (-alpha + beta cl - 

alpha cl epsilon gamma 

+ 

epsilon eta - cl gamma theta 

2 
beta cl epsilon gamma 

epsilon eta - cl gamma theta 

alpha cl gamma theta 

+ 

epsilon eta - cl gamma theta 

2 
beta cl gamma theta 
)) / epsilon, 

epsilon eta - cl gamma theta 
c2 -> (cl (alpha epsilon - beta cl epsilon + 
alpha theta - beta cl theta)) / 
(-(epsilon eta) + cl gamma theta)}} 

Simplify[%] 

cl (-alpha + beta cl) (eta + cl gamma) 
{{c3 -> , 

epsilon eta - cl gamma theta 

cl (alpha - beta cl) (epsilon + theta) 
c2 -> }} 

-(epsilon eta) + cl gamma theta 

• cl is an arbitrary complex number 

4.5.5 Analytical operations 

In Mathematica 

For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Derive (see chapter 4.2.5) Macsyma (see 
chapter 4.3.5) Maple (see chapter 4.4.5) Reduce (see chapter 4.6.5) 

Limits 

In Mathematica 



For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Derive (see chapter 4.2.5) Macsyma (see 
chapter 4.3.5) Maple (see chapter 4.4.5) Reduce (see chapter 4.6.5) 
Limit [ Sin[x]/x, x-j,0] 



Limit[(3 Sin[Pi x] - Sin[3 Pi x])/x~3, x-iO] 

3 
4 Pi 

Limit[ (2x + 5)/ (3x - 2), x - L Infinity ] 

2 
3 



Taylor series 

In Mathematica 

For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Derive (see chapter 4.2.5) Macsyma (see 
chapter 4.3.5) Maple (see chapter 4.4.5) Reduce (see chapter 4.6.5) 

Series[ E~x, x,0,4 ] 

2 3 4 

xxx 5 

1 + x + — + — + — + 0[x] 
2 6 24 



Series[ E~(x+y), x,0,2, y,0,2 ] 

2 2 

y 3 y 3 

1 + y + — + 0[y] + (1 + y + -- + 0[y] ) x + 
2 2 

2 

1 y y 3 2 3 
(- + - + — + 0[y] ) x + 0[x] 

2 2 4 



%"2 

2 3 2 3 

l + 2y + 2y + [y] +(2 + 4y + 4y +0[y])x + 

2 3 2 3 

(2 + 4 y + 4 y + [y] ) x + [x] 



Summation and Products 

In Mathematica 

For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Derive (see chapter 4.2.5) Macsyma (see 
chapter 4.3.5) Maple (see chapter 4.4.5) Reduce (see chapter 4.6.5) 

• modul Algebra' SymbolicSum' dealing with sumation 
Needs[ " Algebra'SymbolicSum'" ] 
Sum[ m~2 x~m, m, 1, n ] 

n n 2 n 1 + n 

(x (-1 -x+x+2nx+nx+x 

1 + n 21+n 22 + n 3 

2nx -2nx +nx ))/ (-1 + x) 



Sum[ Cos[(2m-l) Pi/ (2n+l)], m,l,r ] 

Pi Pi Pi r 

-Cos[ ] + Cos[ ] 

l + 2n l + 2n l + 2n 

Pi Pi Pi r 

Csc[ ] Sin[ + ] 

l + 2n l + 2n l+2n 



Product[ Exp[ Sin[m x]], m, 1, n ] 

Csc[x/2] Sin[(n x)/2] Sin[((l + n) x)/2] 



Sum[ m 2~m/ Factorial[m+2], m,l,n] 

n n 

-6 2 - 2 2 n + Gamma [4 + n] 



Gamma [4 + n] 



Integration 

In Mathematica 

For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Derive (see chapter 4.2.5) Macsyma (see 
chapter 4.3.5) Maple (see chapter 4.4.5) Reduce (see chapter 4.6.5) 

Integrate [ x~2 (a+b x)~p, x] 

3 
p 2 a 

(a + b x) ( 

3 
b (1 + p) (2 + p) (3 + p) 



2 

2 a p x 

+ 

(3 + p) (b + b p) (2 b + b p) 

2 3 
a p x b x 
+ ) 

(3 + p) (2 b + b p) 3 b + b p 



Simplify[%] 



1 + p 2 2 2 

((a +bx) (2a -2abx-2abpx+2b x + 

22222 3 23 

3b px +b p x )) / (b (6 + 11 p + 6 p + p )) 



Integrate[ x~2 Log[x~2 + a~2], x] 

3 x 

2 3 2a ArcTan[-] 3 22 

2 a x 2 x ax Log [a + x ] 

+ 

3 9 3 3 



Simplify[%] 



2 3 3 x 3 2 2 

6a x-2x -6a ArcTan[-] + 3 x Log [a + x ] 

a 



Simplify[Integrate[ x d~x Sin[x], x ]] 

x 2 

(d (-(x Cos[x]) + 2 Cos[x] Log[d] - x Cos [x] Log[d] + 

2 

Sin[x] + x Log[d] Sin[x] - Log[d] Sin[x] + 

3 2 2 

x Log[d] Sin[x])) / (1 + Log[d] ) 



Simplify[Integrate[ x Sqrt[a + b x]~p, x]] 

1 + p/2 
2 (a + b x) (-2 a + 2 b x + b p x) 



2 

b (2 + p) (4 + p) 



Simplify[ Integrate[ 2x Exp[x~2] Log[x]+ Exp[x~2]/x + (Log[x]-2)/(Log[x]~2+x)~2 + ((2/x) Log[x] + 
(l/x) + l)/ (Log[x]~2+x), x]] 

2 
x Log [x] 2 

E Log[x] + Log[x + Log[x] ] 

2 
x + Log[x] 



Ordinary differential equations 

In Mathematica 

For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Derive (see chapter ??) Macsyma (see 
chapter 4.3.5) Maple (see chapter 4.4.5) Reduce (see chapter 4.6.5) 

DSolve[ y'[x] + y[x] Sin[x]/Cos[x] - 1/Cos[x]==0, y[x],x ] 

{{y[i] -> C[l] Cos[x] + Sin[x]}} 



• Bernoulli equation 

Simplify[ DSolve[ x (l-x~2) y'[x] + (2x~2-l) y[x]- x~3 y~3 == 0, y[x], x] ] 

3 2 

{{y[x] -> x (y + Sqrt[l - x ] C[l])}} 



Simplify[ DSolve[ y" [x] + 4y'[x] + 4y[x] - x Exp[x] == 0, y[x], x ] ] 

X X 

-2 E Ex C[l] x C[2] 

«y[x] -> + + + }} 

27 9 2 x 2 x 

E E 



4.5.6 Matrices 

In Mathematica 

For comparison with other CAS choose from: Axiom (see chapter 4.1.6) Derive (see chapter 4.2.6) Macsyma (see 
chapter 4.3.6) Maple (see chapter 4.4.6) Reduce (see chapter 4.6.6) 

xx := all,al2, a21,a22 

yy := yi,y2 

xx // MatrixForm 

all al2 

a21 a22 



Det[ xx ] 



-(al2 a21) + all a22 

zz = Inverse[xx] . yy // MatrixForm 

a22 yl al2 y2 

-(al2 a21) + all a22 -(al2 a21) + all a22 

a21 yl all y2 
_ ( ) + 

-(al2 a21) + all a22 -(al2 a21) + all a22 

Simplify[%] // MatrixForm 

a22 yl - al2 y2 

-(al2 a21) + all a22 

-(a21 yl) + all y2 
-(al2 a21) + all a22 

MatrixPower[ xx, -2 ] // Simplify // MatrixForm 

2 

al2 a21 + a22 al2 (all + a22) 
_ ( ) 

2 2 

(-(al2 a21) + all a22) (al2 a21 - all a22) 

2 

a21 (all + a22) all + al2 a21 
_ ( ) 

2 2 

(al2 a21 - all a22) (-(al2 a21) + all a22) 

Eigenvalues! VM, 0,1,1, -1,1,1 ] 

{1, 1, 2} 



4.5.7 Graphics 

In Mathematica 

For comparison with other CAS choose from: Axiom (see chapter 4.1.7) Derive (see chapter 4.2.7) Macsyma (see 
chapter 4.3.8) Maple (see chapter 4.4.7) Reduce (see chapter 4.6.8) 

2D Graphics 

In Mathematica 

For comparison with other CAS choose from: Axiom (see chapter 4.1.7) Derive (see chapter 4.2.7) Macsyma (see 
chapter 4.3.8) Maple (see chapter 4.4.7) Reduce (see chapter 4.6.8) 



• function of one parameter 

• including labels 

• graph of several functions (Bessel functions J(n,x),n = 0,2,5) 

3D Graphics 

In Mathematica 

For comparison with other CAS choose from: Axiom (see chapter 4.1.7) Derive (see chapter 4.2.7) Macsyma (see 
chapter 4.3.8) Maple (see chapter 4.4.7) Reduce (see chapter 4.6.8) 

• graph of a function with 2 parameters 

• graph of another function with 2 parameters 

Parametric plots 

In Mathematica 

For comparison with other CAS choose from: Axiom (see chapter 4.1.7) Derive (see chapter 4.2.7) Macsyma (see 
chapter 4.3.8) Maple (see chapter 4.4.7) Reduce (see chapter 4.6.8) 

• parametricaly given surface 

• combination of graphical objects 

Contour maps 

In Mathematica 

For comparison with other CAS choose from: Axiom (see chapter ??) Derive (see chapter ??) Macsyma (see 
chapter 4.3.8) Maple (see chapter 4.4.7) Reduce (see chapter 4.6.8) 

• map of a function with 2 parameters 

• map of another function with 2 parameters 

• map of further function with 2 parameters 

• map of transformation of cartesian coordinate grid by a complex function 

Polytopes 

In Mathematica 

For comparison with other CAS choose from: Axiom (see chapter ??) Derive (see chapter ??) Macsyma (see 
chapter ??) Maple (see chapter 4.4.7) Reduce (see chapter ??) 

• graphical presentation of polytopes 

4.5.8 Graphical presentation of formulas 

In Mathematica 

For comparison with other CAS choose from: Axiom (see chapter ??) Derive (see chapter ??) Macsyma (see 
chapter 4.3.9) Maple (see chapter 4.4.8) Reduce (see chapter 4.6.9) 

4.6 Reduce 

• inputs 

outputs 



4.6.1 Number domains 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter 4.1.1) Derive (see chapter 4.2.1) Macsyma (see 
chapter 4.3.1) Maple (see chapter 4.4.1) Mathematica (see chapter 4.5.1) 

Big integers 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter 4.1.1) Derive (see chapter 4.2.1) Macsyma (see 
chapter 4.3.1) Maple (see chapter 4.4.1) Mathematica (see chapter 4.5.1) 

• integers of arbitrary size 
23-12; 

21914624432020321 



factorial 60; 

8320987112741390144276341183223364380754172606361245952449277696409600 
000000000000 



bi:=23~4*37*59*101; 

bi := 61700183203 



• factorization of integers 
on ifactor; 
factorize bi; 

{23,23,23,23,37,59,101} 



bia:=23*ll~6; 

bia := 40745903 



• integer greatest common divisor 
gcd(bi,bia); 

23 



Rational numbers 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter 4.1.1) Derive (see chapter 4.2.1) Macsyma (see 
chapter 4.3.1) Maple (see chapter 4.4.1) Mathematica (see chapter 4.5.1) 

• exact calculation with rational numbers 
1234567890/98765432; 

617283945 

49382716 



rn:=l/2+2/15-64/47; 

- 1027 

rn := 

1410 

Complex numbers 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter 4.1.1) Derive (see chapter 4.2.1) Macsyma (see 
chapter 4.3.1) Maple (see chapter 4.4.1) Mathematica (see chapter 4.5.1) 

• exact calculation with complex numbers 
cn: = (2+3*i)*(15-6*i)+2/(2-4*i); 

63*i - 115 

en := 

2*i - 1 



• basic number domain changed to that of complex numbers 
on complex; 
cn:=cn; 

241 + 167*i 

en := 

5 



Algebraic numbers 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter 4.1.1) Derive (see chapter 4.2.1) Macsyma (see 
chapter 4.3.1) Maple (see chapter ??) Mathematica (see chapter ??) 

• module arnum dealing with algebraic numbers 
load_package arnum; 
defpoly sqrt2**2-2; 
l/(sqrt2+l); 



sqrt2 - 1 

(x**2+2*sqrt2*x+2)/(x+sqrt2); 
x + sqrt2 



on gcd; 

(x**3+(sqrt2-2)*x**2-(2*sqrt2+3)*x-3*sqrt2)/(x**2-2); 

2 
x - 2*x - 3 

x - sqrt2 



off gcd; 

sqrt(x**2-2*sqrt2*x*y+2*y**2); 

x - sqrt2*y 



• multiple algebraic extensions 
off arnum; 
defpoly sqrt5**2-5,cbrt3**3-3; 

*** Defining the polynomial for a primitive element: 

6 4 3 2 

al - 15*al - 6*al + 75*al - 90*al - 116 



cbrt3**3; 

3 

sqrt5**2; 
5 



cbrt3: 



120 5 27 4 2000 3 1170 2 6676 

( *al + *al *al *al + *al 

8243 8243 8243 8243 8243 

6825 
) 

8243 



S qrt(x**2+2*(sqrt5-cbrt3)*x+5-2*sqrt5*cbrt3+cbrt3**2); 

240 5 54 4 4000 3 2340 2 21595 

x + ( *al + *al *al *al + *al 

8243 8243 8243 8243 8243 

13650 
) 

8243 



Big floating point numbers 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter 4.1.1) Derive (see chapter 4.2.1) Macsyma (see 
chapter 4.3.1) Maple (see chapter 4.4.1) Mathematica (see chapter 4.5.1) 

rn:= - 1027/1410; 

- 1027 

rn := 

1410 

cn: = (167*i + 241)/5; 

167*i + 241 

en := 

5 



• computation with floating point numbers 
on rounded: 



rn: 



0.728368794326 



33.4*i + 48.2 

pi; 

3.14159265359 



cos 


pi: 


) 






1 


sin 


i; 






0. 


841470984808 



• computation with an arbitrary number of digits 
precision 50; 

12 



pi; 

3 . 1415926535897932384626433832795028841971693993751 

cos pi; 

- 1 



• should be cos(pi/6) = sqrt(3)/2 
cos(pi/6); 

. 86602540378443864676372317075293618347140262690519 



ws**2; 
0.75 

precision 10; 

50 



• by default, an underflow is converted to zero 
exp(-100000.1**2); 



• it is possible to get very small numbers 
on roundbf; 
exp(-100000.1**2); 

1 . 184683941e-4342953505 



off roundbf; 

• complex functions 
on complex; 



*** Domain mode rounded changed to complex-rounded 

tan(1.0 + 1.0*i); 

0.2717525853 + 1 .083923327*i 

log(1.0 + 1.0*i); 

0.3465735903 + 0.7853981634*i 



4.6.2 Polynomials 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter 4.1.2) Derive (see chapter 4.2.2) Macsyma (see 
chapter 4.3.2) Maple (see chapter 4.4.2) Mathematica (see chapter 4.5.2) 

Basic operations 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter 4.1.2) Derive (see chapter 4.2.2) Macsyma (see 
chapter 4.3.2) Maple (see chapter 4.4.2) Mathematica (see chapter 4.5.2) 

• by default, parentheses are expanded 
pol: = (a+b+c)~4; 

43 3 22 2 22 3 

pol := a + 4*a *b + 4*a *c + 6*a *b + 12*a *b*c + 6*a *c + 4*a*b 

2 2 3 4 3 2 2 

+ 12*a*b *c + 12*a*b*c + 4*a*c + b + 4*b *c + 6*b *c 

3 4 

+ 4*b*c + c 

• differentiation 
dpol:=df(pol,a); 

3 2 2 2 2 3 

dpol := 4* (a + 3*a *b + 3*a *c + 3*a*b + 6*a*b*c + 3*a*c + b 

2 2 3 

+ 3*b *c + 3*b*c + c ) 



df(pol,a,b,2); 

24* (a + b + c) 

• integration 
int(dpol,a); 



3 2 2 2 2 3 2 

a* (a + 4*a *b + 4*a *c + 6*a*b + 12*a*b*c + 6*a*c + 4*b + 12*b *c 

2 3 

+ 12*b*c + 4*c ) 



• verification 
ws-pol; 



4 3 2 2 3 4 

b - 4*b *c - 6*b *c - 4*b*c - c 



• polynomial greatest common divisor 
on gcd,ezgcd; 
(a~2-b~2)/(a~2-2*a*b+b~2); 

a + b 

a - b 



off gcd; 

g := 34*x**19-91*x+70*x**7-25*x**16+20*x**3-86; 

19 16 7 3 

s : = 34*x - 25*x + 70*x + 20*x - 91*x - 86 



fl:=g * (64*x**34-21*x**47-126*x**8-46*x**5-16*x**60-81); 

79 76 67 66 63 61 
fl := - 544*x + 400*x - 1120*x - 714*x + 205*x + 1456*x 

60 54 53 50 48 
+ 1376*x - 1470*x + 2176*x - 2020*x + 1911*x 

47 41 37 35 34 
+ 1806*x + 4480*x + 1280*x - 5824*x - 5504*x 

27 24 21 19 16 

- 4284*x + 1586*x + 1150*x - 2754*x + 2025*x 

15 12 11 9 8 

- 8820*x - 3220*x - 2520*x + 11466*x + 9916*x 

7 6 5 3 

- 5670*x + 4186*x + 3956*x - 1620*x + 7371*x + 6966 



f2:=g * (72*x**60-25*x**25-19*x**23-22*x**39-83*x**52+54*x**10+81); 



79 76 71 68 67 63 
f2 := 2448*x - 1800*x - 2822*x + 2075*x + 5040*x + 1440*x 

61 60 59 58 55 

- 6552*x - 6192*x - 5810*x - 748*x - 1110*x 

53 52 46 44 42 

+ 7553*x + 7138*x - 1540*x - 850*x - 1086*x 

41 40 39 32 30 
+ 625*x + 2002*x + 2367*x - 1750*x - 1330*x 

29 28 26 25 24 

+ 1836*x - 500*x + 545*x + 2150*x + 1729*x 

23 19 17 16 13 

+ 1634*x + 2754*x + 3780*x - 2025*x + 1080*x 

11 10 7 3 

- 4914*x - 4644*x + 5670*x + 1620*x - 7371*x - 6966 



gcd(fl,f2); 



19 16 7 3 
34*x - 25*x + 70*x + 20*x - 91*x - 86 



Factorization 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter 4.1.2) Derive (see chapter 4.2.2) Macsyma (see 
chapter 4.3.2) Maple (see chapter 4.4.2) Mathematica (see chapter 4.5.2) 

• factorization is the transformation of a polynomial into a product of polynomials 
factorize(a~2-b~2); 



{a - b,a + b} 



on complex; 
factorize(a~2+b~2); 

{a - i*b,a + i*b} 



off complex; 

fa:=(x**2*z+y**4*z**2+5)* (x*y**3+z**2)* (-x**3*y+z**2+3)* (x**3*y**4+z**2); 

98 853 7 12 2 78 694 

fa := - x *y *z - x *y *z - x *y *z - 5*x *y - x *y *z 

673 67 652 643 545 

+ x *y *z + 3*x *y *z - 5*x *y *z - x *y *z + x *y *z 



543 5 5 4 11 4 4 11 2 484 
+ 3*x *y *z - x *y*z + x *y *z + 3*x *y *z - x *y *z 

472 47 442 386 384 
+ 5*x *y *z + 15*x *y - 5*x *y *z + x *y *z + 3*x *y *z 

356 344 342 335 333 

- x *y *z + 5*x *y *z + 15*x *y *z + x *y *z + 3*x *y *z 

3427 25 76 74 

- 5*x *y*z + x *z + 3*x *z + x*y *z + 3*x*y *z 

34 3248 46 6 4 

+ 5*x*y *z + 15*x*y *z + y *z + 3*y *z + 5*z + 15*z 



factorize fa; 

3 2 

{x*y + z , 

2 4 2 

x *z + y *z +5, 

3 2 

- x *y + z +3, 

3 4 2 

x *y + z } 



Decomposition 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter ??) Derive (see chapter ??) Macsyma (see 
chapter 4.3.2) Maple (see chapter 4.4.2) Mathematica (see chapter 4.5.2) 

• decomposing a polynomial into simpler polynomials which when composed, produce the original polynomial 
decompose(x**6+9x**5+52x**4+177x**3+435x**2+630x+593); 

3 2 2 

{u + 25*u + 210*u + 593, u=x + 3*x} 



decompose(x**4+2x**3*y + 3x**2*y**2 + 2x*y**3 + y**4 + 2x**2*y +2x*y**2 + 2y**3 + 
5 x**2 + 5*x*y + 6*y**2 + 5y + 9); 

2 2 2 

{u + 5*u + 9,u=x + x*y + y + y} 

Simplification with side relations 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter ??) Derive (see chapter ??) Macsyma (see 
chapter ??) Maple (see chapter ??) Mathematica (see chapter ??) 



• load the module compact for polynomial compactification: simplification of a polynomial in which the variables 
are restricted by side relations 

load_package compact; 

com: = (l-(sin x)**2)**5*(l-(cos x)**2)**5*((sin x)**2+(cos x)**2)**5; 

20 10 20 8 20 6 

com := sin(x) *cos(x) - 5*sin(x) *cos(x) + 10*sin(x) *cos(x) 

20 4 20 2 20 

- 10*sin(x) *cos(x) + 5*sin(x) *cos(x) - sin(x) 

18 12 18 10 

+ 5*sin(x) *cos(x) - 30*sin(x) *cos(x) 

18 8 18 6 

+ 75*sin(x) *cos(x) - 100*sin(x) *cos(x) 

18 4 18 2 18 

+ 75*sin(x) *cos(x) - 30*sin(x) *cos(x) + 5*sin(x) 

16 14 16 12 

+ 10*sin(x) *cos(x) - 75*sin(x) *cos(x) 

16 10 16 8 

+ 235*sin(x) *cos(x) - 400*sin(x) *cos(x) 

16 6 16 4 

+ 400*sin(x) *cos(x) - 235*sin(x) *cos(x) 

16 2 16 14 16 

+ 75*sin(x) *cos(x) - 10*sin(x) + 10*sin(x) *cos(x) 

14 14 14 12 

- 100*sin(x) *cos(x) + 400*sin(x) *cos(x) 

14 10 14 8 

- 860*sin(x) *cos(x) + 1100*sin(x) *cos(x) 

14 6 14 4 

- 860*sin(x) *cos(x) + 400*sin(x) *cos(x) 

14 2 14 12 18 

- 100*sin(x) *cos(x) + 10*sin(x) + 5*sin(x) *cos(x) 

12 16 12 14 

- 75*sin(x) *cos(x) + 400*sin(x) *cos(x) 

12 12 12 10 

- 1100*sin(x) *cos(x) + 1780*sin(x) *cos(x) 

12 8 12 6 

- 1780*sin(x) *cos(x) + 1100*sin(x) *cos(x) 

12 4 12 2 12 

- 400*sin(x) *cos(x) + 75*sin(x) *cos(x) - 5*sin(x) 



10 20 10 18 

+ sin(x) *cos(x) - 30*sin(x) *cos(x) 

10 16 10 14 

+ 235*sin(x) *cos(x) - 860*sin(x) *cos(x) 

10 12 10 10 

+ 1780*sin(x) *cos(x) - 2252*sin(x) *cos(x) 

10 8 10 6 

+ 1780*sin(x) *cos(x) - 860*sin(x) *cos(x) 

10 4 10 2 10 

+ 235*sin(x) *cos(x) - 30*sin(x) *cos(x) + sin(x) 

8 20 8 18 

- 5*sin(x) *cos(x) + 75*sin(x) *cos(x) 

8 16 8 14 

- 400*sin(x) *cos(x) + 1100*sin(x) *cos(x) 

8 12 8 10 

- 1780*sin(x) *cos(x) + 1780*sin(x) *cos(x) 

8 8 8 6 

- 1100*sin(x) *cos(x) + 400*sin(x) *cos(x) 

8 4 8 2 

- 75*sin(x) *cos(x) + 5*sin(x) *cos(x) 

6 20 6 18 

+ 10*sin(x) *cos(x) - 100*sin(x) *cos(x) 

6 16 6 14 

+ 400*sin(x) *cos(x) - 860*sin(x) *cos(x) 

6 12 6 10 

+ 1100*sin(x) *cos(x) - 860*sin(x) *cos(x) 

6 8 6 6 

+ 400*sin(x) *cos(x) - 100*sin(x) *cos(x) 

6 4 4 20 

+ 10*sin(x) *cos(x) - 10*sin(x) *cos(x) 

4 18 4 16 

+ 75*sin(x) *cos(x) - 235*sin(x) *cos(x) 

4 14 4 12 

+ 400*sin(x) *cos(x) - 400*sin(x) *cos(x) 

4 10 4 8 

+ 235*sin(x) *cos(x) - 75*sin(x) *cos(x) 

4 6 2 20 

+ 10*sin(x) *cos(x) + 5*sin(x) *cos(x) 



2 18 2 16 

- 30*sin(x) *cos(x) + 75*sin(x) *cos(x) 

2 14 2 12 

- 100*sin(x) *cos(x) + 75*sin(x) *cos(x) 

2 10 2 8 20 

- 30*sin(x) *cos(x) + 5*sin(x) *cos(x) - cos(x) 

18 16 14 12 

+ 5*cos(x) - 10*cos(x) + 10*cos(x) - 5*cos(x) 

10 
+ cos(x) 



compact (com, cos x~2+sin x~2=l); 

10 10 

sin(x) *cos(x) 



Grobner bases 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter 4.1.2) Derive (see chapter ??) Macsyma (see 
chapter 4.3.2) Maple (see chapter 4.4.2) Mathematica (see chapter 4.5.2) 

• load the module groebner for dealing with Grobner bases 

load_package groebner; 

polys := 45*p + 35*s - 165*b - 36, 35*p + 40*z + 25*t - 27*s, 15*w + 25*p*s + 30*z - 18*t - 
165*b**2, - 9*w + 15*p*t + 20*z*s, w*p + 2*z*t - ll*b**3, 99*w - ll*s*b + 3*b**2, b**2 + 
33/50*b + 2673/10000$ 

vars := w,p,z,t,s,b$ 

groebner (polys, vars); 

{60000*w + 9500*b + 3969, 

1800*p - 3100*b - 1377, 

18000*z + 24500*b + 10287, 

750*t - 1850*b + 81, 

200*s - 500*b - 9, 

2 
10000*b + 6600*b + 2673} 



• solving a system of polynomial equations by Grobner bases 
groesolve(polys,vars); 



• 



- 148*sqrt(ll)*i - 461 
{{t= , 

500 

190*sqrt(ll)*i - 139 

w= , 

10000 

490*sqrt(ll)*i - 367 
z= , 

3000 

- 62*sqrt(ll)*i + 59 
p= 

300 

3*( - 5*sqrt(ll)*i - 13) 
s= , 

50 

3*( - 4*sqrt(ll)*i - 11) 
b= }> 

100 

148*sqrt(ll)*i - 461 
{t= f 

500 

- 190*sqrt(ll)*i - 139 
w= , 

10000 

- 490*sqrt(ll)*i - 367 
z= , 

3000 

62*sqrt(ll)*i + 59 
p= 

300 

3*(5*sqrt(ll)*i - 13) 

s= , 

50 

3*(4*sqrt(ll)*i - 11) 
b= }} 

100 



(see chapter4.6.4) 



4.6.3 Rational functions 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter 4.1.3) Derive (see chapter 4.2.3) Macsyma (see 
chapter 4.3.3) Maple (see chapter 4.4.3) Mathematica (see chapter 4.5.3) 

• integration 



int((3*a*b~2-5*a~2*b)/(a~4-2),a); 

1/4 a 1/4 1/4 

(sqrt(2)*b*( - 10*2 *atan( ) + 5*2 *log(2 + a) 

1/4 
2 

1/4 1/4 1/4 

- 5*2 *log( - 2 + a) + 3*log(2 + a)*b 

1/4 2 

+ 3*log( - 2 + a)*b - 3*log(sqrt(2) + a )*b))/8 



• verification by differentiation 
df(ws,a); 

a*b*( - 5*a + 3*b) 



4 
a - 2 



• partial fraction decomposition 

pf((10x~2-llx-6)/(x~3-x~2-2x),x); 

5 2 3 

i , ,— } 

x + 1 x - 2 x 



4.6.4 Solving equations 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter 4.1.4) Derive (see chapter 4.2.4) Macsyma (see 
chapter 4.3.4) Maple (see chapter 4.4.4) Mathematica (see chapter 4.5.4) 

Linear systems 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter 4.1.4) Derive (see chapter 4.2.4) Macsyma (see 
chapter 4.3.4) Maple (see chapter 4.4.4) Mathematica (see chapter 4.5.4) 

solve(2*xl+x2+3*x3-9,xl-2*x2+x3+2,3*xl+2*x2+2*x3-7, xl,x2,x3); 

{{xl=-l,x2=2,x3=3}} 



Nonlinear equations 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter 4.1.4) Derive (see chapter 4.2.4) Macsyma (see 
chapter 4.3.4) Maple (see chapter 4.4.4) Mathematica (see chapter 4.5.4) 

• using decomposition to solve high degree polynomials 

solve(x**8-8*x**7+34*x**6-92*x**5+175*x**4-236*x**3+226*x**2-140*x+46); 



Unknown : x 

sqrt( - sqrt( - 4*sqrt(3) - 3) - 3)*sqrt(2) + 2 
{x= 

2 

- sqrt( - sqrt( - 4*sqrt(3) - 3) - 3)*sqrt(2) + 2 
x= 

2 

sqrt( - sqrt(4*sqrt(3) - 3) - 3)*sqrt(2) + 2 

x= , 

2 

- sqrt( - sqrt(4*sqrt(3) - 3) - 3)*sqrt(2) + 2 

x= , 

2 

sqrt(sqrt( - 4*sqrt(3) - 3) - 3)*sqrt(2) + 2 

x= , 

2 

- sqrt(sqrt( - 4*sqrt(3) - 3) - 3)*sqrt(2) + 2 

x= , 

2 

sqrt(sqrt(4*sqrt(3) - 3) - 3)*sqrt(2) + 2 

x= , 

2 

- sqrt(sqrt(4*sqrt(3) - 3) - 3)*sqrt(2) + 2 



} 

2 



• multiple use of inversion functions 

solve(log(acos(asin(x**(2/3)-b)-l))+2,x); 

1 3/2 

{x=(sin(cos( ) + 1) + b) 

2 
e 

1 3/2 

x= - (sin(cos( ) + 1) + b) } 

2 
e 



Nonlinear systems 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter 4.1.4) Derive (see chapter ??) Macsyma (see 
chapter ??) Maple (see chapter 4.4.4) Mathematica (see chapter 4.5.4) 



solve polynomial systems by the use of (see chapter2.7) 



solve( ( alpha * cl - beta * cl**2 - gamma*cl*c2 + epsilon*c3, -gamma*cl*c2 + (epsilon+theta)*c3 
-eta *c2, gamma*cl*c2 + eta*c2 - (epsilon+theta) * c3), (c3,c2,cl)); 

{{cl=arbcomplex(12) , 

c2=(cl*( - cl*beta*epsilon - cl*beta*theta + alpha*epsilon 

+ alpha*theta) )/(cl*gamma*theta - epsilon*eta) , 

c3=(cl* 

2 
( - cl *beta*gamma + cl*alpha*gamma - cl*beta*eta + alpha*eta) 

)/(cl*gamma*theta - epsilon*eta)}} 

• the value of the operator arbcomplex is an arbitrary complex number 

• (see chapter4.6.2) 

4.6.5 Analytical operations 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Derive (see chapter 4.2.5) Macsyma (see 
chapter 4.3.5) Maple (see chapter 4.4.5) Mathematica (see chapter 4.5.5) 

Limits 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Derive (see chapter 4.2.5) Macsyma (see 
chapter 4.3.5) Maple (see chapter 4.4.5) Mathematica (see chapter 4.5.5) 

limit (sin (x) /x,x,0) ; 



limit((3*sin(pi*x) - sin(3*pi*x))/x~3,x,0); 

3 

4*pi 

limit((2x+5)/(3x-2),x,infinity); 

2 



Taylor series 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Derive (see chapter 4.2.5) Macsyma (see 
chapter 4.3.5) Maple (see chapter 4.4.5) Mathematica (see chapter 4.5.5) 

• load the module taylor for computing Taylor series 
loacLpackage taylor; 
taylor (e**x, x, 0, 4); 



12 13 14 5 

1 + x + *x + *x + *x + 0(x ) 

2 6 24 



taylor ( e **(x+y), x, 0, 2, y, 0, 2); 

12 3 3 

1 + y + *y + x + y*x + (4 terms) + 0(x ,y ) 

2 



taylorcombine (ws**2); 

2 3 3 

1 + 2*y + 2*y + 2*x + 4*y*x + (4 terms) + 0(x ,y ) 



Summation and Products 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Derive (see chapter 4.2.5) Macsyma (see 
chapter 4.3.5) Maple (see chapter 4.4.5) Mathematica (see chapter 4.5.5) 

sum(m**2*x**m,m,l,n); 

n22 n2 n2 n n n n 

(x*(x *n *x - 2*x *n *x + x *n - 2*x *n*x + 2*x *n + x *x + x - x 

3 2 

- l))/(x - 3*x + 3*x - 1) 



sum(cos((2*m-l)*pi/(2*n+l)),m,l,r); 

2*r*pi 

sin( ) 

2*n + 1 



Pi 

2*sin( ) 

2*n + 1 



prod(e**(sin(m*x)),m,l,n); 

cos (x/2) / (2*sin(x/2) ) 

e 

cos((2*n*x + x)/2)/(2*sin(x/2)) 
e 



for all n,m such that fixp m let factorial(n+m)=if m l then factorial(n+m-l)*(n+m) else facto- 
rial(n+m+l)/(n+m+l); 

sum(m*2**m/factorial(m+2),m,l,n); 



n 2 

2*2 + factorial (n)*n + 3*f actor ial(n)*n + 2*f actor ial(n) 

2 
factorial (n)*(n + 3*n + 2) 



Integration 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Derive (see chapter 4.2.5) Macsyma (see 
chapter 4.3.5) Maple (see chapter 4.4.5) Mathematica (see chapter 4.5.5) 

int(x**2*(a+b*x)**p,x); 

p3 2 222 22323 

((a + b*x) *(2*a - 2*a *b*p*x + a*b *p *x + a*b *p*x + b *p *x 

3 3 3 3 3 3 2 

+ 3*b *p*x + 2*b *x ))/(b *(p + 6*p + ll*p + 6)) 



int(x**2*log(x**2+a**2),x); 

x 3 2 2 3 2 3 

- 6*atan( )*a + 3*log(a + x )*x + 6*a *x - 2*x 

a 



int(x*d**x*sin x,x); 

x 3 2 

(d *(sin(x)*log(d) *x - sin(x)*log(d) + sin(x)*log(d)*x + sin(x) 

2 4 

- cos(x)*log(d) *x + 2*cos(x)*log(d) - cos(x)*x) )/(log(d) 

2 
+ 2*log(d) + 1) 



int(x*sqrt(a+b*x)**p,x); 

p/2 2 2 2 2 2 

2* (a + b*x) *( - 2*a + a*b*p*x + b *p*x + 2*b *x ) 

2 2 

b *(p + 6*p + 8) 



int(2*x*e**(x**2)*log(x)+e**(x**2)/x+(log(x)-2)/(log(x)**2+x)**2+ ((2/x)*log(x) + (l/x) + l)/(log(x)**2+x), X ; 

2 2 

x 3 x 2 2 

(e *log(x) + e *log(x)*x + log(log(x) + x)*log(x) 



2 2 

+ log(log(x) + x)*x - log(x))/(log(x) + x) 



Ordinary differential equations 



In Reduce 

For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Derive (see chapter ??) Macsyma (see 
chapter 4.3.5) Maple (see chapter 4.4.5) Mathematica (see chapter 4.5.5) 

• load the module odesolve for the analytical solution of ordinary differential equations 
load_package odesolve; 
depend y,x; 
odesolve (df(y,x) + y * sin x/cos x - 1/cos x,y,x); 

{y=arbconst (l)*cos(x) + sin(x)} 



• Bernoulli equation 

odesolve(x*(l-x**2)*df(y,x) + (2*x**2 -l)*y - x**3*y**3,y,x); 

5 
1 5*arbconst(2) + 2*x 
{ = } 

2 4 2 

y 5*x - 5*x 



solve(ws,y); 



2 
sqrt(x - l)*sqrt (5)*x 
{y= f 

5 
sqrt(5*arbconst(2) + 2*x ) 

2 

- sqrt(x - l)*sqrt (5)*x 
y= } 

5 
sqrt(5*arbconst(2) + 2*x ) 



odesolve(df(y,x,2)+4*df(y,x)+4*y-x*exp(x),y,x); 

3*x 3*x 
3*e *x - 2*e + 27*arbconst (3)*x + 27*arbconst (3) 
{y= } 

2*x 
27*e 



Substitutions - pattern matching 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter 4.1.5) Derive (see chapter ??) Macsyma (see 
chapter 4.3.5) Maple (see chapter ??) Mathematica (see chapter ??) 

factor cos, sin; 

(al*cos(wt) + a3*cos(3*wt) + bl*sin(wt) + b3*sin(3*wt))**3 where cos(~x)*cos(~y) =1 (cos(x+y)+cos(x- 
y))/2, cos(~x)*sin(~y) = L (sin(x+y)-sin(x-y))/2, sin(~x)*sin(~y) = L (cos(x-y)-cos(x+y))/2, cos(~x)**2 
=i (l+cos(2*x))/2, sin(~x)**2 =i (l-cos(2*x))/2; 

2 2 

(cos(9*wt)*a3*(a3 - 3*b3 ) 

2 2 

+ 3*cos(7*wt)*(al*a3 - al*b3 - 2*a3*bl*b3) + 3*cos(5*wt)* 

2 2 2 2 

(al *a3 + al*a3 - 2*al*bl*b3 - al*b3 - a3*bl + 2*a3*bl*b3) + 

3 2 2 3 2 2 

cos(3*wt)*(al + 6*al *a3 - 3*al*bl + 3*a3 + 6*a3*bl + 3*a3*b3 ) 

+ 3*cos(wt)* 

3 2 2 2 2 2 
(al + al *a3 + 2*al*a3 + al*bl + 2*al*bl*b3 + 2*al*b3 - a3*bl ) 

2 2 

+ sin(9*wt)*b3*(3*a3 - b3 ) 

2 2 

+ 3*sin(7*wt)*(2*al*a3*b3 + a3 *bl - bl*b3 ) + 3*sin(5*wt)* 

2 2 2 2 

(al *b3 + 2*al*a3*bl + 2*al*a3*b3 - a3 *bl - bl *b3 + bl*b3 ) + 

2 2 2 3 2 3 

sin(3*wt)*(3*al *bl + 6*al *b3 + 3*a3 *b3 - bl + 6*bl *b3 + 3*b3 ) 

+ 3*sin(wt)* 

2 2 2 3 2 2 

(al *bl + al *b3 - 2*al*a3*bl + 2*a3 *bl + bl - bl *b3 + 2*bl*b3 ) 

)/4 



operator integrate; 

linear integrate; 

let integrate(~x**~p,x) =l x**(p+l)/(p+l) when df(p,x)=0, integrate(~x,x) =l x**2/2, inte- 
grate(l,~x) =l x$ 

integrate(a"2*b+a"b+3*a-5,a); 

b 2 2 2 

a*(6*a + 2*a *b + 2*a *b + 9*a*b + 9*a - 30*b - 30) 

6*(b + 1) 



integrate (a" (a+1), a); 



integrate (a *a,a) 



4.6.6 Matrices 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter 4.1.6) Derive (see chapter 4.2.6) Macsyma (see 
chapter 4.3.6) Maple (see chapter 4.4.6) Mathematica (see chapter 4.5.6) 

matrix xx,yy; 

let xx= mat((all,al2),(a21,a22)), yy= mat((yl),(y2)); 

xx ; 

[all al2] 
[ ] 

[a21 a22] 



det xx; 

all*a22 - al2*a21 



zz:= xx**(-l)*yy; 





[ - al2*y2 + a22*yl ] 

r i 


zz : 


L J 
[ all*a22 - al2*a21 ] 

= [ ] 
[ all*y2 - a21*yl ] 
r i 




I J 
[ all*a22 - al2*a21 ] 



l/xx**2; 



mat ( ( - 



al2*a21 + a22 



2 2 2 2 

all *a22 - 2*all*al2*a21*a22 + al2 *a21 



al2*(all + a22) 



-), 



2 2 2 2 

all *a22 - 2*all*al2*a21*a22 + al2 *a21 



a21*(all + a22) 



2 2 2 2 

all *a22 - 2*all*al2*a21*a22 + al2 *a21 



2 

all + al2*a21 
)) 

2 2 2 2 

all *a22 - 2*all*al2*a21*a22 + al2 *a21 



mateigen(mat((2,-l,l), (0,1,1), (-l,l,l)),et); 

{{et - 1,2, 

[arbcomplex ( 13) ] 
[ ] 

[arbcomplex ( 13) ] 
[ ] 

[ ] 

}, 
{et - 2,1, 

[ ] 
[ ] 

[arbcomplex ( 14) ] 
[ " ] 
[arbcomplex ( 14) ] 



}} 



4.6.7 Code generation 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter ??) Derive (see chapter ??) Macsyma (see 
chapter 4.3.7) Maple (see chapter ??) Mathematica (see chapter ??) 

• load the module gentran for the generation of numerical programs 

load_package gentran; 

gentranout gentst$ 

matrix m(3,3)$ 

m(l,l) := 18*cos(q3)*cos(q2)*m30*p**2 - 9*sin(q3)**2*p**2*m30 - sin(q3)**2*j30y + sin(q3)**2*j30z 
+ p**2*ml0 + 18*p**2*m30 + jlOy + j30y$ 

m(2,l) := m(l,2) := 9*cos(q3)*cos(q2)*m30*p**2 - sin(q3)**2*j30y + sin(q3)**2*j30z - 9*sin(q3)**2*m30*p !| 
+ j30y + 9*m30*p**2$ 



m 



(3,1) 



m(2,2) 
m(3,2) 



m 



= m(l,3) := -9*sin(q3)*sin(q2)*m30*p**2$ 



= -sin(q3)**2*j30y + sin(q3)**2*j30z - 9*sin(q3)**2* m30*p**2 + j30y + 9*m30*p**2$ 
= m(2,3) := 0$ 
(3,3) := 9*m30*p**2 + j30x$ 
gentranlang!* := 'fortran$ 
fortlinelen!* := 72$ 

gentran literal "c", cr!*, "c", tab!*, "*** compute values for matrix m ***" 5 cr!*, "c", cr!*$ 
for j: = l:3 do for k:=j:3 do gentran m(j,k) ::=: m(j,k)$ 



gentran literal "c", cr!*, "c", tab!*, "*** compute values for inverse matrix ***" 5 cr!*, "c", cr!*$ 

share var$ 

for j: = l:3 do for k:=j:3 do if m(j,k) neq then JJ var := tempvar nil; markvar var; m(j,k) : = 
var; m(kj) := var; gentran eval(var) := m(eval(j),eval(k)) H$ 

• matrix m contains 
m := m; 

[tO tl t2] 

[ ] 

m := [tl t3 ] 

[ ] 

[t2 t4] 



matrix mxinv(3,3)$ 

mxinv := m**(-l)$ 

for j: = l:3 do for k:=j:3 do gentran mxinv(j,k) ::=: mxinv(j,k)$ 

gentran for j: = l:3 do for k:=j + l:3 do jj m(kj) := m(j,k); mxinv(kj) := mxinv(j,k) H$ 

gentranshut gentst; 

Generated FORTRAN program 

Generated program 



c 

c *** compute values for matrix m *** 

c 

m(l,l)=-(9.0*sin(real(q3))**2*p**2*m30)-(sin(real(q3))**2*j30y)+ 
. sin (real (q3) )**2*j30z+18.0*cos (real (q3))*cos (real (q2) )*p**2*m30+ 
. 18.0*p**2*m30+p**2*ml0+j30y+jl0y 

m(l,2)=-(9.0*sin(real(q3))**2*p**2*m30)-(sin(real(q3))**2*j30y)+ 
. sin (real (q3) ) **2* j 30z+9 . 0*cos (real (q3) ) *cos (real (q2) ) *p**2*m30+ 
. 9.0*p**2*m30+j30y 

m(l,3)=-(9.0*sin(real(q3))*sin(real(q2))*p**2*m30) 
m(2,2)=-(9.0*sin(real(q3))**2*p**2*m30)-(sin(real(q3))**2*j30y)+ 
. sin (real (q3) )**2*j30z+9.0*p**2*m30+j30y 
m(2,3)=0.0 
m(3 , 3) =9 . 0*p**2*m30+ j 30x 

c 

c *** compute values for inverse matrix *** 

c 

tO=m(l,l) 

tl=m(l,2) 

t2=m(l,3) 

t3=m(2,2) 

t4=m(3,3) 

mxinv (1 , l)=-(t3*t4)/(tl**2*t4+t2**2*t3-(t0*t3*t4) ) 

mxinv(l,2)=(tl*t4)/(tl**2*t4+t2**2*t3-(t0*t3*t4)) 

mxinv(l,3)=(t2*t3)/(tl**2*t4+t2**2*t3-(t0*t3*t4)) 

mxinv (2 , 2) = (t2**2- (t0*t4) ) / (t I**2*t4+t2**2*t3- (t0*t3*t4) ) 

mxinv(2,3)=-(tl*t2)/(tl**2*t4+t2**2*t3-(t0*t3*t4)) 

mxinv (3 , 3) = (t 1**2- (t0*t3) ) / (t I**2*t4+t2**2*t3- (t0*t3*t4) ) 



do 25001 j=l,3 

do 25002 k=j+l,3 
m(k,j)=m(j ,k) 
mxinv(k, j )=mxinv(j ,k) 
25002 continue 

25001 continue 



4.6.8 Graphics 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter 4.1.7) Derive (see chapter 4.2.7) Macsyma (see 
chapter 4.3.8) Maple (see chapter 4.4.7) Mathematica (see chapter 4.5.7) 

2D Graphics 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter 4.1.7) Derive (see chapter 4.2.7) Macsyma (see 
chapter 4.3.8) Maple (see chapter 4.4.7) Mathematica (see chapter 4.5.7) 

• function of one parameter 
on demo,rounded,numval; 

plot(sin(e~x),x=(0 .. pi), xlabel="Time", ylabel=" Signal"); 

• graph of several functions (Bessel functions J(n,x),n = 0,2,5) 

load_package specfn; 

plot(besselj(0,x),besselj(2,x),besselj(5,x),x=(0 .. 10), title=" Bessel Functions BesselJ(n,x)", yla- 
bel="Value"); 

3D Graphics 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter 4.1.7) Derive (see chapter 4.2.7) Macsyma (see 
chapter 4.3.8) Maple (see chapter 4.4.7) Mathematica (see chapter 4.5.7) 

• graph of a function with 2 parameters 
plot_xmesh:=plot_ymesh:=35; 
plot(sin(pi*sinx+y),x=(-3 .. 3),y=(-3 .. 3),hidden3d,view="30,40"); 

• graph of another function with 2 parameters 
plot_xmesh:=plot_ymesh:=50; 

plot(tan(x*y),x=(-2*pi/3.. 2*pi/3), y=(-2*pi/3 •• 2*pi/3), z=(-5 .. 5), hidden3d,view=" 30,30"); 

Parametric plots 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter 4.1.7) Derive (see chapter 4.2.7) Macsyma (see 
chapter 4.3.8) Maple (see chapter 4.4.7) Mathematica (see chapter 4.5.7) 

• surface defined parametrically 

dd:=pi/10; 

w:=for u:=0 step dd until 2*pi collect for v:=-pi/2 step dd until pi/2 collect sin v, sin(2*v)*sin 
u, sin(2*v)*cos u$ 

plot w; 



Contour maps 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter ??) Derive (see chapter ??) Macsyma (see 
chapter 4.3.8) Maple (see chapter 4.4.7) Mathematica (see chapter 4.5.7) 

• contour map of a function with 2 parameters 
plot_xmesh:=plot_ymesh:=40; 

plot(sin(x+cos y), x=(0 .. 3*pi/2), y=(0 .. 3*pi/2), contour ,nosurface,view=" 0,0"); 

• contour map of another function with 2 parameters 
plot_xmesh:=plot_ymesh:=80; 

plot(tan(x*y), x=(-pi .. pi), y=(-pi .. pi), z=(-3 .. 3), contour ,nosurface,view=" 0,0"); 

• contour map of yet another function with 2 parameters 
plot_xmesh:=plot_ymesh:=30; 

plot(e~(-sqrt(x~2+y~2))*cos(atan(x/y)), x=(-l .. 1), y=(-l .. 1), contour, nosurface,view="0,0"); 



4.6.9 Graphical presentation of formulas 

In Reduce 

For comparison with other CAS choose from: Axiom (see chapter ??) Derive (see chapter ??) Macsyma (see 
chapter 4.3.9) Maple (see chapter 4.4.8) Mathematica (see chapter 4.5.8) 



Chapter 5 

Applications of computer algebra 



5.1 Classical application areas 

• celestial mechanics 

— calculation of orbits, gravitational fields 

— Fourier, Poisson series 



E 



Pi sin ^Qj^j +<2jCOS 5^/3y$j 



— Delaunay's (1867) calculation of the orbit of the Moon took him 20 years (including physical effects like 
non-symmentry of the earth and influence of the Sun); recalculation on a small computer (1980) took 
20 hours of CPU time; the check found only a few mistakes in his hand calculations in some of the high 
order terms 

— systems TRIGMAN (1970), CAMAL (1975) 

• general theory of relativity 

— calculations with various metrics 

— systems CAMAL (1975), SHEEP (1977), general purpose systems 

• quantum electro-dynamics 

• high energy physics 

— interaction of particles, Dirac matrices, Feynman diagrams, calculation of integrals 

— systems REDUCE (1968), SCHOONSHIP (1971) 

5.2 Other application areas 

• physics 

— plasma physics, physics of fluids 

— electron optics, non-linear optics 

— molecular physics 

— electronics 

— mechanics 

• mathematics 
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— number theory 

— group theory 

— computational geometry 

— numerical analysis 

• chemistry 

• biology 

• robotics 

• economy 

5.3 Case study 1. - Perturbation methods 

5.3.1 Celestial mechanics - nonlinear algebraic equations 

• from J. P. Fitch, International Conference on Computer Algebra and its Applications in Theoretical Physics, 
p. 262-275, Dubna, (1985). 

• solution of a series in a small parameter 

y 2 = i - e 

• zeroth order solution yo = ±1, n-th order solution y n 

y = y n + v + 0(s n+2 ) 

r] = e n+ H 

(y n +r, + 0(e n+2 )) 2 = l-e 

y 2 n +r, 2 +2r,y n +0(e n+2 ) = l-e 

y 2 n + 2r 1 y + O(e n+2 ) = l-e 

r,= \(l-e-y 2 n )+0(e n+2 ) 

• resulting iteration scheme 



Vo = 1 

y n +i = y n + ^-e-yl) 

• another example: the Kepler equation 

E = u + e sin E 

• zeroth order solution E = u,n-th order solution E n 

E = E n +r] + 0(e n+2 ) 

E n +n = u + e sin(£„ + t}) + 0{e n+2 ) 

E n + r] = u + e(sin E n cos r] + cos E n sin r]) + 0{e n+2 ) 

E n +rj = u + £ sin E n + 0(s n+2 ) 

E n =u + A n 



A n+1 = e sin(u + A n ) 
A„_|_i = e(sin u cos A n + cos u sin A n ) 

Ao = 
A n+1 = e sin u f 1 - -^ + -^ - . . 

+£ COS U f A„ - -^ + -^ 

• final solution in the form of Fourier series 

n 

A n = ^P i (e)sin(ju) 

• Fourier series are used by specialized computer algebra systems such as CAMAL and TRIGMAN 

5.3.2 Mechanics - nonlinear ordinary differential equations 

• from D.M. Klimov, V.M. Rudenko, Computer Algebra Methods in Mechanics, Nauka, (1989). 

• nonlinear periodic movement, nonlinear ordinary differential equation 

U" +ulU =ef(U',U,e) 

• new time variable r = cot, to is an unknown frequency 

to 2 U" +ulU = ef(uU',U,e) 

• Poincare perturbation method — expansion into a power series 



• example: the van der Pol equation 



u = j2 u ^y 

00 
i=0 



U" + U -e U'(l - U 2 ) =0 
to 2 U" + U -etoU'(l-U 2 ) =0 



• after substituting power series in 



e° 



e 1 



e 2 



cu 2 U^ + U = 



lo 2 U[' + U 1 = -2U00J1UH + loqU^I - U 2 ) 
WqUL' + U 2 = -(2lo lo 2 + u\)UZ - 2lu coiU[' 
+wiE^(1 - U 2 ) - ZuoUtfJiUo + u U[(l - U 2 ) 



• initial condition 

• zeroth order solution 

• first order solution 



U" + Ui= 2to 1 A cos t- A sin r(l - A 2 Q cos r) 



• after Fourier expansion 



U" + Ui = 2u; 1 A cos r + A sin r (-^ - 1 ) 



A 3 
+ -^sin(3r) 
4 



• no secular terms in the solution require zero coefficients of sinand coson the right hand side 

2u 1 A = , A (^ - l) = 



• solution 



Uq = 2 cos r, loq = 1 

C/i = Ai cos r + Bi sin r — - sin(3r) , lo\ = 



• using the initial condition 

• second order equation 



first order solution 



• second order solution 



1 2 

U 2 ' +U 2 = 2Ai sin r + (4w 2 + -) cos r - - cos(3r) 

5 
+3^4i sin(3r) + - cos(5r) 

16 



3 1 . , 

U\ = - sin t sin(3r) , to\ = 

4 4 



3 5 

C/2 = ^2 cos r + B2 sin r + — cos(3r) — — cos(5r) 



• etc. for higher orders 



5.3.3 Quantum mechanics - eigenvalue problem 

• from T.C. Scott, R.A. Moore, M.B. Monagan, G.J. Fee, E.R. Vrscay, J. of Comp. Phys., Vol. 87, No. 2, 
p.366-395, (1990). 

• Rayleigh-Schrodinger perturbation theory 

• Hamiltonian, solvable part i^and perturbative part Hi 

H = H + XHi 

• eigenvalues (energies) and eigenfunctions are expanded into power series in A 



e = J2 XPR p 

oo 

$ = ^A p $ p 

p=0 



eigenvalue problem 



hierarchy of equations 



H$ = E$ 



X° 
A 1 
A 2 



H $° = E $° 



Ho® 1 + #i$° = Eo® 1 + £i$° 

H $ 2 +H 1 ^> 1 = E $ 2 +E 1 ^> 1 + E 2 $° 



• normalizations 



X j : H & + H-L&- 1 = ^ Ei&'i 



i=0 



< $|$ >= 1 

A°: <$°|$°>=1 

A 1 : < $°|$ 1 > + < $ 1 |$° >=0 



X j : J2 <$ i \$ j ~ i >=1 

8=0 



• separable problems, second order ordinary differential equations (ODEs), inhomogenous ODE 

»w = f(t) 



^ + p «i7 + « ( " 



• one solution j/i(i)of the homogeneous equation and another linearly independent solution y2(t)'m which the 
degree decreases 



1/2 (*) = 2/i (*) / — ^dt 

W(t) = W(y u y 2 )=exp(- J P{t)dt\ 



• particular solution of the inhomogenous equation, variation of parameters 



m(t)f(t)^ „^_ /" 1/1 (*)/(*) dt 



2/p = u{t)yi{t) + v(t)y 2 (t) 

u() = ~J ^w~ dt ' v{t) = J W 

• the general solution of the inhomogenous equation is 

y = u p + C 1 y 1 + C 2 y 2 

• one particle Dirac equation 

„2 



• 



Hamiltonian 



for hydrogen-like atoms 



(c<3 • p + /3m c + IV)$ = E$ 



H = H + Ai?i 
#o = c<3-p + /3m c 2 + -AI + P) V 

Hi = ^{I-(3)V, X = a 2 



jp 2 y2 Rn L , ZW 

E = m c - Z —5- 1 + 



n 2 V 4n 2 



• two particle Dirac equation 

• necessary algebraic operations 



— indefinite integrals 

— definite integrals 

— pattern matching, zero recognition 

— manipulation and simplification of large sums 

5.4 Case study 2. - General theory of relativity 

5.4.1 Basic notions 

• Riemann geometry, metric tensor <?,j, covariant components gij 

ds 2 = y^gjjdx z dx-' 

• contravariant components g^oi metric tensor 

i 

• Christoffel symbols of the first kind 

p, 1 (dg ik , d g a d g it 

i-ikl — 7i \ ~~^ — T + 



2 V d x l d x k d x k 



• Christoffel symbols of the second kind 



-kl 

3 



• Riemann curvature tensor 



• Ricci tensor 



d r*. d r*. 

pi 3 1 jk , \ s -pi pm \ s pi pm 

jM ~ 'dx k ~ ~ ~dx r 2^ L ™k L 3l ~ Z^ L ™l L 3k 



k 

• for a general metric in 4D, there are 9 990 diagonal and 13 280 off-diagonal terms 

• Ricci scalar 

R = Y.a ij Ri3 

i,3 

• Einstein vacuum equations 

itij — ^-Qij 

• system of second order nonlinear partial differential equations with cosmological constant A 

5.4.2 Examples 

• diagonal metrics, e.g., spherically symmetric metrics 

ds 2 = -exp(p(r))dr 2 -r 2 d@ 2 

-r 2 sin 2 9de 2 + exp(q(r))d t 2 

• non-diagonal metrics, e.g., Bondi metrics 

2 Vexp(2/3) 2 2 ,„ . , 2 

ds = U r exp(27Jdu 

r 

+2 exp(2/3)d udr + 2Ur 2 exp(2j)d ud@ 
-r 2 (exp(2j)d 6 2 + exp(-27) sin 2 6 d $) 

• verification of metric tensor solution g, B. Nielsen and H. Pedersen, SIGSAM Bull., Vol. 22, 1, Issue 83, 
p. 7-11, (1988) 

• coordinate system (ui,U2,U3,U4) 



ii 2 
ds 2 = ^(apl+a-'pl+PpD+p-'dul 

p\ = cos u\ d iii + sin u\ sinh ui d u^ 
pi = — sin u\ d iii + cos u\ sinh u-i d U3 
ps = d iii + cosh u-i d us 

a = v£ _1 v£ _1 



6 4 



• solution of Einstein vacuum equation with 



Rij — 



5.4.3 Other problems and references 

• finding a solution 

• calculation of properties of any solutions discovered (symmetry) 

• problem of equivalence, classification of geometries 

• equations including mass (energy, momentum) 

• references 

- M.A.H. MacCallum, EUROCAL'87, p.34-43, Springer- Verlag,(1989). 

- M.A.H. MacCallum, Computer Algebra in Physical Research, p. 278-287, World Scientific, (1991). 

5.5 Case study 3. - Collision integrals in plasma physics 

• collision integrals in plasma fluid equations 

5.5.1 Basic notions 

• Fokker-Planck equation for distribution function f s 

dfs, d f s F si d f s (d f s 



dt 8 Xi m s dvi \ dt 

Coulomb collisions with a small angle 



df s \ _ dioifs) + ld 2 (b ij f s ) 



dt J „ d Vi 2 dvi dvj 



• Rosenbluth potentials 



a A v ) = Z^ C st § — d°u 

t J 9 

bij{vl = j2Ds t f^ ){Sijg2 - gigj) ^ 



, - 2 3 



g = u - v , g=\g\ 
macroscopic quantities: density N s , velocity u s , pressure p s , viscous pressure P s , heat flow q s 



1 f 

—m s 

2 S J 



c 2 sC si fs d 3 v 



-> -> 2 -e 

V -U s , C=C. 



• collision terms 

R.u = m s Jf s a t d 3 v 
Rs, = m s J M c si a i + l/2 bii )d* V 

ttsdij = Wis I JsyCsi&j """ Csj&i ^f ^ij^si^i 

+bij - l/3Sijb k k)d 3 v 
Rsu = l/2m s jfMa i+ 2c si c sk a k 
+c S ibkk + 2c sk b ik )d z v 

• common form of all collision terms 

• 13-th moment approximation of the distribution function 

f s (x,v,t) = f s0 (x,v,t)(l + $ s (x,v,t)) 



«u 



fsO ~ 


- Ns c-p( C 'l 


^a'^H "2 J 


a 2 ~- 


2fcT, 




m s 


$ s = 


p s o? s 5p s aj 


fsft = 


= fsofto(l + *s + *t + *s*t) 



5.5.2 Analytical calculation of collision integrals 

• from A. Salat, Plasma Physics, Vol. 17, p. 589-607, (1975). 

• common form of the integrals 



exp ( Jl^l _ (*z*T) i^§« 



• transformation of variables, new form 



// 



exp( _^_fc^£§£W,A, 



• sum of 288 terms up to the 9-th degree in v, one term of the 9-th degree in ogives after substitutions 19 683 
terms, in one part of one collision integral is 259 584 terms (one term is a product appearing in a sum) 

• integration over d 3 Cderived from 

exp(— x 2 )dx = \pK 



/< 



integration over (fig; rotation and transformation into spherical coordinates g,9,(j>(z = cos 9) 



Lrr^M-^)'"- 



cji = gz— + g\Jl - z 2 (e 2i cos <f> + e 3i sin <f>) 
u 

integration over d<f>, products cos<f), sin<^, elimination of e 2 i, e-a 

UiUj 



integration over dg 



g exp 



2j T C3JC3J — 


VlJ 


u 2 


g 2 - 2guz\ 
a 2 ) 


dg = 




a 2 \ n d n 


exp ( 


'u 2 z 2 


2z) du n 





erfc 



(-t) 



• integration over dz 



J —a 



W 
h 
h 



w + exp(w )erfc(— w)dw 



uz u 

— , a= - 
a a 

exp(a 2 )erf (a) — a 
a 2k exp(a 2 )erf (a) - 



■2k+l 



2k + 1 



kh-i 



• final result contains only elementary functions and the error function 

• algebraic program in REDUCE was unable to calculate all the collision integrals and so had to be rewritten 
into a more efficient symbolic program 

• speed up of the symbolic program compared to the algebraic program depends on the size of a collision integral 



Size 


small medium large 


Speed up 


3x 7x ? 



5.6 Case study 4. - Numerical solving of partial differential equations 

5.6.1 References 

• weighted residual method, replacement of the solution by its approximation, minimization of the residual error 

R. Zippel. Symbolic/numeric techniques in modeling and simulation. In B. Donald, D. Kapur, and J. Mundy, 
editors, Symbolic and Numerical Computation in Artificial Intelligence. Academic Press, 1992. 

• finite element methods 

P. S. Wang. FINGER: A symbolic system for automatic generation of numerical programs in finite element 
analysis. J. Symb. Comp., 2(3):305-316, 1986. 

• finite difference methods 

M. C. Wirth. On the Automation of Computational Physics. PhD thesis, Lawrence Livermore National 
Laboratory, Livermore, 1980. UCRL-52996. 



S. Steinberg and P. J. Roache. Symbolic manipulation and computational fluid dynamics. J. Comp. Phys., 
57:251-284, 1985. 

G. O. Cook Jr. ALPAL: a program to generate physics simulation codes from natural descriptions. Int. J. 
Mod. Phys. C, 1(1):1-51, 1990. 

R. Liska and L. Drska. FIDE: A REDUCE package for automation of Finite difference method for solving 
pDE. In S. Watanabe and M. Nagata, editors, Proceedings of ISSAC'90, 169-176, New York, 1990. ACM 
Press, Addison Wesley. 

• numerical code generation 

B. L. Gates. A numerical code generation facility for REDUCE. In B. W. Char, editor, Proceedings of SYMSAC 
'86, 94-99, Waterloo, 1986. ACM. 

N. Sharma and P. Wang. Symbolic derivation and automatic generation of parallel routines for finite element 
analysis. P. Gianni, editor, Proceedings of ISSAC'88, 33-56, Berlin, 1988. Springer- Verlag. 

• IRENA-interface to the numerical library (NAG, REDUCE) 

M. C. Dewar and M. G. Richardson. Reconciling symbolic and numeric computation in a practical setting. In 
A. Miola, editor, Proceedings of DISCO'90, 195-204, Berlin, 1990. Springer- Verlag. 



5.7 Survey articles on applications 



• W.S. Brown and A.C. Hearn. Applications of Symbolic Mathematical Computations. Comput. Phys. Com- 
mun. Vol. 17, 207-215, 1979. 

• J. Calmet and J. A. van Hulzen. Computer algebra applications. In B. Buchberger, G.E. Collins and R. Loos, 
editors, Computer Algebra, Symbolic and Algebraic Computation, 245-258, Berlin, 1983. Springer- Verlag. 

• A. M. Cohen, editor. The SCAFI Papers, Proceedings of the 1991 SCAFI, Seminar on Studies in Computer 
Algebra for Industry, CAN, CWI Amsterdam, 1991. 

• H.I. Cohen and J. P. Fitch. Uses Made of Computer Algebra in Physics, J. Symbolic Computation, Vol. 11, 
No. 3, 291-305, 1991. 

• J. P. Fitch. The Application of Symbolic Algebra in Physics - a Case of Creeping Flow. In E. W. Ng, editor, 
EUROSAM'79, 30-41, 1979. Springer- Verlag. 

• J. P. Fitch. Applying Computer Algebra. In International Conference on Computer Algebra and its Applica- 
tions in Theoretical Physics, 262-275, 1985. Dubna. 

• E.W. Ng. Symbolic-Numeric Interface: A Review. In E.W. Ng, editor, EUROSAM'79, 330-345, 1979. 
Springer- Verlag. 



Chapter 6 

Another sources of study 



6.1 Basic references 

• B. Buchberger, G.E. Collins, and R. Loos, editors. Computer Algebra, Symbolic and Algebraic Computation. 
Springer- Verlag, Wien, 1983. first book on computer algebra algorithms 

• J.H. Davenport, Y. Siret, and E. Tournier. Computer Algebra, Systems and Algorithms for Algebraic Com- 
putation, Academic Press, London, 1988. the best introduction to computer algebra algorithms and systems, 
introduction to REDUCE 

• K.O. Geddes, S.R. Czapor and G. Labahn. Algorithms For Computer Algebra. Kluwer Academic Publishers, 
Boston, 1992. the best reference book on computer algebra algorithms 

• D. Harper, C. Wooff, and D. Hodgkinson. A Guide to Computer Algebra Systems. John Wiley & Sons, 
Chichester, 1991. comparision of the systems Derive, MACSYMA, Maple, Mathematica, REDUCE 

6.2 Other references 

• algorithms 

R. Zippel. Effective Polynomial Computation. Kluwer Academic Publishers, Boston, 1993. 

M. Mignotte. Mathematics for Computer Algebra. Springer- Verlag, Berlin, 1992. 

F. Winkler. Polynomial Algorithms in Computer Algebra. Springer- Verlag, Berlin, 1996. ISBN 3-211-82759-5 

• design of symbolic computing systems 

A. Miola, M. Temperini (eds.). Advances in the Design of SC Systems. Springer- Verlag, Berlin, 1996. 

• Derive 

W. Ellis, and E. Lodi. Derive for the Calculus Student: A Tutorial. Brooks/Cole, Pacific Grove, 1991. 
D.C. Arney. Exploring Calculus with DERIVE. Addison-Wesley, 1992. 

• Macsyma 

Richard Pavelle, editor. Applications of Computer Algebra, Kluwer Academic Publishers, 1985. 

Barbara Heller. Macsyma for Statisticians. John Wiley and Sons, 1991. 

Chiang C. Mei. Mathematical Analysis in Engineering. Cambridge University Press, 1994. 

• Maple 

A. Heck. Introduction to Maple. Springer- Verlag, Berlin 1993. 
D. Redfern. Maple Handbook. Springer- Verlag, Berlin 1993. 
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• Mathematica 

R. Maeder. Programming in Mathematica. Addison- Wesley, Redwood City, CA, 1990. 
R.E. Crandall. Mathematica for the Sciences. Addison- Wesley, Redwood City, CA, 1991. 
R. Maeder. The Mathematica Programmer. Academic Press Professional, 1994. 

• REDUCE 

F. Brackx. Computer Algebra with LISP and REDUCE. Kluwer Academic Publishers, Boston, 1991. ISBN: 
0-7923-1441-7, 300 p. 

M.A.H. MacCallum, and F.J. Wright. Algebraic Computing with REDUCE, Clarendon Press, Oxford, 1991. 

G. Rayna. REDUCE Software for Algebraic Computation. Springer- Verlag, Berlin, 1987. 

6.3 Journals 

• Journal of Symbolic Computation , Academic Press , monthly, 2 volumes per year, since 1985, basic journal 
for the theory of computer algebra, editor Bob Caviness, University of Delaware 

• SIGSAM Bulletin, bulletin of the special interest group SIGSAM , ACM Press, quarterly, since 1967, editor 
Robert Corless 

• Applicable Algebra in Engineering, Communication and Computing, Springer International, quarterly, since 
1990, editor J. Calmet, Karlsruhe 

• Macsyma Newsletter, Macsyma Inc. 

• MapleTech, The Maple Technical Newsletter, Birkhauser, Boston, biannually, since 1989, editor T. Scott, 
Waterloo 

• Maple Roots Report, The Newsletter from Waterloo Maple Software 

• MathUser, The Wolfram Research Newsletter for Mathematica Users, biannually, Wolfram Research, Inc., 
mathuser @ wri . com 

• Mathematica Journal, Addison- Wesley, quarterly, since 1990 

6.4 Electronic information sources 

6.4.1 General electronic information sources 

• news group sci. math. symbolic, much on Mathematica and Maple 

• information on the special interest group ACM SIGSAM, Association for Computing Machinery, Special In- 
terest Group on Symbolic and Algebraic Manipulation, 1515 Broadway, New York, NY 10036 

information about ACM 

information about SIGSAM 

• WWW Computer algebra servers 

SymbolicNet Symbolic Mathematical Computation Information Center, Kent State University, includes mailing 
list for announcements 

CAIN Europe Computer Algebra Information Network, CAN, Computer Algebra Netherlands 

RISC , Research Institute for Symbolic Computation, Linz 

• overview of computer algebra systems and packages, including both commercial and public systems, can be 
found at CAIN or at the collection of symbolic software at the University of Berkeley 



6.4.2 Electronic resourses related to particular systems 

Axiom electronic information sources 

• WWW site 

• bibliography 

Derive electronic information sources 

• WWW site 

• email: swh@aloha.com 

Macsyma electronic information sources 

• WWW site 

• for information about Macsyma 
e-mail:info@macsyma.com 

• for service on Macsyma 
e-mail:service@macsyma.com 

• ftp site:fermat. macsyma. com for extra applications packages, patches and demos 

Maple electronic information sources 

• WWW site 

• e-mail 

— information and sale:info@maplesoft.on.ca 

— technical support:support@maplesoft.on.ca 

• Maple Share Library - library of users programs 

— anonymous ftp://ftp.maplesoft.com 
anonymous ftp://ftp.can.nl/pub/maple-ftplib/ 

— e-mail:maple-netlib@can.nl ("send info" in the message body) 

— library contributions: Michael Monagan 

• Maple Users Group 

for information on how to subscribe to the mailing list please send email to majordomo@daisy.uwaterloo.ca 
with the command "info maple-list" in the body of the message. 

Mathematica electronic information sources 

• WWW site 

• e-mail 

— general information, sale 

— general information, sale in Europe 

— customer service 

— user registration 

— technical questions and support , send bugs in Mathematica here 

— technical questions and support in Europe 



— suggestions 

— MathUser Newsletter for users 

• Mathsource 

— library of Mathematica packages, notebooks, technical reports, examples, news and information 

— e-mail:mathsource@wri.com ("help intro" in the message body) 

— ftp 

— support 

REDUCE electronic information sources 

• WWW site at Koln 
WWW site at ZIB Berlin 

• REDUCE secretary 

• REDUCE Network Library - library of user written REDUCE packages 

— WWW site 

— e-mail US:reduce-netlib@rand.org 

— e-mail Europe:elib@elib.zib-berlin.de ("help" or "send index" in the message body) 

• REDUCE Forum - discussion group 

— contributions 
subscribtion 

6.5 References for CAS Comparisons 

Asl94 Helmer Aslaksen, Multiple-valued complex functions and computer algebra, Research Report No. 631, De- 
partment of Mathematics, National University of Singapore, October 1994. 

Ber96 Laurent Bernardin, "A Review of Symbolic Solvers", SIGSAM Bulletin, Volume 30, Number 1, March 1996, 
9-20. 

Fau95 J. C. Faugere, "GB: State of GB + Tutorial, Proceedings of the POSSO Workshop on Software, draft version, 
edited by Jean-Charles Faugere, Joel Marchand and Renaud Rioboo, Paris, March 1-4, 1995, 55-60. 

Gra95 Hans-Gert Grabe, "On Factorized Grobner Bases", Computer Algebra in Science and Engineering, edited 
by Fleischer, Grabmeier, Hehl and Kuchlin, World Scientific Singapore 1995, 77-89. 

Har91 David Harper, Chris Wooff and David Hodgkinson, A Guide to COMPUTER ALGEBRA SYSTEMS, John 
Wiley & Sons, 1991. 

Her94 W. Hereman, "Review of Symbolic Software for the Computation of Lie Symmetries of Differential Equa- 
tions", Euromath Bulletin, Volume 1, Number 2, 1994, 45-79. 

Her95 Willy Hereman, "Visual data analysis: maths made easy", Physics World, Volume 8, Number 4, April 1995, 
49-53. 

Her96 Willy Hereman, "Computer algebra: lightening the load", Physics World, Volume 9, Number 3, March 1996, 
47-52. 

Pos96 Frank Postel and Paul Zimmermann, "A review of the ODE solvers of Axiom, Derive, Macsyma, Maple, 
Mathematica, MuPAD and Reduce" , submitted to the 5th Rhine Workshop on Computer Algebra to be held 
in Saint-Louis, France, April 1-3, 1996. 



Rob93 Nicolas Robidoux, "Does Axiom Solve Systems of O.D.E.'s Like Mathematica?" , LA-UR-93-2235, Los 
Alamos National Laboratory, Los Alamos, New Mexico. 

Sim92 Barry Simon, "Comparative CAS Reviews", Notices of the American Mathematical Society, Volume 39, 
Number 7, September 1992, 700-710. 

Sto91 David R. Stoutemyer, "Crimes and Misdemeanors in the Computer Algebra Trade", Notices of the American 
Mathematical Society, Volume 38, Number 7, September 1991, 778-785. 

Wes94 Michael Wester, "A Review of CAS Mathematical Capabilities", Computer Algebra Nederland Nieuwsbrief, 
Number 13, December 1994, ISSN 1380-1260, 41-48 (newer version of the paper below). 

Wes95 Michael Wester, "A Review of CAS Mathematical Capabilities", Applied Mechanics in the Americas, Vol- 
ume III, edited by Luis A. Godoy, Sergio R. Idelsohn, Patricio A. A. Laura and Dean T. Mook, American 
Academy of Mechanics and Asociacion Argentina de Mecanica Computacional, Santa Fe, Argentina, 1995, 
450-455. 

Zim95 Paul Zimmermann, "Wester's test suite in MuPAD 1.2.2", Computer Algebra Nederland Nieuwsbrief, Num- 
ber 14, April 1995, ISSN 1380-1260, 53-64. 

6.6 Computer Algebra Conferences 

6.6.1 International 

ISSAC 

International Society of Symbolic and Algebraic Computation (ISSAC) Conferences 

ISSAC'89 1989 Portland, Oregon, USA 

ISSAC '90 1990 Tokyo, Japan 

ISSAC'91 1991 Bonn, Germany 

ISSAC '92 1992 University of California, Berkeley, California, USA 

ISSAC'93 July 6-8, 1993 Kiev, Ukraine 

ISSAC '94 August 25-27, 1994 London, United Kingdom 

ISSAC'95 July 10-12, 1995 Concordia University, Montreal, Canada 

ISSAC'96 July 24-26, 1996 ETH, Zurich, Switzerland 

ISSAC'97 July 21-23, 1997 Maui, Hawaii, USA 

ISSAC is the major yearly international computer algebra conference. 

IMACS-ACA 

IMACS Conferences on Applications of Computer Algebra ( ACA ) 
ACA'95 May 16-19, 1995 University of New Mexico, Albuquerque, New Mexico, USA 
ACA'96 July 17-20, 1996 RISC-Linz, Hagenberg, Austria 
ACA'97 July 24-26, 1997 Maui, Hawaii, USA 

The primary goal of these conferences is to promote the interaction of users of computer algebra, in particular, 
scientists, engineers and educators. 

6.6.2 Systems meetings 

Axiom 

(Thanks to Grant Keady.) 

Jacques Calmet will have details. 



International AXIOM Meeting 1996? Karlsruhe, Germany 



Derive 



(Thanks to Bernhard Kutzler.) 



Derive User Group (DUG) Meetings 



1st UK 
2nd UK 
1st German 
2nd German 
1st US 
2nd US 



April 14, 1992 
September 20, 1993 
April 24, 1993 
April 14, 1995 
November 20, 1994 
November 19, 1995 



Nottingham, UK 
Birmingham, UK 
Schweinbach, Germany 
Nuernberg, Germany 
Orlando, Florida, USA 
Houston, Texas, USA 



Derive Conferences 



1st Scandinavian DERIVE Conference 
1st International DERIVE Symposium 
2nd International DERIVE Symposium 
3rd International DERIVE Symposium 
1st International DERIVE Conference 
2nd Int. DERIVE/TI-92 Conference 
DERIVE Days Duesseldorf 
DERIVE Days Leeds 



October 1-2, 1993 
April 26-30, 1992 
September 26-30, 1993 
July 30-August 3, 1995 
July 11-15, 1994 
July 2-6, 1996 
April 19-21, 1995 
April 13-15, 1996 



Kungsbacka, Sweden 
Krems, Austria 
Krems, Austria 
Honolulu, Hawaii, USA 
Plymouth, UK 
Bonn, Germany 
Duesseldorf, Germany 
Leeds, UK 



Macsyma 

(Thanks to Jeff Golden.) 



MACSYMA Users' Conferences (MUCs) 



MUC I July 27-29, 1977 University of California at Berkeley, USA 

Proceedings published as NASA CP-2012 
MUC II June 20-22, 1979 Washington, D.C., USA 

MUC III July 23-25, 1984 General Electric, Schenectady, New York, USA 
Macsyma and PDEase in Undergraduate Education 

August 13, 1996 University of Washington, Seattle, USA 



Maple 

(Thanks to Michael Monagan.) 

The Maple retreats were held from 1982 to 1994 every year in June at Sparrow Lake in Southern Ontario, Canada. 
This was a quiet setting where people would relax in a pleasant away from work atmosphere. The meetings were 
attended by Maple developers, local Maple usrs, and typically 4 invited guests, some of whom gave a presentation. 
The original purpose of the meetings was to have a small group get together for brain storming sessions which would 
be uninterrupted over the course of 3 days. Later, as the meetings became larger, and more popular, and people 
from the Maple company also took part, the Maple retreats became more like a scientific workshop meething with 
prepared presentations from many speakers, and large group discussions. These were useful for finding out what 
people were doing and what needed to be done, but they did not provide a good atmosphere where people would 
work on the design of Maple. The Maple retreat meetings were stopped in 1994. Since then a smaller group of 
typically 6 to 10 people have met for 2 day and 1 day meetings to focus on design issues. 



Mathematica 






(Thanks to Michael Trott.) 




1st 


January 11- 


13, 1990 


2nd 


January 12- 


-15, 1991 


3rd 


May 27-31, 


1992 


4th 


September ! 


2-4, 1992 


5th 


April 21-23 


, 1994 


1st Australian 


July 8-10, ] 


L995 



Mathematica Conferences 



Redwood City, California, USA 

San Francisco, California, USA 

Boston, Massachusetts, USA 

Rotterdam, Netherlands 

Champaign, Illinois, USA 

University of Tasmania, Hobart, Australia 



Mathematica Developer Conferences 



1st May 6-8, 1993 
2nd October 6-8, 1995 



Champaign, Illinois, USA 
Champaign, Illinois, USA 



6.7 Manuals 

• AXIOM 

R.D. Jenks, and R.S. Sutor. AXIOM, the Scientific Computation System. Springer- Verlag, 1992. 

Axiom Release 2.0 Companion Guide, The Numerical Algorithms Group Limited, Oxford, March 1995. 

S.M. Watt, P.A. Broadbery, S.S. Dooley, P. Iglio, S.C. Morrison, J.M. Steinbach, R.S. Sutor, IBM Thomas J. 
Watson Research Center. AXIOM Library Compiler, User Guide, The Numerical Algorithms Group Limited, 
Oxford, November 1994. 

• Derive 

User Manual, DERIVE, A Mathematical Assistant for Your Personal Computer, Version 2. Soft Warehouse, 
Inc., Honululu, 1990. 

Albert Rich, Joan Rich, Theresa Shelby and David Stoutemyer. User Manual DERIVE Version 3: A Mathe- 
matical Assistant for Your Personal Computer, Soft Warehouse, Inc., September 1994. 

• Macsyma 

Macsyma Mathematics and System Reference Manual, 15th edition, Macsyma Inc., 1995 
Macsyma User's Guide: A Tutorial Introduction, Second Edition, Macsyma Inc., 1995. 

• Maple 

— new manuals 

K. M. Heal, M. L. Hansen and K. M. Rickard, Maple V Learning Guide, Springer- Verlag, 1996. 

M. B. Monagan, K. O. Geddes, K. M. Heal, G. Labahn and S. M. Vorkoetter, Maple V Programming 

Guide, Springer- Verlag, 1996. 

Darren Redfern, The Maple Handbook: Maple V Release 4, Springer- Verlag, 1996. 

— older manuals 

BW. Char, K.O. Geddes, G.H. Gonnet, B.L. Leong, M.B. Monagan, and S.M. Watt. First Leaves: A 

Tutorial Introduction to Maple V. Springer- Verlag, New York, 1992. 

BW. Char, K.O. Geddes, G.H. Gonnet, B.L. Leong, M.B. Monagan, and S.M. Watt: Maple V Library 

Reference Manual. Springer- Verlag, New York, 1991. 

B.W. Char, K.O. Geddes, G.H. Gonnet, B.L. Leong, M.B. Monagan, and S.M. Watt: Maple V Language 

Reference Manual. Springer- Verlag, New York, 1991. 

• Mathematica 

S. Wolfram. Mathematica, A system for Doing Mathematics by Computer. Addison- Wesley, Redwood City, 
CA, 1991. 

S. Wolfram. The Mathematica Book, Third Edition, Mathematica Verison 3, Wolfram Media/Cambridge 
University Press, 1996. 

• REDUCE 

A.C. Hearn and J.P. Fitch (ed.), REDUCE User's Manual 3.6, RAND Publication CP78 (Rev. 7/95), Rand, 
Santa Monica, CA, 1995. 

• for other computer algebra systems and packages, including also public systems, consult the overview at CAIN 
Netherlands or the collection of symbolic software at the University of Berkeley 

6.8 Distributors addresses 

• AXIOM: NAG Ltd. Wilkinson House, Jordan Hill Road, OXFORD, 0X2 8DR, United Kingdom 

• Derive: Soft Warehouse, Inc. 3615 Harding Avenue, Suite 505, Honolulu, Hawaii 96816-3735, U.S.A. email: 
swh@aloha.com 



• Macsyma: Macsyma Inc., 20 Academy Street, Arlington, Massachusetts 02174-6436, U.S.A., e-mail:info- 
macsyma@macsyma.co, 

Scientific Software Service, Attention: Mr. Gregory Kapsias, Niddastrasse 108, D-60329 Frankfurt /a Main 1, 
Germany, tel: (49) 69-252255, fax: (49) 69-232464, e-mail:100144.347@compuserve.com, 

• Maple: Waterloo Maple Software, 450 Phillip Street, Waterloo, Ontario, N2L 5J2, Canada, e-mail:info@maplesoft.on.ca 
Waterloo Maplo Software GmbH, Tiergartenstrasse 17, W-6900, Heidelberg, Germany 

Czech Software First CS 1, Jiri Hrebicek, Brno, tel/fax: 05-74 1248, Czech Republic 

• Mathematica: Wolfram Research, Inc. 100 Trade Center Drive, Champaign, IL 61820-7237, U.S.A., e- 
mail:info@wri.com 

Wolfram Research Europe Ltd., Evenlode Court, Main Road, Long Hanborough, Oxon, OX8 2LA, United 
Kingdom , e-mail : info-euro@ wri . com 

ELKAN s.r.o. V tunich 12, 120 00 Praha 2, tel. 02-23 55 473, Czech Republic 

• REDUCE: REDUCE secretary, RAND, 1700 Main Street, P.O. Box 2138, Santa Monica, CA 90407-2138, 
U.S.A., e-mail:reduce@rand.org 

Herbert Melenk, Symbolik, Konrad Zuse Zentrum fur Informationstechnik, ZIB Berlin, Heibronner Str. 10, 
D-10711 Berlin- Wilmersdorf, Germany, e-mail:melenk@cs.zib-berlin.de 

Codemist Ltd. " Alta" , Horsecombe Vale, Combe Down, Bath BA2 5QR, United Kingdom, e-mail:jpff@maths. bath. ac.uk 

Richard Liska, Faculty of Nuclear Sciences and Physical Engineering, Czech Technical University, Brehova 7, 
115 19 Praha 1, Czech Republic, e-mail:liska@siduri. fjfi.cvut.cz 

• for other computer algebra systems and packages, including also public systems, consult the overview at CAIN 
Netherlands or the collection of symbolic software at the University of Berkeley 



